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A  Thinned  High  Frequency  Linear  Antenna  Array 
to  Study  Ionospheric  Structure 


ABSTRACT 

This  study  reports  on  the  design,  modeling  and  performance 
measurements  of  a  high  frequency  (HF)  linear  antenna  array  with  36 
sensors.  The  array  was  designed  to  achieve  a  narrow  azimuthal 
beamwidth  while  maintaining  grating  lobes  10  dB  below  the  main  beam 
during  a  +/-30°  scan  in  azimuth.  The  configuration  chosen  utilizes 
two  active  vertical  monopole  elements  and  two  parasitic  backpoles  to 
form  a  subarray.  The  subarrays,  or  sensors,  are  spaced  at  distances 
greater  than  half  a  wavelength  to  provide  a  large  effective  array 
aperture  while  the  elemental  radiation  pattern,  provided  by  the 
subarrays,  suppresses  the  grating  lobe  as  the  array  is  scanned. 
Radiation  patterns  for  the  array  were  determined  using  three 
independent  techniques;  theoretical  calculation,  computer  modeling 
using  the  Numerical  Electromagnetics  Code  (NEC) ,  and  measurement  of 
the  fielded  antenna.  Results  showed  close  agreement  in  antenna 
performance  among  the  three  methods  of  pattern  determination.  The 
three  step  process  of  theory,  numerical  modeling  and  measurement 
appears  to  be  an  optimum  approach  to  antenna  design. 


2 


Table  of  Contents 


1.  Introduction 

1.1  RADC  High  Frequency  Test  Facility 

2 .  Linear  Array  Antenna 

2.1  Design  objectives 

2.2  Verona  array  description 

3.  Theoretical  Methods  and  Considerations 

3.1  Calculation  of  theoretical  radiation  patterns 

3.1.1  Simple  vertical  monopole  element. 

3.2  Two  Element  Subarray. 

3.2.1  Two  active  elements  phased  broadside 

3.2.2  Active  element  and  backpole  phased  endfire 

3.3  Design  of  Subarray  Dimensions. 

3.4  Calculation  of  the  array  factor. 

3.5  Pattern  multiplication  of  array  factor  and  subarray  pattern. 

3.6  Reducing  grating  lobes. 

4.  Calculations  using  Numerical  Electromagnetics  Code  (NEC) 

4.1  Introduction  to  NEC. 

4.2  Far-field  Patterns  Calculated  with  NEC. 

4.2.1  Vertical  monopole  over  perfect  ground 

4.2.2  Two  active  elements. 

4.2.3  Monopole  with  backpole. 

4.2.4  Full  subarray 

4.2.5  Model  of  the  full  array. 

5.  Field  Measurements 

5.1  Introduction  and  techniques 

5.2  Measurement  of  the  subarray  pattern 

5.3  Measurement  of  the  array  pattern 

6.  Array  Description 

7.  Discussion 


References 


3 


Illustrations : 

1.  Verona  linear  array  coverage 

2.  Subarray  configuration  for  the  Verona  linear  array 

3.  Verona  linear  antenna  array  setup 

4.  Theoretical  elevation  patterns  for  a  vertical  monopole. 

5.  Theoretical  azimuthal  far-field  pattern  for  two 
vertical  monopoles  phased  broad-side. 

6.  Theoretical  pattern  for  one  active  element  and  one 
backpole  phased  end-fire. 

7.  Theoretical  azimuthal  far-field  pattern  for  a  subarray 
with  two  active  elements  and  two  backpoles. 

3.  Array  factor,  subarray  pattern,  and  resultant 

radiation  pattern  from  pattern  multiplication  of  array 
factor  with  the  subarray  pattern. 

9.  Array  factor,  subarray  pattern,  and  resultant  radiation 
pattern  for  linear  array  scanned  20  deg.  off  broadside. 

10.  Elevation  radiation  pattern  for  a  simple  vertical  monopole 
over  perfect  and  finite  ground  planes  calculated  with  the  NEC 
modeling  program. 

11.  Azimuthal  pattern  for  a  two  element  subarray  phased 
broadside  calculated  with  NEC. 

12.  Azimuthal  radiation  pattern  for  a  two  element  subarray 
phased  end-fire  calculated  using  the  NEC  routine. 

13.  NEC  azimuthal  radiation  pattern  for  a  four  element 
subarray  including  two  active  elements  and  two  backpoles. 

14.  NEC  radiation  pattern  for  a  linear  array  of  15  subarrays, 
including  two  active  elements  and  two  backpoles.  Zero  degree 
scan  angle. 

15.  NEC  radiation  pattern  for  a  linear  array  of  15  subarrays, 
with  a  32  dB  cosine  squared  weighting  taper. 

16.  Measured  azimuthal  radiation  pattern  for  the  fielded  four 
element  subarray. 

17.  Measured  36  element  linear  array  radiation  pattern  using  the 
Ava  Tx.  site  as  a  source  with  a  32  dB  cosine  square  weighting 

13.  Schematic  of  one  active  element. 

19.  Subarray  schematic. 


4 


1. 


INTRODUCTION 


1.1  RADC  High  Frequency  Test  Facility 

RADC/EECP  has  assembled  a  high  frequency  (HF)  radar  and 
communication  test  facility  at  the  RADC  Verona  Test  Annex,  Verona 
NY.  The  system  consists  of  a  set  of  36  digital  HF  radio  receivers 
controlled  by  a  DEC  Micro-Vax  II  computer  and  an  extended  aperture 
linear  antenna  array.  The  mission  of  the  facility  is  two-fold. 
First,  investigate  propagation  mechanisms  and  limits  imposed  on  HF 
systems  by  the  ionosphere.  Second,  explore  digital  processing 
techniques  such  as  adaptive  sidelobe  cancellation,  narrowband 
noise  excision,  and  time  domain  noise  excision  that  could  lead  to 
improvement  of  current  and  future  Air  Force  surveillance  and 
communications  systems  which  use  the  HF  frequency  band. 

The  radar  system  at  Verona  converts  the  RF  signals  from  the  36 
antenna  sensors  of  the  linear  array  to  a  digital  format  using  A/D 
converters  in  each  of  36  receivers.  Antenna  calibration, 
beamforming,  and  all  radar  processing  is  carried  out  by  the  system 
computer  during  post  processing.  This  paper  details  only  the 
linear  antenna  array.  Information  concerning  the  receivers, 
operational  methods,  and  data  processing  will  be  presented  in  a 
separate  paper. 

The  purpose  of  equipping  the  test  facility  with  an  extended 
aperture  linear  array  was  to  create  a  high  resolution  probe  to 
study  ionospheric  structure  and  use  this  knowledge  to  improve 
backscatter  rada1-  techniques. 
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2. 


LINEAR  ANTENNA  ARRAY 


2.1  Design  Objectives 

The  objective  of  the  HF  test  facility  at  Verona  is  to  measure 
ionospheric  backscatter  signals  with  high  resolution  and  record 
the  receiver  baseband  digitally.  There  are  many  sources  of 
ionospheric  clutter;  radio  aurora,  f-region  irregularities, 
equatorial  irregularities  and  others.  All  these  mechanisms 
contribute  interference  that  is  known  to  degrade  HF  propagation 
performance.  Measurements  using  the  linear  antenna  array  will 
help  provide  the  data  needed  to  evaluate  the  effectiveness  of 
digital  processing  techniques  to  reduce  interference  caused  by  the 
irregularities  at  the  high  latitudes.  The  backscatter  data  will 
also  be  used  to  characterize  or  map  ionospheric  phenomena  in  the 
high  latitudes.  This  mapping  will  help  determine  the  spatial  and 
temporal  extent  of  ionospheric  irregularities. 

The  frequency  range  of  the  receiver  system  is  5-30  MHz  which 
encompasses  most  of  the  HF  band.  However,  high  intensity 
phenomena  are  usually  present  during  nighttime  operation  when  the 
ionosphere  dictates  that  the  most  effective  operating  frequencies 
are  at  the  low  end  of  the  HF  band.  Therefore,  6-12  MHz  was  used 
as  the  design  operating  band  for  the  linear  array.  While  the 
array  will  certainly  work  to  some  extent  over  a  much  larger 
frequency  band,  the  array  was  tuned  to  perform  best  in  the  6-12 
MHz  region.  With  this  criterion  in  mind,  the  36  element  linear 
array  was  configured  to  obtain  maximum  resolution  in  azimuth  while 
still  maintaining  low  sidelobe  control  over  the  operating 
frequency  band  of  6-12  MHz. 

A  random  array  of  elements  spaced  over  a  large  aperture  is 
one  way  to  approach  this  problem.  However,  it  was  decided  not  to 
use  this  method  for  several  reasons.  First,  increased  signal 


6 


processing  is  necessary  to  reduce  the  sidelobes  to  levels 
comparable  to  a  periodic  array.  Secondly,  the  mainbeam  pattern 
and  sidelobe  levels  are  difficult  to  predict  during  all  operating 
situations  and  finally  the  array  was  to  retain  as  much  commonality 
as  possible  with  other  existing  HF  antenna  systems. 

A  linear  array  of  widely  spaced  (greater  than  half  wavelength) 
directional  elements  was  chosen  as  a  practical  solution  to  the 
problem.  Such  an  array  provides  a  wide  aperture  that  results  in  a 
narrow  beamwidth  in  azimuth  and  also  provides  good  sidelobe 
control.  A  wider  aperture  is  possible  using  subarrays  of  monopole 
elements  that  control  grating  lobes  as  the  array  is  scanned. 
Backpoles  placed  a  quarter  wavelength  behind  each  active  element 
provide  increased  directivity  of  the  array  in  the  forward 
direction.  While  a  backscreen  is  more  effective,  environmental 
concerns  and  colocation  with  other  experiments  on  the  selected 
site  u^rently  prohibit  use  of  a  large  backscreen.  The  geographic 
coverage  area  of  the  array  is  shown  in  Figure  1.  Note  that  the 
coverage  area  includes  a  large  overlap  of  the  area  associated  with 
the  nighttime  auroral  oval. 


2.2  Verona  Array  Description 

The  linear  HF  antenna  array  at  Verona  NY  is  composed  of  36 
subarrays  with  a  spacing  of  20  meters  between  centers  of  each 
subarray.  This  elemental  spacing  corresponds  to  0.8  wavelengths 
at  the  upper  design  frequency  of  12  MHz.  The  dimensions  and 
components  of  the  subarrays  are  shown  in  Figure  2.  The  main  idea 
in  the  subarray  design  process  was  to  use  a  given  number  of 
monopoles  as  the  basic  components  to  create  a  directional  element. 
Two  monopole  elements,  placed  a  half  wavelength  apart  (D2=12.5  m) , 
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in  a  broadside  configuration  provided  forward  directivity  and 
reduced  the  energy  received  from  signals  arriving  from  the 
direction  parallel  to  the  array.  The  6  meter  height  of  these  two 
elements  was  chosen  to  keep  the  element  length  near  a  quarter 
wavelength  of  the  operating  frequencies.  A  12  meter  backpole  was 
placed  a  quarter  wavelength  (D1  =  6.25  m)  behind  each  of  the  6 
meter  elements.  A  height  of  12  meters  was  chosen  to  keep  the 
backpoies  greater  than  a  quarter  wavelength  even  at  the  low  end  of 
the  operating  band.  The  four  monopoles  acting  together  as  a 
subarray  provide  forward  directivity  and  wide  nulls  at  angles  of 
t90’  and  130°.  The  spacing  between  the  subarrays  (D3)  was 
determined  from  the  constraint  of  reducing  the  grating  lobe  by  10 
dB  when  the  array  is  scanned  30’  at  a  frequency  of  12  MHz. 

A  schematic  of  the  linear  array  setup  is  shown  in  Figure  3. 
Several  parameters  are  included  in  this  figure  to  point  out  the 
major  features  of  the  array  including  its  length,  spatial 
orientation,  operating  frequency  band,  scan  limits,  and  minimum 

beamwidth . 
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3.  THEORETICAL  METHODS  AND  CONSIDERATIONS 

3.1  Calculation  of  Theoretical  Radiation  Patterns 

Far-field  azimuthal  radiation  patterns  were  calculated  at 
several  stages  of  the  antenna  array  development.  These  stages 
included  a  simple  monopole  element,  two  elements  phased 
broadside,  two  elements  phased  endfire,  a  four  element  subarray, 
and  the  36  element  linear  array. 

A  great  deal  of  consideration  was  given  to  the  effect  of  the 
ground  upon  the  vertical  field  patterns.  While  a  finite  ground 
might  have  little  effect  on  the  input  impedance  of  a  vertical 
monopole,  the  field  pattern  can  be  greatly  influenced  by  the 
ground  conductivity.  Fortunately  for  the  antenna,  (and 
unfortunately  for  the  experimenters) ,  the  land  at  Verona  is  for 
all  intents  and  purposes  a  swamp.  Ground  conductivity  for  this 
area  is  very  high.  A  conductivity  value  of  a  =  0.1  mhos/meter  and 
a  relative  dielectric  constant  of  e  =  30  were  used  for  all 
calculations.  These  values  were  chosen  from  literature  as  being 
consistent  for  a  marshy  ground1.  While  this  is  by  no  means  a 
perfectly  reflective  ground  it  improves  the  antenna  performance  at 
receive  angles  close  to  the  horizon.  In  addition,  a  ground  screen 
was  installed  in  front  of  the  array  to  ensure  a  stable  impedance 
for  the  antenna.  The  ground  screen  consisted  of  two  22  gauge 
aluminum  clad  wires  extending  75  meters  in  front  of  each  element. 
Field  calculations  were  made  with  the  assumption  of  an  infinite 
and  perfect  ground  plane. 

3.2  Simple  vertical  moncpcle  element. 

The  first  antenna  structure  investigated  was  a  vertical 
monopole  over  a  perfect  ground  plane.  This  antenna  has  a  radiation 
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pattern  similar  to  a  vertical  dipole  in  free  space  and  is 
therefore  omnidirectional  in  the  azimuthal  plane.  The  elevation 
field  pattern  for  a  vertical  monopole  has  been  extensively 
calculated  in  the  past  and  verified  for  nearly  all  situations2. 
Parameters  that  influence  the  pattern  are  ground  conductivity, 
element  length,  element  radius,  and  operating  frequency.  Elevation 
patterns  are  shown  in  Figure  4.  for  a  thin,  6  meter  monopole  over 
a  perfect  ground  plane  and  also  with  a  =  0.1  and  e  =  30.  Note 
that  the  loss  in  gain  for  the  antenna  over  a  finite  conducting 
ground  becomes  critical  only  at  angles  close  to  the  horizon, 
(greater  than  85°).  This  strongly  suggests  good  low  angle 
performance  for  the  Verona  linear  antenna,  which  allows  for  long 
range  OTH  radar  operation. 

3.3  Two  Element  Subarrays 

3.3.1  Two  Active  Elements  Phased  Broadside 

The  next  step  in  the  pattern  shaping  process  was  to  add  a 
second  monopole  and  phase  the  resultant  pair  to  increase  the  gain 
broadside  (0°)  while  placing  a  null  in  the  endfire  (90°, 270°) 
direction.  The  azimuth  pattern  for  a  two  element  broadside  array 
of  point  sources  with  D2  =  0.5  \  is  described  by  Eq.  1. 


E(<*>)  =  2  •  Sin  (k-D.*- Cos  (<*>)) 

EQ:  1 

“  2 

k  =  27T 

A 

The  pattern  resulting  from  Eq.  1  is  shown  in  Figure  5.  Note  the 


nulls  in  the  ±90°  directions  and  maximum  directivity  to  the  front 
and  rear  of  the  array. 


3.3.2  Active  Element  and  Backpole  Phased  Endfire 


10 


The  next  step  in  the  pattern  shaping  process  was  the  addition 
of  parasitic  backpoles  placed  a  quarter  wavelength  (D1  =  6.25  m) 
behind  each  active  element.  One  active  element  and  one  parasitic 


backpole  were  considered  as  an  endfire  array.  A  backscreen 
greater  than  a  quarter  wavelength  in  height  (at  6  MHz)  and  placed 
one  quarter  wavelength  behind  the  active  radiators  acts  as  a 
reflector.  This  reflector  increases  gain  in  the  forward  direction 
by  several  decibels  and  provides  nulling  of  any  signals  arriving 
from  behind  the  array.  Such  a  backscreen  would  be  over  12  meters 
in  height  and  700  meters  long.  Environmental  constraints  and 
other  on-going  experiments  conducted  at  the  Verona  test  site 
prevented  building  such  a  structure.  However,  the  antenna  gained 
significant  benefit  by  placing  a  passive  or  parasitic  element  one 
quarter  wavelength  behind  each  active  element.  These  "backpole" 
elements  are  connected  only  to  the  ground  screen  and  to  a  copper 
grounding  stake  driven  1.2  meters  deep  into  the  earth  beside  each 
backpole.  The  increase  in  electric  field  intensity  as  a  function 
of  azimuth  angle  <*>,  of  one  active  monopole  with  a  parasitic 


Rll  =  self-resistance  of  a  vertical  antenna  element 

R1L  =  effective  loss  resistance  of  a  single  element 

r12  +  lxi2  =  z12  =  Mutual  impedance,  elements  1  and  2. 

r22  +  3X22  =  z22  =  Backpoles  self-impedance. 

rm  =  arctan  X12  /  T2  =  arctan  x22 
R12  R22 
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!  dr  -  spacing  between  elements 

Note  in  Eq.  2,  that  as  the  impedance  of  the  back  pole  element 
( Z22 )  is  increased  the  effect  of  the  backpole  decreases  and  the 
gain  approaches  that  of  a  single  active  element  (G(0)=1).  The 
effectiveness  of  the  backpole  therefore  depends  to  a  large  extent 
upon  its  length.  Also,  it  can  be  shown  from  Eq.  2  that  when  the 
backpole  is  made  longer  than  X/4  it  acts  as  a  reflector,  and  when 
it  is  shorter  than  ,\/4  it  becomes  capacitive  and  behaves  as  a 
director.  Changing  frequency  will  alter  the  performance  of  the 
array  since  this  corresponds  to  changing  the  dimensions  of  the 
array  in  terms  of  wavelength.  Complete  cancellation  of  the  signal 
arriving  from  180 1  is  possible  only  if  dimensions  and  physical 
orientation  of  the  array  are  perfect.  Signals  arriving  from 
angles  not  exactly  180“  are  never  completely  cancelled.  This 
arrangement  does 

nowever,  predict  a  significant  increase  in  the  forward  (0°) 
directivity  and  a  decrease  in  directivity  to  the  rear  (180°).  The 

field  pattern  for  one  active  element  and  one  backpole  in  a  two 

element  endfire  configuration  array  of  point  sources  with  (D1  - 

,4)  is  described  by  Eq.  3  and  illustrated  in  Figure  6.  The 

pattern  shows  the  largest  gain  in  the  forward  direction  and 

indicates  a  null  in  the  ±180°  direction. 


Ex(0)  =  2 • Eo • Cos ( k • • Sin ( 0 ) ) 


EQ:  3 


0 

k 


Distance  in  wavelength  between  elements. 
Azimuthal  scan  angle  from  broadside. 

2  IT 


A 
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3.4  Design  of  Subarray  Dimensions 

The  last  step  in  elemental  beam  shaping  is  to  combine  the 
broad-side  and  end-fire  configurations  to  create  the  final 
subarray.  Combination  of  the  broadside  and  end-fire  patterns 
through  pattern  multiplication  results  in  the  subarray  pattern 
formed  by  the  two  active  elements  and  the  two  passive  backpoles. 
The  pattern  for  this  setup  is  shown  in  Figure  7.  The  main 
element  gain  is  shown  to  be  in  the  forward  direction.  Nulls  are 
located  in  the  ±  90°  and  180°  directions.  This  pattern  is 
considered  to  be  the  theoretical  element  pattern. 


3.5  Calculation  of  the  Array  Factor 

The  far-field  array  factor  for  a  linear  array  of  36  isotropic 
point  sources  is  described  by  the  square  of  Eq.  4.  The  field 
pattern  for  the  array  factor  with  N  =  36,  D3  =  0.8  and  <pQ  =  0 
is  shown  by  the  plot  in  Figure  8.  This  pattern  is  for  a  uniform 
array  where  all  36  elements  are  weighted  equally.  The  peak  of  the 
first  side-lobe  is  13.2  dB  below  the  main  beam.  The  3  dB 
beamwidth  is  approximately  2.5’. 
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3.6  Pattern  Multiplication  of  Array  Factor  and  Subarray 
Pattern 

A  broadside  (O’  scan)  far-field  array  pattern  for  the  thinned 
linear  antenna  array  is  also  shown  in  Figure  8.  It  was  calculated 
by  multiplying  the  subarray  (element)  pattern  (also  shown  in 
Figure  8.)  by  the  array  factor.  Uniform  amplitude  weighting  was 
used  in  calculating  this  pattern.  Characteristic  13  dB  sidelobes 
’^e  seen  on  either  side  of  the  main  beam.  However,  now  a  broad 
subarray  pattern  has  been  placed  over  the  array  factor.  This 
results  in  higher  directivity  in  the  forward  direction  and  nulls 

at  angles  of  ±90°,and  180°. 

3.7  Grating  Lobe  Reductions 

A  plot  of  the  array  factor  when  the  array  is  scanned  off 
broadside  by  20°  is  shown  in  Figure  9.  Note  the  grating  lobes  that 
have  moved  into  the  pattern.  They  are  located  at  -120°,  -80°  and 
160 ’  and  are  equal  in  amplitude  to  the  main  beam  of  the  pattern 
now  at  20’.  Another  plot  in  Figure  9  shows  the  total  array 
pattern.  Notice  that  the  amplitude  of  the  grating  lobes  relative 
to  the  main  beam  has  been  lowered  significantly  by  the  element 
pattern  with  the  nearest  one  at  -80°  now  lower  by  16  dB.  This 
reduction  in  grating  lobes  of  the  array  by  using  directive 
subarrays  has  permitted  the  aperture  of  the  array  to  be  increased 
from  437.5  meters  for  the  usual  //2  spacing  between  elements  to 
700  meters  or  0.8  wavelengths  between  elements.  This  corresponds 
to  a  decrease  in  beamwidth  from  4°  to  2.5’  at  the  12  MHz  design 
f  requency . 
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4.  Calculations  Using  Numerical  Electromagnetics  Code  (NEC) 

4.1  Introduction  to  NEC 

The  Numerical  Electromagnetics  Code  (NEC) 4  is  a  computer 
program  that  uses  the  Method  of  Moments  technique  to  predict  the 
electromagnetic  response  of  antennas  and  other  metal  structures. 
The  physical  structure  to  be  analyzed  is  modeled  using  wire 
segments  to  approximate  the  correct  shape  and  electrical 
characteristics.  The  required  integral  equations  are  then  solved 
to  determine  the  currents  on  an  electromagnetically  excited  wire 
or  structure.  Similar  to  the  theoretical  calculations  a  step  by 
step  process  using  NEC  was  employed  to  model  various  monopole 
combinations.  These  structures  were  then  combined  to  obtain  the 
subarray  and  finally  the  complete  array.  All  computations  were 
carried  out  on  a  Digital  Equipment  Corp.  Micro-Vax  II  computer. 


4.2  Far-field  Patterns  Calculated  with  NEC 
4.2.1  Vertical  Monopole  Over  perfect  Ground 

A  simple  vertical  monopole  antenna  was  modeled  first  using 
NEC.  The  height  of  the  antenna  is  6  meters  and  the  element  is 
driven  by  a  1  volt  source  close  to  the  ground.  A  calculation 
assuming  a  perfect  ground  plane  was  made  to  create  a  baseline  and 
verify  that  the  model  is  accurate.  Next  an  element  over  a  finite 
ground  plane  was  modeled  using  NEC.  An  option  available  with  the 
NEC  program  to  use  Sommer f eld  Integrals  was  employed.  This  method 
allows  simulation  of  a  structure  very  close  to  the  ground. 

Various  ground  plane  configurations  and  different  frequencies  were 
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examined.  Resultant  patterns  are  shown  in  Figure  10.  for  both  the 
perfectly  conducting  ground  and  a  ground  with  finite  conductivity 
representative  of  the  earth  at  the  Verona  test  site  (a=0.1 
mhos/meter,  e=30) . 

4.2.2  Two  Active  Elements 

A  simulation  of  two  active  elements  phased  for  a  broad-side 
configuration  was  computed  and  the  field  patterns  calculated. 

This  arrangement  generates  the  common  " figure-eight"  pattern  with 
nulls  at  ±90°  in  relation  to  broadside.  The  resultant  azimuth 
pattern  is  shown  in  Figure  11. 

4.2.3  Monopole  with  Backpole 

Similar  to  the  process  using  theoretical  methods,  a  backpole 
was  paired  with  an  active  element  to  reduce  the  gain  in  the  back 
half  plane.  The  backpole  in  the  NEC  model  is  again  located  a 
quarter  wavelength  behind  the  active  element.  Ground  wires  are 
connected  from  the  front  element  to  the  backpole.  The  backpole  is 
not  driven  with  a  voltage  but  purely  reacts  in  the  passive  manner 
as  described  earlier  in  Section  3.3.2.  This  setup  as  in  the 
earlier  theoretical  case  results  in  an  endfire  configuration.  The 
resultant  pattern  is  shown  in  Figure  12.  Note  that  nulls 
approximately  6  dB  below  the  forward  direction  are  formed  at  180°. 
In  this  and  the  following  simulations  the  a  perfectly  conducting 
ground  plane  model  was  used  since  only  the  azimuthal  variation  is 
under  consideration. 

4.2.4  Full  Subarray 

The  full  subarray  of  two  active  elements  and  two  backpoles  was 
modeled  next.  Both  active  elements  are  driven  in  phase  by  1  volt 
sources.  The  resultant  pattern  for  azimuthal  scan  is  shown  in 
Figure  13.  This  is  the  element  pattern  for  the  array.  It  is  this 
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pattern  that  can  be  multiplied  by  the  array  factor  to  create  the 
full  array  pattern.  The  subarray  pattern  calculated  using  NEC 
shows  nearly  15  dB  of  increased  directivity  in  the  forward 
direction  (O’)  compared  to  behind  the  array  (130’). 

4.2.5  Model  of  the  Full  Array 

The  final  modeling  effort  is  naturally  the  full  array  of  36 
subarrays  and  calculation  of  the  composite  element-array  far-field 
patterns.  Unfortunately,  our  version  of  the  NEC  code  cannot 
accommodate  structures  as  large  as  the  full  array  of  36  elements. 
The  largest  array  that  could  be  handled  in  the  program  was  15 
subarrays.  The  limitation  of  15  subarrays  should  only  decrease  the 
gain,  beamwidth,  and  number  of  lobes  in  the  patterns.  The  overall 
response  of  the  15  element  linear  array  model  should  still  be  a 
good  representation  of  the  36  element  array. 

All  active  elements  were  fed  with  a  one  volt  source  with  no 
phase  difference  between  elements  causing  the  main  beam  of  the 
pattern  to  occur  in  the  broadside  direction.  The  result  is 
identical  to  a  plane  wave  arriving  that  is  oriented  parallel  to 
the  array.  The  resultant  azimuthal  radiation  pattern  is  shown  in 
Figure  14.  In  order  to  determine  the  expected  side  lobe  behavior 
of  the  array,  an  amplitude  weighting  taper  was  used  during 
calculation  of  the  field  pattern.  This  is  done  by  factoring  each 
of  the  voltage  sources  with  a  32  dB  Cosine  Squared  amplitude 
weight.  The  cosine  squared  function  is  only  one  of  many  possible 
weightings  and  was  chosen  because  it  does  not  broaden  the  main 
beam  as  much  as  other  comparable  weighting  functions.  The  pattern 
for  this  simulation  is  shown  in  Figure  15.  Notice  that  the 
sidelobe  levels  have  been  reduced  from  the  uniform  level  of  13.2 
dB  shown  in  Figure  14  to  32  dB  below  the  main  beam.  The  main  beam 
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has  broadened  by  1  degree  and  about  3  dB  of  main  beam  amplitude 
loss  has  occurred.  This  loss  and  beam  spread  will  cause  a 
reduction  in  array  resolution  but  it  does  not  prevent  adequate 
system  performance. 

5.  Field  Measurements 

5.1  Introduction  and  Techniques 

The  HF  test  facility  using  the  linear  antenna  array  at  Verona 
NY  is  currently  in  operation.  Calibration  and  testing  are 
complete  and  some  preliminary  clutter  data  have  been  collected. 
The  pertinent  pattern  measurements  made  were  the  subarray  pattern 
and  the  full  36  subarray,  linear  antenna  pattern.  The  essential 
antenna  information  gained  from  these  measurements  includes  the 
gain,  average  sidelobe  level,  beamwidth,  grating  lobe  response, 
front  to  back  ratio  and  scan  direction  accuracy. 


5.2  Measurement  of  the  Subarray  Pattern 

The  subarray  pattern  was  determined  by  moving  a  continuous 
wave  probe  across  the  subarray  aperture  and  plotting  the  received 

signal  strength  as  a  function  of  azimuth.  Measurements  were  also 
made  behind  the  array  to  determine  the  element  directivity  to  the 
rear  (180°).  Figure  16  shows  the  azimuthal  element  pattern  for  all 
the  elements  at  a  frequency  of  10.2  05  MHz.  Although  not  shown  on 
this  plot  the  front  to  back  ratio  was  also  measured  and  is 
approximately  12  dB.  Note  that  the  subarray  or  element  gain  falls 
off  rapidly  for  signals  arriving  at  angles  greater  than  25°  from 
broadside.  This  is  the  mechanism  that  permits  the  array  to  scan 
±30°  with  the  grating  lobes  reduced  to  at  least  10  dB  below  the 
main  beam. 

5.3  Measurement  of  the  Array  Pattern 
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The  field  pattern  for  the  full  array  was  measured  by  scanning 
the  array  main  beam  past  a  CW  signal  transmitted  from  Ava  NY.  Ava 
is  located  21  miles  away  at  a  bearing  from  Verona  of  31°  east  of 
true  north  which  corresponds  to  21°  east  of  the  array  boresight. 
The  power  of  the  transmitted  signal  was  ten  kilowatts  and  the 
antenna  used  was  a  horizontal  log  periodic.  The  resultant  full 
array  pattern  is  shown  in  Figure  17.  This  pattern  was  obtained 
using  measured  data  and  applying  a  32  dB  cosine  squared  amplitude 
weighting  function.  Sidelobe  response  to  the  test  signal  is  less 
than  32  dB  RMS.  Note  that  the  beamwidth  is  greater  than  the  2.5 
degrees  available  without  the  weighting  taper  but  is  still  only 
approximately  4  degrees.  This  method  of  course  does  not  account 
for  the  response  of  the  subarray  pattern  but  only  the  subarray 
pattern  at  a  single  bearing.  Therefore,  the  actual  gain  will  fall 
off  as  the  scan  angles  become  greater  than  20°  as  shown  in  Figure 

16  by  the  subarray  pattern. 

6.  Array  Description 

The  elements  and  pre-amplifiers  for  the  linear  antenna  array 
were  designed  and  built  at  RADC  "in-house".  The  elements  are 
constructed  from  electrical  supply  house  materials  including,  1 
inch  diameter  rigid  aluminum  conduit,  pvc  pipe,  pvc  flanges,  and 
pvc  junction  boxes.  A  diagram  of  the  element  is  shown  in  Figure 
18  . 

The  signal  output  of  each  active  monopole  element  is  fed  to  a 
power  combiner  located  halfway  between  the  elements.  The  combined 
output  is  then  band-pass  filtered  and  sent  to  the  pre-amplifier. 
Pre-amplifiers  were  deemed  essential  to  overcome  the  sum  of  losses 
from  cables,  filters,  and  power  combiners,  which  was  calculated  to 
be  15-20  dB.  The  gain  of  the  pre-amps  was  chosen  to  be  20  dB.  The 
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output  of  the  pre-amplifier  is  carried  through  RG/215  coaxial 
cable  to  the  receivers  located  in  a  central  receive  building.  The 
building  is  located  slightly  off  the  middle  of  the  array  and  30 
meters  behind  the  array.  The  RG/215-U  used  for  the  co-axial  feed 
cable  is  low  loss  double  shielded  cable  with  an  outside  armor 
sheath.  A  unique  aspect  of  this  array  is  that  eacn  one  of  the 
cables,  connecting  a  subarray  back  to  one  of  the  receivers,  are  of 
different  length.  The  longest  cable  is  approximately  1600  feet, 
while  the  shortest  is  150  feet.  Phase  and  amplitude  differences 
for  all  channels  are  accounted  for  and  adjusted  digitally  during 
the  beamforming  process.  Figure  19  details  the  subarray 
configuration  including  the  elements,  power  combiners,  filters, 
and  pre-amplifiers. 

7.  Discussion 

Once  the  required  operating  properties  of  the  antenna  ie. 
gain,  beamwidth,  frequency  band,  etc.  were  specified,  the  antenna 
design  effort  and  testing  for  this  system  consisted  of  a  three 
stage  process.  First,  antenna  theory  was  used  to  design  an 
antenna  that  would  meet  these  specifications.  Second,  numerical 
methods  were  used  to  test  the  theoretical  designs  and  ensure  that 
the  designs  met  the  performance  specifications.  Finally,  pattern 
measurements  of  the  fielded  antenna  were  made  to  verify  that  the 
array  performed  as  designed.  Table  1.  contains  some  of  the  major 
required  specifications  and  the  results  obtained  from  each  of  the 
three  techniques  used.  Although  there  are  some  minor  differences, 
in  general,  the  agreement  between  the  numbers  is  quite  good  and 
shows  that  the  use  of  theory  and  modeling  techniques  leads  to  the 
fabrication  of  an  antenna  with  the  properties  needed  for 
conducting  meaningful  experiments. 
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Table  1.  Parameter  Values  from  Each  of  Three  Techniques 


Parameter 

Theoretical 

NEC 

FIELD 

Calculations 

Calculations 

Measurements 

26 

24 

23 

BW  (dB) 

2 . 5 

2 

«  2.5 

AVE.  SLL  (dB) 

32 

s  30 

32 

F/B  ratio  (dB) 

>  20 

15 

12-14 
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Theoretical  pattern  for  one  active  element  and  one 
backpole  phased  end-fire. 
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Figure  1C>.  Measured  azimuthal  radiation  pattern  for  the  fielded  four 
element  subarray. 
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Abstract 


HF  skywave  signals  are  propagated  beyond  the  Earth's 
curvature  by  refraction  within  the  ionosphere.  As  a 
result,  over-the-horizon  radars  have  an  advantage  in 
providing  advanced  notice  of  targets  of  interest  as  well 
as  a  much  enhanced  surveillance  area  when  compared  to 
conventional  line  of  sight  radars. 

After  passing  through  the  ionosphere  signal  wavefronts 
may  experience  perturbations  in  both  amplitude  and  phase. 
The  corrupted  wavefront  propagates  and  impinges  on  an 
antenna  aperture.  If  we  consider  the  ionosphere  and 
antenna  as  a  system,  these  ionospherically  induced  errors 
will  negatively  influence  the  antenna's  response  and 
consequently  affect  radar  performance. 

We  have  begun  a  measurement  program  to  quantify 
statistically  the  effects  of  ionospherically  induced 
errors  on  apertures  up  to  10  km  in  length.  A  complex 
correlation  of  the  signals  received  along  spaced  coherent 
receivers  is  used  to  form  a  mutual  coherence  function 
which  is  useful  in  evaluating  antenna  performance.  We 
examine  limitations  on  the  mutual  coherence  function  in 
two  dimensions  (laterally  in  space,  and  time)  for  four 
different  propagation  paths.  The  various  propagation 
modes  are  easily  separated  using  a  12.5  KHz  bandwidth 
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swept  frequency  continuous  waveform  (SFCW) ;  emphasis  is 
placed  on  F2  low  ray  propagation. 


1.  Introduction 

Temporal  coherence  affects  the  spectral  width  of  a 
signal  source  and  spatial  coherence  affects  the  apparent 
size  of  a  signal  source.  In  our  experiment,  the  source, 
which  simulates  a  target  return,  is  coherent  in  both 
these  domains  during  a  measurement  period.  The 
intervening  medium,  namely  the  ionosphere,  is  responsible 
for  degrading  the  image  of  the  source  as  seen  by  a  six 
element  10-km  receive  aperture.  Ideally,  the  ionosphere 
would  reflect  an  exact  image  of  this  perfect  source. 
However,  variations  in  the  spatial  distribution  of 
electron  density  cause  the  image  to  appear  to  have  come 
from  and  extended  source.  Also,  the  time  variability  of 
this  medium  produces  effects  similar  to  those  produced  by 
a  source  of  finite  spectral  width. 

If  the  errors  induced  by  the  ionosphere  are  considered 
as  antenna  errors,  the  aperture  becomes  only  partially 
coherent.  Techniques  developed  earlier111  may  be  used  to 
analyze  the  performance  of  antennas  with  partially 
coherent  illumination.  The  corrupted  wavefront  is  then 
described  in  terms  of  a  mutual  coherence  function. 
Several  authors,  for  example,  Altshuler  et.  al.t2]  and 
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Mitchell131,  have  shown  the  degradation  in  the  performance 
of  antennas  with  such  illumination. 

In  figure  1,  a  computer  simulation  shows  the  ideal  and 
corrupted  antenna  patterns  of  a  100  element  10-km  linear 
array  of  isotropic  radiators.  Interelement  spacing  is  2 
wavelenghts  at  6  MHz,  and  a  Hamming  amplitude  taper  is 
applied.  The  ideal  pattern  shows  the  expected  0.3-degree 
beamwidth  and  -43  dB  side  lobe  level.  In  the  case  of  the 
corrupted  pattern,  a  0.7-radian  rms  phase  deviation  and 
a  1  km  phase  error  correlation  length  were  used.  These 
errors  are  assumed  to  have  originated  in  the  ionosphere, 
however  we  consider  the  effective  antenna  performance  as 
though  they  were  truely  antenna  errors.  The  averaged 
power  pattern  of  this  antenna/ ionosphere  system 
experiences  gain  loss,  beam  broadening,  and  an  increase 
in  side  lobe  level. 

Data  have  been  recorded  to  quantify  the  spatial  and 
temporal  variability  of  the  ionosphere  with  diurnal, 
seasonal,  and  latitude  changes.  This  ongoing 
experimental  program  is  known  as  the  Extended  Coherent 
Aperture  Project  (ECAP) .  In  this  paper,  we  present  four 
measured  normalized  mutual  coherence  functions  for  winter 
nighttime  conditions  using  F2  low  ray  propagation.  These 
data  are  considered  typical  for  the  conditions  stated, 
although  they  are  snapshots  or  samples  of  this  function. 


2 .  Mutual  Coherence  Function 

The  mutual  coherence  function  is  a  complex  correlation 
of  the  signals  across  an  aperture,  which  can  be  written 
after  Born  and  Wolf[4]  as: 

Y{x1.x2,x)  =r12  (t)  =<V^  ( t+x )  y2*  ( t)  >  (1) 

Vvhere  V  refers  to  a  scalar  field  component,  x,  and  x2 
refer  to  position,  t  refers  to  time,  and  t  represents 
a  time  delay.  The  sharp  brackets  indicate  an  infinite 
time  average  and  *  denotes  complex  conjugation.  A  time 
stationary  ergodic  process  is  assumed.  In  optics,  it  is 
convenient  to  define  a  time  average  since  practical 
devices  can  only  measure  such  an  average.  However,  at 
HF,  measurements  of  the  instantaneous  fields  are  easily 
made  so  that  an  ensemble  average  can  be  computed.  We 
choose  to  rewrite  (1)  as 

r(6,t)  =<V(x+b,  t+x)  V(x,  t)  >  (2) 

where  6  refers  to  a  spatial  delay  and  the  sharp 
brackets  now  indicate  an  ensemble  average.  That  is,  the 
average  is  performed  over  a  group  of  measured  wavefronts. 
Spatial  and  temporal  stationarity  are  implied  in  (2) 
since  T(6,t)  does  not  depend  on  space  or  time  origin. 
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Each  computed  value  of  r(6,x)  its  normalized  by  the 
average  value  of  r{0,0)  for  six  individual  receiving 
elements . 

Spatial  coherence  effects  are  described  by  r(6,0)  and 
may  be  measured  perpendicular  to  the  direction  of 
propagation  with  an  array  of  antennas.  F(6,0)  is  also 
called  the  lateral  coherence  function153  and  can  be 
related  to  the  angular  spectrum  of  incident  plane  waves. 
Temporal  coherence  effects  are  described  by  r(0,x)  and 
can  be  related  to  the  power  spectrum  of  the  fluctuations 
(Weiner-Khinchine  theorem) .  We  are  primarily  interested 
in  the  lateral  coherence  function  r(6,0)  because  it  is 
used  to  predict  the  average  pattern  of  an  arbitrary 
antenna113.  However,  the  two  dimensional  mutual  coherence 
function  r(6,x)  is  useful  in  revealing  limitations  on 
coherent  integration  time  for  operational  radars. 

Figure  2  shows  how  the  various  quantities  of  interest 
are  interrelated.  The  heavy  two-way  arrow  represents  a 
Fourier  transform  pair  and  the  one-way  arrow  represents 
a  nonreversible  path.  Nonreversible  implies,  for 
example,  that  the  aperture  field  A(x)  cannot  be  uniquely 
determined  from  knowledge  of  the  autocorrelation  alone. 

With  no  wavefront  perturbation,  T(6,0)  is  unity  and 
the  effective  antenna  transfer  function  x  (x)  is  the 
unaltered  autocorrelation  of  the  field  distribution.  The 
average  power  pattern  P(Q)  (  0  being  the  azimuthal 
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angle)  is  then  simply  the  Fourier  transform  of  the 
coherent  transfer  function.  Alternately,  the  Fourier 
transform  of  the  lateral  coherence  function  T(v,0)  (  v 

is  an  angular  coordinate  and  the  A  denotes  a  Fourier 
transform)  is  a  dirac  delta  function  which  must  be 
convolved  with  |l?(0)  |2  to  produce  the  average  power 
pattern.  In  this  case,  the  average  power  pattern  is  the 
expected  coherent  power  pattern. 

With  wavefront  perturbation,  P(0)  can  also  be 
considered  in  two  ways.  First,  the  average  power  pattern 
is  the  Fourier  transform  of  the  product  of  T(6,0)  and 
<A(x) A* (x+6) >  .  Or  equivalently,  the  average  power 
pattern  is  the  convolution  of  the  ideal  power  pattern 
with  f(v,0)  .  Knowledge  of  r(6,0)  and  its  statistical 
variation  is  necessary  to  predict  average  antenna 
performance . 

To  determine  the  mutual  coherence  function  for  an 
ionospheric  propagation  path,  one  can  transmit  a  known 
coherent  (space,  time)  signal  and  directly  measure  the 
complex  correlation  of  the  reflected  skyw’ave  signal  along 
a  series  of  spaced  aerials.  We  then  assign  these 
measured  errors  to  the  aperture,  and  treat  the  aperture 
as  being  partially  coherent.  One  could,  in  principle, 
measure  the  angular  power  spectrum  directly  and  Fourier 
transform  this  spectrum  to  arrive  at  the  lateral 
coherence  function.  For  an  aperture  approaching  10  km  in 
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size,  the  former  represents  a  more  tractable  approach  mj 
is  the  one  used  here. 

3 .  Experiment 

As  shown  in  Figure  3,  and  east-looking  six  element 
receive  aperture  was  formed  by  the  addition  of  two 
portable  elements  to  four  fixed  elements  of  the  existing 
Wide  Aperture  Research  Facility  (WARF)  located  in  Los 
Banos,  CA.  The  portable  receive  elements  are 
relocateable  and  convey  demodulated  SFCW  signals  at 
baseband  to  the  WARF  control  center  by  FM  modulation  of 
a  UHF  carrier.  The  fixed  elements  are  wired  directly  to 
the  WARF.  This  array  is  easily  reconfigured  for  north- 
look  measurements  using  four  additional  fixed  elements 
and  by  moving  the  two  portable  receive  elements,  as  shown 
in  Figure  3.  Each  monopole  element  is  a  5.5-m  whip  fed 
against  a  set  of  eight  15.2-m  radial  wires.  The  location 
of  each  whip  was  surveyed  to  a  relative  accuracy  of  ±5 
cm. 

Assuming  that  the  spatial  statistics  depend  only  on  the 
relative  positions  of  the  elements  and  not  on  their 
absolute  position  (spatial  stat ionarity ) ,  we  generate 
correlations  for  15  separations  using  six  coherent 
receivers.  Note,  in  figure  3,  that  the  array  was 
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desiqned  to  have  some  redundant  spacings  in  order  to 
verify  the  above  assumption. 

The  measurement  process  includes  data  collection  for 
three  two-week  periods  throughout  the  year. 
Specifically,  data  were  collected  for  an  autumnal  equinox 
(September  1989) ,  a  winter  solstice  (January  1990) ,  and 
a  summer  solstice  (June  1990) .  Only  data  recorded  for 
the  first  two  periods  have  been  processed  at  the  time  of 
this  writing.  In  each  period,  the  first  week  is 
dedicated  to  east-look  measurements,  followed  by  a  week 
of  north-look  measurements.  A  complete  diurnal  cycle  is 
examined.  Prior  to  data  recording,  a  wide  sweep  (6-30 
MHz)  ionogram  is  recorded.  This  recording  is  a  useful 
presentation  of  range  or  time  delay  versus  frequency;  it 
reveals  the  various  ionospheric  layers  or  modes  that  are 
present . 

During  the  measurements,  a  ground  wave  calibrator  (GWC) 
is  situated  approximately  5  km  in  front  of  the  receive 
array.  It  transmits  an  SFCW  signal  and  is  received  in 
parallel  with  a  similar  over-the-horizon  (OTH)  signal. 
The  GWC  signal  is  common  to  all  array  elements  and  is 
used  as  a  coherent  reference  for  the  array. 

Figure  4  shows  the  geographical  locations  of  the  WARF 
receiver  and  four  remote  OTH  signal  sites.  For  the  east- 
look  measurements,  a  transmitter  in  Lost  Hill,  CA  excites 
either  a  transponder  in  White  Sands  Missile  Range,  NM,  or 
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a  transponder  in  Little  Rock  Air  Force  Base,  AR.  The 
ranges  to  the  WARF  are  1416  and  2585  km  respectively. 
For  the  north-look  measurements,  a  remote  transmitter  is 
situated  in  Langley,  BC  and  Yellowknife  NWT,  both  in 
Canada.  The  ranges  to  the  WARF  are  1336  and  2851  km 
respectively. 

A  12.5  KHz  bandwidth  SFCW  signal  was  used  at  each  of 
the  four  remote  sites  providing  a  range  resolution  of  12 
km.  This  resolution  is  quite  adequate  to  separate  the 
various  time  delayed  multipath  components  that  may  exist. 

Doppler  resolution  is  achieved  by  Fourier  transforming 
across  32  consecutive  sweeps.  A  waveform  repetition 
frequency  of  20  Hz  was  selected.  This  implies  1.6 
seconds  of  data  are  used  to  separate  the  GWC  and  OTH 
signals  in  a  range-Doppler  (RD)  map.  The  GWC  and  OTH 
signals  are  offset  in  both  time  and  frequency  to  allow 
easy  detection  in  the  RD  map.  Hanning  weights  are 
applied  in  both  range  and  Doppler  to  suppress  processing 
sidelobes.  The  amplitude  and  phase  of  the  OTH  and  GWC 
signals  are  interpolated  within  the  RD  map  for  a  better 
estimate  of  their  values. 

4.  Results 

The  results  presented  here  represent  four  snapshots  of 
the  two  dimensional  mutual  coherence  function  during  the 
winter  solstice  (January  1990) ,  one  snapshot 
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corresponding  to  each  of  the  four  remote  sites.  Each 
function  was  computed  during  nighttime  conditions  using 
an  F2  low  ray  propagation  path.  An  average  Kp  index  of 
2  was  reported  at  the  time,  which  is  indicative  of 
relatively  quiet  ionospheric  conditions. 

Figures  5,  6,  and  7  display  measured  mutual  coherence 
functions  for  mid-latitude  propagation  paths,  whereas 
figure  8  displays  this  function  for  an  auroral  path. 
Each  of  these  curves  are  computed  from  a  file  containing 
approximately  13  minutes  of  data  during  the  universal 
times  indicated  in  the  figures.  The  averaging  was 
performed  on  an  ensemble  of  wavefronts,  each  corrected 
for  spherical  curvature  and  slowly  moving  ionospheric 
tilts. 

For  each  sampled  wavefront  a  coherence  factor  was 
defined.  This  factor  is  the  ratio  of  actual  power  that 
the  six  element  array  would  receive  to  the  maximum  power 
if  the  errors  were  absent.  Only  the  wavefronts  that  had 
a  coherence  factor  greater  than  the  median  for  the  file 
were  used  to  produce  these  curves.  In  this  way,  the 
malignant  effects  of  Faraday  rotation  or  low  signal  to 
noise  ratio  were  removed  form  the  results. 

Most  evident  in  the  figures  is  that  coherence  degrades 
with  range  for  both  the  north  and  east  look  directions. 
The  coherence  is  reduced  significantly  when  looking 
towards  the  long  range  north-look  site,  possibly  due  to 
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interaction  with  the  auroral  trough.  There  were  several 
occasions  when  coherence  from  the  Yellowknife  site 
increased  noticeably  during  the  late  afternoon. 

The  mutual  coherence  functions  in  figures  5  through  8 
show  the  temporal  properties  of  ionospheric  errors  as 
well.  In  general,  temporal  coherence  appears  unrelated 
to  spatial  coherence.  Limitations  on  coherent 
integration  time  may  be  inferred.  A  comprehensive 
analysis  of  the  seasonal  and  diurnal  variation  of  the 
mutual  coherence  function  and  its  effects  on  radar 
antenna  performance  will  follow. 
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ANGLE  (DEGREES) 


Figure  1.  Comparison  between  coherent  and  partially 

coherent  antenna  patterns  for  a  100  element 
10  km  linear  array. 
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Figure  2.  Relationships  between  various  parameters 
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Figure  3.  Receive  element  and  ground  wave  calibrator 
locations . 
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re  4.  Geographical  locations  of  WARF  receiver  and 
remote  sites. 
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Figure  7.  Mutual  coherence  function  computed  from 
signal  transmitted  via  Langley,  BC. 


Yellowknife  (MKT)  Trans.,  7.93  KHz,  19  Jan  90,  06:50:32  UT 


Figure  8.  Mutual  coherence  function  computed  from 

signal  transmitted  via  Yellowknife,  NWT  site 


ADAPTIVE  NULLING  OF  TRANSIENT  ATMOSPHERIC  NOISE 
RECEIVED  BY  AN  HF  ANTENNA  ARRAY 
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The  MITRE  Corporation 
Bedford,  MA  01730 


Summary 

Impulsive  noise  transients  at  HF.  presumably  caused  by  atmospheric  lightning 
discharge,  can  be  a  severe  source  of  interference  that  limits  the  sensitivity  of  an  HF 
receiver  array  that  may  be  used  in  direction-finding  or  radar  applications.  To  the  extent 
that  such  energy  subtends  a  small  angle  at  the  receiving  array  aperture  and  is  received 
in  the  sidelobe  region,  it  would  be  useful  to  suppress  this  source  of  noise  by  adaptive 
spatial  nulling.  We  have  recently  examined  records  of  elemental  signals  recorded  from 
a  linear  HF  antenna  array  in  which  the  occurrence  of  transient  atmospheric  noise  is 
evident.  Closer  examination  of  some  of  the  stronger  transients  reveals  that  they  are 
highly  directional  and  provide  a  naturally  occurring  source  of  energy  that  is  essentially 
phase-coherent  across  the  array.  These  observations  are  based  on  the  use  of  a  very 
simple,  iterative,  retrodirective  beamforming  algorithm  that  calculates  an  estimate  of  the 
phase  steering  vector  corresponding  to  the  angular  direction  of  the  observed  transient. 
This  simple  technique  is  useful  not  only  for  establishing  the  array  weight  vector  to  be 
used  for  nulling  the  interference,  but  also  for  checking  the  phase  coherence  and 
possibly  compensating  the  array  for  calibration  or  propagation  path  phase 
perturbations.  Results  of  application  of  this  technique  to  the  recorded  field  test  data  will 
be  presented  and  the  performance  of  the  retrodirective  nulling  technique  will  be 
compared  with  an  optimal  baseline  established  by  the  computationally  more  costly 
method  of  adaptive  sidelobe  cancellation  using  a  constrained  least-squares  algorithm 
employing  reduced-rank  principal  components  analysis.  Results  will  be  compared  in 
terms  of  the  amount  of  suppression  of  the  transient  noise  and  by  the  impact  on  the 
adapted  array  pattern. 
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ADAPTIVE  NULLING  OF  TRANSIENT  ATMOSPHERIC  NOISE 
RECEIVED  BY  AN  HF  ANTENNA  ARRAY 


Dean  O.  Carhoun 
The  MITRE  Corporation 
Bedford,  MA  01730 

Abstract 

Impulsive  noise  transients  at  HF,  presumably  caused  by  atmospheric 
lightning  discharge,  can  be  a  severe  source  of  interference  that  limits  the  sensitivity 
of  an  HF  receiver  array  that  may  be  used  in  direction-finding  or  radar  applications. 
To  the  extent  that  such  energy  subtends  a  small  angle  at  the  receiving  array  aperture 
and  is  received  in  the  sidelobe  region,  it  would  be  useful  to  suppress  this  source  of 
noise  by  adaptive  spatial  nulling.  We  have  recently  examined  records  of  elemental 
signals  recorded  from  a  linear  HF  antenna  array  in  which  the  occurrence  of  transient 
atmospheric  noise  is  evident.  Closer  examination  of  some  of  the  stronger  transients 
reveals  that  they  are  highly  directional  and  provide  a  naturally  occurring  source  of 
energy  that  is  essentially  phase-coherent  across  the  array.  These  observations  are 
based  on  the  use  of  a  very  simple,  iterative,  retrodirective  beamforming  algorithm 
that  calculates  an  estimate  of  the  phase  steering  vector  corresponding  to  the  angular 
direction  of  the  observed  transient.  This  simple  technique  is  useful  not  only  for 
establishing  the  array  weight  vector  to  be  used  for  nulling  the  interference,  but  also 
for  checking  the  phase  coherence  and  possibly  compensating  the  array  for 
calibration  or  propagation  path  phase  perturbations.  Results  of  application  of  this 
technique  to  the  recorded  field  test  data  will  be  presented  in  terms  of  the  amount  of 
suppression  of  the  transient  noise  and  the  performance  of  the  retrodirective  nulling 
technique  will  be  compared  with  an  optimal  baseline  established  by  the 
computationally  more  costly  method  of  adaptive  sidelobe  cancellation  using  a 
constrained  least-squares  algorithm  employing  reduced-rank  principal  components 
analysis. 
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1.0  Adaptive  Nulling  Fundamentals 


Assume  a  linear  array  of  N  equally  spaced  omnidirectional  antenna  elements, 
with  one-half  wavelength  element  spacing,  that  is  used  to  form  a  narrow  beam  in  a 
designated  direction.  Assume  that  T  complex  digital  samples,  representing  an 
analytic  signal,  are  recorded  from  each  array  element  and  let  X  denote  a  T  x  N 
matrix  of  elemental  array  data  in  which  the  j^1  column,  xj,  consists  of  T  samples 
recorded  from  the  j^1  element  of  the  array.  Denote  by  wq  the  N  x  1  quiescent 
weight  vector  used  to  form  and  point  a  beam  in  the  desired  direction.  Assume  that 
M  <  N  array  elements  are  used  as  the  auxiliary  elements  of  an  adaptive  sidelobe 
canceller  used  for  adaptive  nulling.  Denote  the  auxiliary  array  data  matrix  by  A  and 
its  adaptive  weight  vector,  to  be  determined,  as  w.  Denoting  the  mainbeam  signal 
as  y  =  Xwq,  the  sidelobe  canceller  output,  z,  is  expressed  as  z  =  y-Aw  where  w  is 
obtained  by  minimizing  z^z  (H  denotes  conjugate  transpose).  The  adaptive  weight 
vector  is  determined  as  (see  [1]) 


w 


=  (AHAr!AHy 


(1) 


from  which  it  follows  that  the  adapted  output  is 


z  =  y-A(A^A)  lA^y  =  PA-*-y- 


(2) 


Pa^-  =  I-A(A^A)-^  A^  is  a  projection  operator  that  projects  the  mainbeam  vector 
y  onto  the  orthogonal  complement  of  the  vector  subspace  spanned  by  the  columns 
of  A  to  produce  the  adapted  output  z.  By  assumption,  the  auxiliary  array  elements 
receive  negligible  desired-signal  energy  and  are  dominated  by  the  interference.  The 
adaptive  sidelobe  canceller  output  can  be  viewed  geometrically  as  the  projection  of 
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the  array  output  onto  a  subspace  that  is  orthogonal  to  the  subspace  determined  by 
the  interference.  When  the  weak  signal  assumption  is  violated,  the  adaptive  weight 
vector  may  need  to  be  modified  in  accordance  with  linear  constraints  imposed  to 
prevent  desired-signal  suppression  [2]. 

In  order  to  preserve  the  main  features  desired  from  the  antenna  pattern,  while 
maximizing  the  signal-to-interference-plus-noise  ratio  (SINR)  when  there  are  more 
auxiliary  elements  than  interfering  sources,  we  have  found  it  useful  to  use  only  the 
dominant  principal  components  obtained  from  a  singular  value  decomposition 
(SVD)  of  A  in  order  to  determine  the  adaptive  weight  vector  [3].  Let  A  =  US  V  be 
the  singular  value  decomposition  of  A  and  assume  there  are  p  principal  components 
distinguishing  the  p-dimensional  interference  subspace  from  the  M-p  dimensional 
noise  subspace,  as  determined  from  analysis  of  the  singular  values  [3,4],  Then  the 
adapted  output  may  be  expressed  as 

z  =  (I-UUH)y  (3) 

or  as 

zp  =  (I-UpUpH)y  (4) 

if  the  principal  components  method  is  used,  where  Up  denotes  the  first  p  columns 
of  the  matrix  of  left  singular  vectors,  U,  corresponding  to  the  p  dominant  principal 
components. 

It  follows  directly  from  ( 1  >  and  the  SVD  of  A  that  the  adapted  auxiliary  weight 
vector  can  be  determined  as  the  projection  of  the  quiescent  weight  vector  onto  the 


vector  subspace  spanned  by  the  columns  of  V,  which  are  determined  principally  by 
the  interference.  Mathematically, 


w  =  VVHw0 


(5) 


where  w©  is  the  vector  of  quiescent  auxiliary  weights,  or  if  the  principal 
components  method  is  used, 


wp  =  (VpvpH)w0 


(6) 


where  Vp  denotes  the  first  p  columns  of  the  matrix  of  right  singular  vectors  V  (the 
eigenvectors  of  the  sample  covariance  matrix  (AHA)/T),  corresponding  to  the  p 
dominant  principal  components. 

Observe  that  the  principal  components  method  allows  us  to  remove  the 
distinction  between  an  adaptive  sidelobe  canceller  and  a  fully  adaptive  array.  We 
may  use  all  the  array  elements  for  adaptation,  using  the  phase  information  obtained 
from  all  the  elements,  and  choose  the  correct  number  of  degrees  of  freedom  as 
determined  by  the  observed  number  of  dominant  principal  components.  In  this  case 
we  may  equate  A  =  X  and  w0  =  Wq  and  we  see  that  the  adapted  output  can  be 
expressed  as 


zp  =  X(I-VpVpH)wq. 


(7) 


This  allows  us  to  interpret  the  projection  operator,  (I-VpVp^),  not  only  as  an 
operator  used  to  modify  the  quiescent  weight  vector,  but  also  as  a  spatial  filter  that 
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could  be  applied  to  the  array  data  matrix  X  to  remove  the  interference  components 
prior  to  beamforming.  It  is  possible,  therefore,  to  observe  the  effect  of  the  adaptive 
spatial  filter  on  the  elemental  signals. 

2.0  The  Complex  Coherence  Function 

In  the  data  analysis  described  later  we  have  applied  the  principal  components 
method  of  adaptation  to  a  fully  adaptive  array,  operating  at  HF,  for  the  purpose  of 
suppressing  observed  transients.  While  this  provides  an  optimum  performance 
baseline,  we  are  also  interested  in  simplifying  the  computations.  We  will  describe  a 
computationally  simple  algorithm  that,  in  effect,  extracts  the  largest  principal 
component  and  we  will  show  that  it  is  effective  in  suppressing  large  transients  of 
short  duration.  This  algorithm  is  based  on  determining  a  normalized  retrodirective 
steering  vector  vs  from  a  measurement  of  the  complex  coherence  function  of  the 
received  energy  across  the  array,  and  using  the  measured  vector  in  a  projection 
operator  vsv$H  that  can  be  used  to  project  onto  a  one  dimensional  subspace  or  its 
orthogonal  complement  using  the  operator  (I-vsv$H). 

Fundamental  to  all  adaptive  phased  array  processing  is  the  measurement  of  the 
phase  distribution  across  the  array  resulting  from  one  or  more  incident  wavefronts. 
While  this  phase  information  is  contained  in  the  eigenvectors  (or  right  singular 
vectors)  that  we  employ  in  the  principal  components  method,  it  can  also  be 
measured  directly,  particularly  from  a  computation  of  the  complex  coherence 
function  [5].  Let  one  of  the  array  elements  be  used  as  a  reference  and  let  its  analytic 
signal  vector  be  denoted  as  xr  =  where  xr  and  ar  are  complex  vectors  of  T 
samples  and  is  a  scalar  constant  dependent  on  the  spatial  phase  of  the 
wavefront  measured  at  the  reference  sensor.  In  our  notation  the  phasor  e^,  which 
we  assume  to  be  constant  for  the  duration  of  T  samples,  multiplies  each  component 
of  the  complex  vector  ar  which  contains  the  amplitude  and  temporal  phase  variation. 

Similarly  xfc  =  a^e^k  denotes  the  analytic  signal  sampled  at  the  k^1  element.  If  we 


form  the  zero-lag  cross-correlation  between  the  reference  signal  and  each  of  the 
other  elemental  signals  we  obtain  the  complex  coherence  function 


x^xk  =  arHake’^(f)k-<Pr)  k  =  1,...,  N  (8) 


The  real  part  of  (8),  2ar^akcos (cp^-tp^ ,  measures  the  mutual  coherence  across  the 

array  and  the  complex  exponential  e^9k-9r)  describes  the  phase  distribution  across 
the  array.  Notice  that  the  measurement  is  taken  over  T  samples  of  the  data  to  reduce 
the  effects  of  background  noise  that  is  essentially  uncorrelated,  or  phase  incoherent, 
from  one  sensor  to  another.  In  our  analysis  we  have  ignored  the  noise  contribution 
since  we  are  concerned  with  nulling  interference  that  exhibits,  by  implication,  a 
high  per-element  interference-to-noise  ratio. 

The  complex  phasor  e^r-^)  can  be  used  as  a  retrodirective  steering  vector 
to  focus  the  array  to  match,  or  cancel,  the  phase  of  the  impinging  wave  fronts. 

When  there  is  only  one  incident  wavefront,  the  phasor  is  the  steering  vector 
estimate  for  that  source.  A  better  estimate  of  the  incident  phase  can  be  obtained  by 
forming  a  retrodirective  beam  and  using  the  beamformed  output  as  the  reference  in  a 
second  measurement  of  the  complex  coherence  function  intended  to  improve  the 
phase  estimate.  The  process  can  be  iterated  several  times  if  necessary.  Finally,  the 
steering  vector  estimated  after  several  iterations  can  be  used  in  a  one-dimensional 
orthogonal  projection  operator  for  adaptive  nulling  or  spatial  filtering  of  the  array 
data.  The  method  is  somewhat  equivalent  to  extracting  the  largest  principal 
component  of  the  array  data,  but  without  the  necessity  of  computing  a  singular 
value  decomposition. 

When  more  than  a  single  degree  of  freedom  is  required  because  of  multiple 
sources  incident  from  distinctly  different  directions  the  simple  algorithm  described 
will  not  be  adequate,  although  it  may  suppress  a  clearly  dominant  source  of 
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interference.  The  technique  can  be  useful  in  suppressing  short  transient  bursts  of 
interference  when  only  one  source  is  clearly  dominant  in  the  observation  interval. 

3.0  Adaptive  Nulling  of  Transient  Interference 

Figure  1  shows  an  8.5  second  record  (16,384  complex  samples)  of  the 
envelope  of  the  analytic  signal  recorded  from  a  single  element  (number  41)  of  an 
82-element  narrowband  linear  HF  antenna  array  tuned  to  the  30-meter  band. 
Transient  noise  just  before  2  seconds  and  between  6  and  7  seconds  is  clearly 
evident.  The  behavior  is  suggestive  of  that  observed  from  atmospheric  lightning 
discharge  [6],  but  there  is  no  corroborating  evidence  to  firmly  establish  that  these 
transients  are  due  to  lightning  discharge.  Further  analysis  will  show,  however,  that 
the  largest  spikes  are  spatially  coherent  and  can  be  adaptively  nulled  by  the 
procedures  described  above. 


Time  (seconds) 

Figure  1.  Detected  Output  from  Element  41 


We  analyzed  a  64-sample  record  containing  the  large  transient  at  1.9  seconds. 
Figure  2  shows  the  mutual  coherence  measured  across  the  array  using  element  41 
as  a  reference.  It  is  evident  that  the  mutual  coherence  falls  gradually  across  the 
array  and  the  low  spatial  frequency  suggests  a  small  angle  of  incidence  with  respect 
to  boresight  Figure  3  shows  the  phase  variation  across  the  array  (modulo  2k) 
calculated  from  the  complex  coherence  function  after  one  iteration.  The  dashed  line 
plots  the  ideal  phase  distribution,  calculated  for  a  plane  wave  incident  at  4.25 
degrees,  which  closely  matches  the  observed  phase.  In  this  measurement  element  9 
was  inoperative,  producing  only  receiver  noise,  accounting  for  the  anomalous  value 
at  element  9. 


Figure  2.  Mutual  Coherence  (4.25-Degree  Transient) 
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Element  Number 


Figure  3.  Element  Phase  Variation 


Figure  4  shows  the  output  envelope  of  the  signal  from  a  beam  formed  at  4.25 
degrees,  further  supporting  the  conjecture  that  the  transient  is  a  spatially  coherent 
source  incident  at  4.25  degrees.  Closer  analysis  reveals  that  there  are  actually  two, 
and  perhaps  three,  distinct  sources  contained  in  the  record.  Figure  5  is  a  composite 
plot  of  the  beam  scans  formed  from  the  first  three  principal  components  of  the  array 
data;  there  appear  to  be  distinct  sources  at  2.5  and  4  degrees  and  possibly  a  third 
source  at  5  degrees. 
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Figure  4.  Receiver  Output  for  Beam  Formed  at  4.25  Degrees 


Figure  5.  Beam  Scan  for  First  Three  Principal  Components 
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A  similar  analysis  was  performed  using  64  samples  enclosing  the  large 
transient  at  6.5  seconds;  the  mutual  coherence  function  is  shown  in  figure  6. 
Again,  spatial  coherence  is  evident  and  the  spatial  frequency  suggests  a  larger  angle 
of  incidence.  Fine-grained  principal  components  analysis  revealed  a  strong 
transient  at  49.25  degrees  with  an  underlying  narrow-band  interferer  at  46  degrees. 


Figure  6.  Mutual  Coherence  (49.25-Degree  Transient) 
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Let  the  normalized  retrodirective  steering  vector  calculated  for  the  4-degree 
transient  be  denoted  as  V4  and  let  the  retrodirective  vector  for  the  49-degree  source 
be  denoted  as  V49.  Then  a  product  of  orthogonal  projections  can  be  applied  as  a 
spatial  filter  to  the  array  data  matrix  in  order  to  suppress  the  transient  events 


X'  =  X(I-v4V4H)(I-V49V49H).  (9) 

Figure  7  shows  the  spatially  filtered  output  for  element  41.  Although  each  of 
the  components  of  the  spatial  filter  were  determined  from  64  samples,  the  spatial 
filter  was  applied  to  all  16,384  samples  of  the  array  data.  Suppression  of  the 
transient  events  is  evident. 


Figure  7.  Spatially  Filtered  Output  of  Element  4 1 
(Retrodirective  Steering  Vector  Method) 
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Figure  8  shows  the  spatially  filtered  output  of  element  41  when  the  64-sample 
records  surrounding  each  transient  are  subjected  to  a  principal  components  analysis, 
using  a  singular  value  decomposition  of  the  data,  in  order  to  determine  the  adaptive 
spatial  filter  parameters.  In  this  case  five  principal  components  were  used  in 
forming  the  orthogonal  projection  for  each  transient  record.  Comparison  of  figure 
7  with  figure  8  clearly  shows  that  the  simple  algorithm  based  on  estimating  a 
retrodirective  steering  vector  is  suboptimum.  It  does,  however,  produce  significant 
transient  suppression  and  is  computationally  much  less  costly  than  the  SVD- 
dependent  method  of  principal  components. 


Figure  8.  Spatially  Filtered  Output  of  Element  41 
(Principal  Components  Method) 
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4.0  Conclusion 


In  this  paper  we  have  described  a  computationally  simple  algorithm,  useful 
for  adaptive  nulling  of  spatially  coherent  transients  received  by  a  phased  antenna 
array,  based  on  measurement  of  the  complex  coherence  function  to  extract  the  phase 
distribution  from  which  a  retrodirective  steering  vector  is  derived  and  used  in  an 
orthogonal  projection.  The  effectiveness  in  suppressing  transients  was 
demonstrated  by  an  actual  example  using  data  recorded  from  a  narrowband  linear 
HF  array.  The  method  was  compared  with  an  optimal  procedure  based  on  principal 
components  analysis  and  found  to  be  nearly  optimal  for  the  example  used. 
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Abstract 

An  experiment,  is  described,  to  determine  the 
effects  of  mutual  coupling,  near  end  fire,  on  the  array 
gain  of  a  large,  random,  planar  array  of  monopoles.  An 
HF  transmit  antenna  array  2.5  km  in  diameter  with  4,000 
randomly  spaced  elements  is  frequency  scaled  resulting 
in  a  1.8  m  diameter  model  antenna  array.  The 
mathematical  technique  and  c^puter  based  numerical 
calculations  used  to  approximate  the  expected  array  gain 
loss  are  presented.  The  approximate  array  gain  loss  is 
shown  to  be  as  much  as  13  dB  at  0°  elevation  angle. 
Calculation  of  the  actual  array  gain  is  prohibited  by 
the  enormous  amount  of  computer  time  and  memory 
requirements,  for  example  calculation  of  a  single 
element  pattern  using  a  MicroVAX  II  takes  approximately 
47  hours  of  CPU  time.  Therefore,  an  experiment  as 
described  offers  a  practical  approach  to  determine  the 
array  gain  loss. 
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1.  Introduction 


A  large,  (approximately  2.5  km  in  diameter) 
circular,  planar,  transmit  antenna  array  which  operates 
in  the  High  Frequency  (HF)  band  (6-30  MHz)  and  consists 
of  4,000  randomly  spaced  sleeve  dipoles  with  array 
element  electrical  spacing  ranging  from  a  half 
wavelength,  at  the  lowest  frequency  (6  MHz) ,  to  5/2  of  a 
wavelength  at  the  highest  frequency  (30  MHz)  could  be 
useful  for  Over-The-Horizon  (OTH)  radar  applications. 
Such  large  planar  arrays  of  closely  spaced  elements, 
present  several  significant  technical  problems,  when 
phased  to  operate  at  low  elevation  angles  (6el)  near  end 
fire.  One  problem  is  due  to  the  mutual  coupling  between 
antenna  elements  that  can  cause  a  serious  reduction  in 
array  gain.  The  coupling  effects  become  especially 
critical  at  the  low  end  of  the  HF  operating  band  (for 
wideband  arrays)  and  at  very  low  elevation  angles  near 
end  fire.  The  degree  of  mutual  coupling  increases  as 
the  electrical  spacing  of  the  elements  decreases  and 
when  the  elements  are  oriented  such  that  the  element's 
free-space  pattern  exhibits  a  maximum  in  the  direction 
of  its  neighboring  elements.  The  low  HF  frequency  and 
very  low  elevation  angle  capability  is  extremely 
important  for  an  OTH  radar  system.  These  systems  must 
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operate  at  low  elevation  angle  beam  positions  (5°  to  30° 
or  corresponding  polar  angles  85°  to  60°),  to  detect 
targets  at  3000  km  ranges,  and  at  low  HF  frequencies  due 
to  ionospheric  propagation  limitations,  especially  when 
looking  north  and  during  a  significant  portion  of  the 
winter  night. 

The  coupling  effect  is  maximum  at  the  low  HF 
frequency  end  of  the  radar  system's  operating  frequency 
range  where  the  transmit  antenna  elements  are  most 
closely  spaced  electrically.  Because  of  this  potential 
serious  degradation  in  OTH  radar  system  performance,  it 
is  important  that  mutual  coupling  effects,  near  end 
fire,  for  random  planar  arrays  be  identified  and 
quantified. 

2.  Mathematical  Analysis 
2.1  Introduction 

Mutual  coupling  is  the  result  of  the  complex 
interaction  between  antenna  elements  in  electrically 
close  proximity.  The  interaction  between  elements 
changes  the  current  distribution  and  thus  the  element 
pattern  of  an  array  element  from  that  which  it  would 
have  if  it  were  isolated  in  free  space.  The  element 
pattern  of  an  isolated  element  will  in  general  be  quite 
different  from  the  element  pattern  of  an  array  element. 
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Furthermore,  the  element  pattern  of  each  element  in  a 
randomly  spaced  array  will  be  unique.  The  element 
pattern  is  of  interest  because  the  superposition  of  all 
the  array  element  patterns  yields  the  active  array 
radiation  pattern.  [1]  [2] 

The  array  under  consideration  consists  of  4,000 
identical,  parallel,  randomly  located,  matched-loaded, 
quarter  wavelength  monopoles  on  an  infinite  ground 
plane.  This  array  is  a  good  approximation  to  the  4,000 
element  sleeve  dipole  array  described  above.  The 
analytical  method  used  to  analyze  this  array  is  a 
modification  of  the  method  of  moments  used  by  Fenn  [3] 
to  analyze  periodic  arrays  of  monopoles.  The  theory  of 
images  was  used  to  replace  the  array  of  monopoles  over 
an  infinite  ground  plane  by  an  array  of  dipoles  in  free 
space.  The  only  factors  that  must  then  be  accounted  for 
are  that  the  equivalent  dipole  has  twice  the  input 
impedance  and  half  the  gain  of  the  monopole.  This 
dipole  image  approach  was  adopted  because  the  mutual 
impedance  between  parallel  dipoles  in  free  space  can  be 
easily  calculated  using  an  equation  given  by  Castello 
and  Hunk  [ 4 ] . 
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2.3  Mutual  Coupling  Calculations 

Fenn  of  Lincoln  Laboratories  conducted  a  study  to 
estimate  the  mutual  coupling  effect  in  large 
periodically  spaced  arrays  [3],  Franchi  (RADC)  (private 
communication)  using  Fenn's  data  derived  two  empirical 
equations,  one  to  estimate  the  depth  of  the  gain  null  at 
zero  degrees  elevation  angle  as  a  function  of  the  number 
of  array  elements  and  the  element  density  and  the  other 
to  estimate  the  3  dB  width  of  the  null.  The  depth  of 
the  null  is  given  by 

f(Ne,p)  =  1  /  [1  +  p  (0.65  -  1)]  (1) 

and  the  3  dB  width  of  the  null  is  approximated  by 

w  (Ne,p)  =  126 °  [4p/Ne]V4  (2) 

where  p  is  the  element  density  with  respect  to  a  unit 
cell  size  of  >i/2  by  y/2  and  Ne  is  the  total  number  of 
elements.  Predicted  scan  losses  using  this  empirical 
formula  for  an  array  with  Ne  =  4,000  and  p  =  1/2  and  w 
(Ne,p)  «  19 ",  at  6  MHz  are  13.2  dB  and  3  dB  for 
elevation  angles  of  0°  and  9.5°  respectively.  See 
Figure  1  for  a  plot  of  predicted  scan  loss  as  a  function 
of  elevation  angle  and  frequency  for  a  4,000  element 
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periodically  spaced  antenna  array  using  Equations  (1) 
and  (2) .  Note  that  the  scan  loss  increases  as  the 
frequency  and  elevation  angle  decrease. 

3 .  Computer  Based  Numerical  Methods 

3 . 1  Introduction 

The  usefulness  of  computer-based  numerical  methods, 
such  as  Moment  Method,  to  determine  the  mutual  coupling 
effect  in  a  large  array  is  limited  by  the  amount  of 
available  processor  (CPU)  time  and  computer  memory. 
Therefore,  when  an  array  is  very  large,  on  the  order  of 
a  few  thousand  elements,  an  experiment  is  the  only 
practical  method  to  measure  the  expected  reduction  in 
array  gain. 

3.2  Numerical  Method  Limitations 

Although  calculation  of  all  4,000  element  array 
gain  patterns  is  prohibitive,  it  is  useful  to  calculate 
several  element  patterns  to  approximate  the  mutual 
coupling  effect  on  individual  array  elements  and 
calculate  an  average  element  pattern.  To  calculate 
these  element  patterns,  we  must  solve  for  the  current 
vector  (I)  in  Equation  3  either  by  directly  inverting 
the  previously  calculated  impedance  matrix  [Z]  then 
multiplying  it  by  the  applied  voltage  vector  (V)  or  by 
some  other  matrix  solution  technique. 
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(3) 


(i)  =  [zr1  (v) 

The  matrix  inversion  method  is  itself  time  consuming; 
Path  and  Terzuoli  in  1978  [5]  reported  inverting  an 
impedance  matrix  of  (3752  by  3752)  using  the 
Electromagnetic  Surface  Patch  (ESP)  Code  on  a 
VAX-11/780.  The  inversion  took  120  hours  of  CPU  time. 
Rather  than  resorting  to  inverting  the  matrix  to  achieve 
a  matrix  solution,  an  iterative  method,  the  Conjugate 
Gradient  Method,  was  chosen.  This  method  always 
converges  to  the  solution  in  a  finite  number  of 
iterations  given  an  arbitrary  initial  guess,  provided 
the  matrix  is  nonsingular  and  round-off  errors  are  kept 
negligible.  Also,  Conjugate  Gradient  Method  requires 
much  less  memory  («5N)  than  Moment  Method  («N2) .  In 
comparison  to  the  matrix  inversion  method  above  the 
solution  to  Equation  3  where  [Z]  is  a  4,000  by  4,000 
array  would  only  take  approximately  35  hours  of  CPU  time 
on  a  VAX-11/780  [6], [7].  Calculation  of  the  mutual 
impedance  matrix  [Z],  solving  for  the  current  vector  (I) 
and  computing  a  single  phi  cut  element  gain  pattern  for 
a  4,000  element  random  array  such  as  the  pattern  shown 
in  Figure  3  (a)  takes  approximately  46  CPU  hours  on  a 
MicroVAX  II. 
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3.3  Numerical  Results 

3.3.1  Element  Gain  Patterns 

The  element  gain  pattern  of  an  isolated  monopole  on 
an  infinite  ground  plane  is  shown  in  Figure  2  (a)  and 
its  3-D  equivalent  plot  is  shown  in  Figure  2  (b) .  Since 
an  isolated  monopole  pattern  does  not  vary  as  a  function 
of  $,  Figure  2  (b)  is  simply  the  revolution  of  Figure  2 
(a)  about  its  Z  axis.  Note  in  Figure  2  (a)  that  the 
polar  angles  (0)  between  90°  to  70°  and  -90°  to  -70° 
correspond  to  the  near  end  fire,  elevation  angles  (0el) 
between  0°  to  20°.  Figu  3  (a)  is  the  i-0  cut,  element 
gain  pattern  for  element  #  1,  located  near  the  center  of 
the  full  array  of  4,000  elements.  The  gain  pattern  of 
an  array  element  depends  on  the  spacing  between  itself 
and  its  neighboring  elements  and  since  the  element 
spacing  is  random,  the  gain  pattern  varies  as  a  function 
of  the  azimuth  angle  $  as  well  as  the  polar  angle  0. 

This  effect  is  illustrated  in  the  3-D  display  in  Figure 
3  (b) .  Further  calculation  showed  that  there  is  a 
significant  difference  between  gain  patterns  of 
different  array  elements. 

3.3.2  Array  Gain  Loss 

The  array  gain  loss  is  approximated  by  calculating 
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the  difference  between  the  isolated  element  gain  pattern 
and  the  average  element  pattern  of  several  elements 
located  near  the  center  of  the  array. 

A  close  look  at  Figure  4  (a)  will  show  the 
significant  gain  loss  at  the  low  elevation  angles, 
between  the  isolated  element  shown  in  Figure  2  (a)  and 
element  #  1  shown  in  Figure  3  (a) .  Using  the  gain 
patterns  of  several  elements  located  near  the  center  of 
the  array,  an  average  element  pattern  was  calculated  and 
is  shown  in  Figure  4  (b) .  Note  that  the  average  element 
pattern,  unlike  the  individual  element  patterns,  is  not 
a  strong  function  of  $,  therefore,  the  element  pattern 
is  completely  described  by  the  2-D  plot  shown  in  Figure 
4  (b) .  The  array  gain  loss  can  be  approximated  by 
determining  the  difference  between  the  average  element 
pattern  and  the  isolated  element  pattern.  The 
comparison  between  the  isolated  and  average  element 
patterns  is  shown  in  Figure  5  (a)  and  the  array  gain 
loss  as  a  function  of  elevation  angle  (0el)  in  Fi<?ure  5 
(b) .  Since  only  relatively  few  element  patterns  were 
used  to  calculate  the  average  element  pattern  the  actual 
gain  loss  may  be  significantly  different  from  that  shown 
in  Figure  5  (b) .  The  approximation  for  the  gain  loss 
can  be  improved  by  using  an  average  element  pattern 
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based  on  a  larger  number  of  element  patterns,  however, 
calculating  a  large  number  of  element  patterns  is 
impractical,  therefore,  an  experiment  to  measure  the 
actual  array  gain  loss  provides  a  practical  approach. 

4.  Experiment  Design 

4 . 1  Background 

The  physical  dimensions  of  the  transmit  array, 
described  in  the  introduction,  and  the  prohibitive  use 
of  computer  time  and  memory  reguirements  for  calculating 
the  predicted  mutual  coupling  effects,  led  RADC  to 
design  and  conduct  an  experiment  using  a  frequency 
scaled  model.  As  mentioned  previously,  the  effects  of 
mutual  coupling  are  more  serious  at  low  HF  frequencies, 
therefore,  the  model  array  is  scaled  in  frequency  from  6 
MHz  to  8  GHz.  Measurements  where  taken  at  8  GHz. 

4 . 2  Antenna  Array  Model  and  Ground  Screen 

The  antenna  array  model  constructed  is  a  1.875  m 
diameter  array  with  4,000  randomly  spaced  monopole 
elements.  This  is  scaled  in  frequency  for  a  2.5  km 
diameter  array.  The  4,000  monopole  elements  are 
specially  manufactured  Straight  Bulkhead  Jack  Receptacle 
connectors  with  an  extended  center  conductor  equal  to  a 
quarter  wavelength  at  8  GHz  as  shown  in  Figure  6.  The 
monopoles  were  installed  in  the  array  base  plate  and 
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match  terminated  in  50  n  loads  as  shown  in  Figures  7  and 
8  respectively. 

It  is  a  well  documented  fact  that  the  element 
pattern  of  a  single  monopole  over  a  ground  plane  of 
finite  radius  is  affected  by  the  size  of  the  ground 
plane  [8],  Since  we  wish  to  measure  the  array  element 
gain  patterns  for  near  end  fire,  elevation  angles  ee^ 
from  5°  to  20°  the  ground  plane  must  be  sufficiently 
large  such  that  the  isolated  element  gain  pattern 
maximum  occurs  at  polar  angles  ©max  >85°.  Element  gain 
patterns  for  a  thin  quarter  wavelength  monopole  at  the 
center  of  a  circular  ground  plane  of  normalized  radii 
(/3a)  were  calculated  and  the  angles  of  peak  directivity 
plotted.  A  least  square  approximation  was  applied  to 
the  data  to  fit  the  curve 

©max  =  90  -  B  (l/^a)1/2)  (4) 

where  B  is  a  constant.  A  plot  of  the  polar  angle  of 
peak  directivity  versus  /3a  for  B  =  118.505  is  shown  in 
Figure  9.  Using  the  results  illustrated  in  Figure  9  an 
approximate  normalized  ground  plane  radius  of  /3a  «  2,500 
is  required  so  that  ©max  ~  87.5°.  Therefore,  the 
effects  of  the  finite  ground  plane  on  the  element 
pattern  are  negligible  for  ©  <  ©max. 
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Since  the  element  gain  pattern  measurements  will  be 
made  for  one  azimuth  angle,  the  ground  plane  does  not 
need  to  be  circular,  however  the  ground  plane  length  in 
the  azimuthal  direction  of  measurement  must  be 
approximately  equal  to  the  normalized  radius  of  /3a  = 
2,500  and  the  width  sufficiently  large  so  as  to  contain 
several  Fresnel  zones.  Therefore,  to  minimize  the 
effects  of  a  finite  ground  plane  the  model  array  is 
placed  on  a  4.88  m  X  23.16  m  solid  aluminum  ground  plane 
such  that  the  center  of  the  array  is  17.07  m  from  the 
forward  edge  and  6.09  m  from  the  rear  edge.  The  ground 
plane,  shown  in  Figure  10,  is  constructed  from  thirty 
eight,  1.22  m  X  2.44  m  sheets  of  0.64  cm  thick  aluminum 
plates  bolted  together  with  countersunk  bolts  and 
mounted  on  top  of  2.13  m  high  platform.  Electrical 
continuity  between  each  sheet  is  maintained  by  inserting 
an  0.64  cm  X  0.32  cm  aluminum  knitted  wire  mesh  with 
mesh  core  gasket  between  the  plates  and  covering  the 
seam  with  conductive  adhesive  aluminum  tape.  The  ground 
plane  was  leveled  to  a  tolerance  of  ±  0.32  cm  at  each  of 
the  48  adjustment  points  using  a  transit. 

4.3  Experiment  Configuration 

Far  field  measurements  require  that  the  transmit 
and  receive  antennas  be  separated  by  a  minimum  range  of 
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*  =  2  D2/X 

where  D  is  the  diameter  of  the  array  and  R  is  the  range. 
For  D  equal  to  the  diameter  of  the  array  and  a 
wavelength  of  2.5  cm  at  12  GHz  the  range  length  (R) 
required  equals  290  m  An  antenna  range  of  this  size  is 
not  readily  available,  therefore,  an  approximation  for  R 
was  made.  Since  element  patterns  will  be  measured 
individually  with  all  other  elements  match  loaded  the 
coupling  between  elements  should  become  negligible 
within  half  the  diameter  of  the  array.  For  this  case 
then  D  is  equal  to  D/2  and  the  range  R  becomes  R/4  or 
approximately  72  m.  The  elevation  angles  of  interest 
were  determined  by  using  Figure  1  and  noting  that  the 
scan  loss  becomes  negligible  at  elevation  angles  greater 
than  16°.  However,  for  completeness  measurements  should 
be  made  at  as  high  an  elevation  angle  as  the  antenna 
range  permits.  Using  an  existing  RADC  tower  and 
locating  the  array  72  m  from  the  tower,  elevation  angles 
up  to  19 *  can  be  measured.  Figure  11  is  a  sketch  of  the 
Experiment  Configuration.  Since  the  reciprocity  theorem 
applies  for  antenna  radiation  patterns  [9],  the  transmit 
antenna  model  and  its  ground  plane  is  set  up  on  the 
ground  and  operated  in  the  receive  mode  while  a  transmit 
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antenna  horn  is  moved  up  and  down  the  tower.  A  standard 
gain  monopole  is  mounted  on  the  ground  plane  next  to  the 
array  and  4,000  monopole  element  gain  patterns  are  made 
by  taking  the  difference  in  dB  between  the  array  test 
element  and  the  standard  gain  element  received  signal 
levels. 

4.4  Measurement  Equipment  Description 

The  equipment  for  this  experiment  consists  of  the 
transmit,  receive  and  antenna  positioning  equipment 
controlled  by  a  Hewlett  Packard  computer  as  shown  in  the 
block  diagram  illustrated  in  Figure  12. 

4.4.1  Transmit  Equipment 

The  transmit  equipment  consists  of  a  Backward  Wave 
Oscillator  (BWO)  source,  10  W  Traveling  Wave  Tube  (TWT) 
amplifier,  100  ft  of  low  loss  coaxial  cable  and  X-Band 
horn  antenna.  The  X-band  horn  antenna  has  elevation  and 
azimuthal  patterns  such  that  the  gain  (towards  the 
receive  element  under  test)  variation  is  less  than  1  dB 
as  the  horn  travels  the  height  of  the  tower. 

4.4.2  Receive  Equipment 

The  receive  equipment  consists  of  an  SA  1711  receiver, 
ratiometer,  RF  switch,  and  two  mixers.  The  receive 
equipment  is  set  up  in  the  shared  mixer  mode 
configuration  to  provide  synchronous  switching  in  order 
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to  time  share  the  mixer  between  the  standard  gain 
monopole  and  the  array  test  monopole  element  while 
maintaining  phase  lock  of  the  RF  signal  through  the 
Automatic  Phase  Control  (APC)  mixer.  Receiver  channel  A 
is  connected  to  the  standard  gain  monopole  and  channel  B 
to  the  array  element  under  test.  The  ratiometer 
displays  the  difference  in  dB  between  these  two 
channels.  Ten  samples  are  averaged  per  guarter  degree 
of  elevation  angle  and  this  average  is  recorded  on  disk 
for  later  processing.  Signal  averaging  is  used  to 
minimize  the  local  background  environmental  noise 
effects. 

4.4.3  Antenna  Positioning  Equipment 

The  antenna  positioning  equipment  consists  of  a 
controller  card,  microstep  indexer  card,  driver  card  and 
stepper  motor  configured  as  shown  in  Figure  13.  The 
system  controller  card  is  operated  through  the  IEEE-488 
bus  by  a  program  run  on  the  HP  computer.  Signals  from 
the  controller  are  sent  to  a  microstep  indexer  which 
provides  the  driver  card  with  the  correct  number  of 
properly  sequenced  phase  control  signals  to  accelerate, 
drive  and  decelerate  the  stepper  motor  as  well  as 
correct  for  backlash  in  the  10:1  gear  reducer.  The 
driver  card  supplies  power  to  a  2-phase  stepper  motor  to 
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turn  a  10:1  gear  reducer  and  drive  the  drum  and  cable 
assembly  that  is  connected  to  the  X-band  horn  cart.  On 
command  from  the  HP  computer  the  stepper  motor  drives 
the  transmitting  X-band  horn  cart  along  the  rail  system 
attached  to  the  tower  with  a  positional  accuracy  of 
greater  than  0.08  cm. 

5.  Element  Gain  Pattern  Measurements 

5.1  Reference  Monopole  Element  Gain  Pattern 

The  reference  monopole  element  gain  pattern  was 
measured  in  the  conventional  way  using  a  standard  gain 
horn  and  the  gain  transfer  measurement  technique. 

The  reference  monopole  element  is  identical  to  the  4,000 
elements  in  the  array  and  is  mounted  on  the  ground  plane 
1.73  m  (46  wavelengths)  away  from  the  center  of  the 
array.  The  standard  gain  horn  is  located  0.5  m  from  the 
reference  monopole  element  next  to  the  ground  plane  and 
at  the  same  distance  from  the  transmitting  horn  as  the 
reference  element.  The  reference  element  and  the 
standard  gain  horn  are  used  in  the  receive  mode  with  the 
transmitting  horn  located  on  the  tower  at  an  elevation 
angle  of  10'  with  respect  to  the  ground  plane.  Then  a 
complete  reference  element  pattern  measurement  is  made 
from  O'  to  19"  at  1/4°  increments.  Transmitter  power 
fluctuations  due  to  the  10  Watt  TWT  output  and  cable 
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flexing  as  the  transmitting  horn  travels  up  the  tower 
are  reduced  by  averaging  100  measurements  at  every  1/4* 
increment.  Also,  a  correction  is  made  to  the  reference 

t 

element  pattern  measurement  to  account  for  the  gain  loss 
as  the  transmitting  horn  travels  straight  up  the  tower 
rather  than  at  a  constant  distance  from  the  reference 
element.  The  resulting  reference  monopole  element  gain 
pattern  is  shown  in  Figure  11  (a)  compared  to  the 
calculated  monopole  gain  pattern  of  a  monopole  on  an 
infinite  ground  plane.  Notice  the  peak  of  the  measured 
element  pattern  is  below  0el  =  2.5”,  therefore, 
measurements  made  at  elevation  angles  greater  than  2.5” 
are  only  effected  by  the  ±  1/2  dB  ripple  on  the  element 
pattern.  The  reference  monopole  element  pattern  is  not 
greatly  effected  by  the  finite  size  of  the  ground  plane 
or  the  presence  of  the  4,000  element  array. 

t  5.2  Array  Monopole  Element  Gain  Patterns 

The  difference  in  dB  between  an  array  element  and  the 

■*  reference  element  received  signals  are  measured 

simultaneously  and  averaged  over  ten  measurements  for 
each  1/4”  of  elevation  angle  from  0”  to  19°.  Averaging 
is  used  to  minimize  the  local  background  environment 
noise  effects.  Then,  the  difference  must  be  added  to 
the  previously  measured  reference  element  gain  pattern 
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to  produce  an  array  element  gain  pattern  in  dBi.  The 
above  procedure  is  repeated  4,000  times,  manually 
switching  elements  between  measurements.  A  comparison 
between  the  measured  array  element  gain  pattern  for 
element  #  1  located  near  the  center  of  the  array  and 
calculated  gain  for  element  #  1  is  shown  in  Figure  14 
(b)  .  The  difference  between  the  measured  and  the 
calculated  pattern  may  be  due  to  several  approximations 
in  the  theory.  Assuming  a  single  sinusoidal  current 
distribution  on  the  wire  is  probably  the  largest 
contributor  to  the  error,  especially  since  the  actual 
elements  have  a  length  to  diameter  ratio  (1/d)  of  7.38. 
The  delta-gap  voltage  source  model  used  for  the  feed 
region  may  also  have  some  effect.  Also,  the  ripple  on 
the  measured  element  pattern  is  due  to  the  finite  size 
ground  plane  where  the  theory  assumes  an  infinite  ground 
plane. 

5.3  Array  Gain  Pattern 

The  array  gain  is  calculated  as  the  superposition  of 
all  4,000  array  element  patterns.  The  array  gain 
patterns  calculated  without  mutual  coupling,  predicted 
with  mutual  coupling  and  measured  with  mutual  coupling 
for  ♦  =  O’  and  5°  <  ee^  <  19 0  are  shown  in  Figure  15. 

In  this  figure,  the  calculated  array  gain  without  mutual 
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coupling  effects  is  the  top  curve.  The  middle  curve  is 
based  on  the  predicted  gain  loss  due  to  mutual  coupling 
given  by  equations  (1)  and  (2)  for  a  periodically  spaced 
array  with  Ne  =  4,000  and  p  =  1/2.  The  bottom  curve  is 
the  measured  array  gain  pattern  including  mutual 
coupling.  Note  that  the  measured  array  gain  (bottom 
curve)  is  lower  than  predicted  for  a  periodically  spaced 
array  (middle  curve) .  The  array  gain  loss  between  the 
calculated  without  mutual  coupling  (top)  curve  and  the 
measured  with  mutual  coupling  (bottom)  curve  at  ©el  =  5° 
is  approximately  8  dB  and  is  still  as  much  as  2  dB  at 
0ei  =  19- . 

6 .  Summary 

The  analytical  technique  used  to  analyze  the 
expected  performance  of  a  circular,  randomly  spaced 
array  of  monopoles  over  an  infinite  ground  plane  was 
presented.  This  analysis  is  valid  for  antenna  elements 
with  sinusiodal  current  distributions  such  as  monopoles 
over  an  infinite  ground  plane.  Excessive  CPU  time  and 
memory  requirements  prohibited  using  this  method  to 
calculate  the  expected  array  gain  based  on  4,000  element 
patterns.  However,  an  approximate  array  gain  loss  was 
calculated  using  an  average  element  pattern  based  on 
several  elements  located  near  the  center  of  the  array. 
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Due  to  the  computational  difficulty  of  calculating  the 
array  gain  and  the  physical  dimensions  of  the  proposed 
transmit  antenna  array,  RADC  designed  an  experiment 
using  a  scaled  antenna  array  model.  The  antenna  array 
model  and  ground  plane  were  constructed  at  the  RADC 
Hanscom  AFB  antenna  test  range.  The  measurement  of  all 
4,000  transmit  antenna  array  element  patterns  was 
completed  and  the  array  gain  pattern  calculated.  The 
result  of  the  measurements  was  compared  to  the  predicted 
array  gain  and  the  calculated  array  gain  without  mutual 
coupling.  The  array  gain  loss  between  the  calculated 
without  mutual  coupling  and  the  measured  with  mutual 
coupling  is  approximately  8  dB  at  0el  =  5°  and  2  dB  at 
©ei  =  19°.  Therefore  the  array  gain  loss,  near  end 
fire,  due  to  mutual  coupling  must  be  taken  into  account 
when  designing  the  transmit  antenna  array.  Future 
efforts  to  reduce  the  gain  loss  at  low  elevation  angles 
will  include  low  coupling  array  element  designs  and 
terminating  array  elements  with  the  impedance  of  an 
average  element's  input  impedance  when  scanned  near  end 
fire. 
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Figure  1.  Predicted  Scan  Loss  in  dB  for  a  4,000  Element 
Periodically  Spaced  Array  as  a  Function  of  Frequency  and 
Elevation  Angle. 
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Figure  2.  (a)  Element  Gain  Pattern  for  an  Isolated 

Monopole  on  an  Infinite  Ground  Plane.  (b)  3-D  view. 
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Figure  3.  (a)  Gain  Pattern  for  Element  #  1  in  the  4,000 

Element  Array  Environment.  (b)  3-D  view. 
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Figure  4.  (a)  Gain  Comparison  Element  Between  an 

Isolated  Monopole  and  a  Monopole  Near  the  Center  of  the 
Array.  (b)  Average  Element  Pattern. 
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Figure  5.  (a)  Comparison  Between  Isolated  Element 

Pattern  and  Average  Element  Pattern.  (b)  Array  Gain 
Loss  as  a  Function  of  Elevation  Angle  (©ei) • 


Figure  6.  Monopole  Element. 
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Figure  7.  Top  View  of  AOTH  Transmit  Antenna  Array  Model 
with  4,000  Elements  Installed. 


Figure  8.  Bottom  View  of  the  Transmit  Antenna  Array 
Model  with  4,000  Loads  Attached. 
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NORMALIZED  GROUND  PLANE  RADIUS 


Figure  9.  Least  Square  Approximation  of  Polar  Angle  of 
Peak  Directivity  Versus  Normalized  Ground  Plane  Radius 
for  15  <  0a  <  3600. 


I 


3 


Figure  10.  Ground  Plane  with  the  Model  AOTH  Transmit 
Antenna  Array  Installed. 
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ANTENNA  MODEL  SPECIFICATIONS 


DIAMETER 
ELEMENTS 
ELEMENT  SPACING 
FREQUENCY  RANGE 
GROUND  PLANE 


1.875  m 

4000,  1/4  WAVELENGTH  MONOPOLES 
RANDOM.  6/2  WAVELENGTH  EXCLUSION  ZONES 
8-12  GHz  (6-9  MHz) 

4.88  -  23.16  m 


Figure  11,  Experiment  Configuration 
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Figure  13.  Antenna  Positioning  Equipment  Block  Diagram. 
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Figure  14.  (a)  Gain  Comparison  Between  a  Monopole  on  an 

Ifinite  Ground  Plane  and  the  Reference  Monopole.  (b) 
Comparison  Between  the  Measured  and  the  Calculated  Gain 
for  Element  #  1. 
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ANTENNA  DESIGNS  FOR  THE 
AN/FPS-118  OTH  RADAR 

K.  JOHN  SCOTT 

GENERAL  ELECTRIC  COMPANY 
SYRACUSE,  NEW  YORK 


ABSTRACT 

The  AN/FPS-1 18  Radar  is  a  Frequency  Modulation/  Continuous  Wave  (FM/CW) 
bi-static  system  having  a  transmit  site,  a  receive  site,  and  an  operations  site.  The  transmit 
and  receive  sites  are  generally  separated  by  approximately  100  miles  to  provide  sufficient 
isolation  between  the  transmit  and  receive  functions. 

The  transmit  subsystem  utilizes  elemental  transmitters  and  a  low-level 
beamformer  which  has  amplitude  and  phase  control  for  each  transmitter.  The  elemental 
transmitters  drive  a  twelve-element  dipole  array  for  each  of  six  bands  with  a  backscreen 
and  groundscreen  reflector.  A  calibration  coupler  is  provided  at  each  transmitter  output  to 
sample  amplitude  and  phase.  An  iterative  calibration  procedure  is  used  to  establish  the 
aperture  illumination  needed  for  the  desired  beamwidth,  scan  angle,  and  frequency. 
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The  receive  subsystem  utilizes  elemental  receivers  and  a  digital  beamfonmer.  The 
elemental  receivers  are  connected  to  a  monopole  array  of  246  elements  for  three 
interleaved  bands  with  a  backscreen  and  groundscreen  reflector.  A  test  signal  is  inserted 
at  each  receiver  input  to  provide  amplitude  and  phase  corrections  for  each  receiver  path. 

An  ionospheric  backscatter  sounder  system  provides  the  data  needed  to  direct  the 
operation  of  the  radar.  The  sounder  uses  a  transmit  antenna  having  a  beamwidth 
nominally  equal  to  the  60-degree  scan  coverage  of  each  radar  segment.  Eight  narrow 
receive  beams,  spaced  in  7.5-degree  increments,  are  provided  by  an  analog  beamformer 
which  utilizes  a  section  of  the  radar  receive  antenna.  This  paper  describes  the  transmit  and 
receive  antenna  designs  which  evolved  from  the  system  performance  requirements. 

1.0.  OVERVIEW 

The  site  ground  conditions  and  the  surrounding  terrain  are  integral  parts  of  a  High 
Frequency  (HF)  antenna  design.  At  HF  0«ng  wavelengths),  antennas  are  physically  very 
large  and  site  groundslope,  blockage  in  the  "look"  direction,  and  soil  characteristics 
(dielectric  constani/conductivity)  must  be  included  in  antenna  performance  tradeoffs.  The 
choice  of  polarization  is  dictated  by  the  elevation  angle  coverage  requirements  and  antenna 
height  constraints.  The  AN/FPS-1 18  antennas  use  primarily  vertical  polarization  and  long 
groundscrecos  to  enhance  pattern  response  at  low  elevation  angles. 

The  CONUS  OTH  System  will  consist  of  an  East  Coast  Radar  (ECRS),  A  West 
Coast  Radar  (WCRS),  an  Alaskan  Radar  (ARS),  and  a  Central  Radar  System  (CRS). 
Figure  1  shows  the  total  coverage  of  the  CONUS  OTH  System. 
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Figure  1.  ANIFPS-118  Radar  System  Coverage 


Six  antennas  are  provided  for  the  ECRS.  Three  broadside  transmit  arrays  and 
three  broadside  receive  arrays  located  in  the  state  of  Maine  are  used  to  provide  180-degree 
coverage  along  the  Atlantic  Coast  of  the  North  American  continent  Figures  2  and  3  show 
the  transmit  and  receive  antennas,  respectively.  A  transmit/receive  antenna  pair  covers  a 
60-degree  segment  in  eight  sequential  7.5-degree  increments.  The  transmit  beam  for  any 
increment  is  overlaid  with  three  2.5-degree  receive  beams. 


Figure  3.  Receive  Antenna 
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A  transmit  antenna  array  consists  of  six  sub-bands  designated  A  through  F  in 
ascending  frequency  order.  Each  sub-band  is  a  12-element  broadside  linear  array  of 
dipoles,  mounted  in  front  of  a  backscreen  and  above  a  groundplane.  The  frequency 
bandwidth  ratio  of  each  sub-band  is  1.348,  providing  a  5-to  30-MHz  total  bandwidth. 

Bands  A  through  D  have  the  radiating  element  canted  at  a  45-  degree  angle  relative 
to  the  vertical.  This  angle  is  chosen  to  provide  coverage  at  both  high  and  low  elevation 
angles  using  an  elliptical  polarization,  which  is  primarily  vertical  at  low  angles  and 
horizontal  at  high  angles.  Bands  E  &  F  use  a  vertically-oriented/polarized  element  to 
provide  only  low  angle  coverage.  Subarray  placement,  element  cant  angle,  and  inter¬ 
subarray  spacings  are  chosen  to  minimize  coupling  between  subarrays  and  minimize 
transmission  line  loss.  This  results  in  a  band  A,  C,  E,  D,  F,  B  arrangement  of  subarrays. 

The  receive  antenna  is  a  broadside  linear  array  of  246  equally  spaced  monopole 
elements  placed  in  front  of  a  backscreen  and  driven  against  a  groundplane.  Each  of  the 
246  elements  covers  a  5-  to  30-MHz  frequency  band.  The  antenna  array  is  used  in  three 
bands  designated  as  high,  mid,  and  low  band.  Each  band  has  82  active  elements,  with  the 
high  band  using  the  82  contiguous  centrally  located  elements;  mid-band  using  the  82 
centrally  located,  even  numbered  elements;  and  the  low  band  utilizing  every  third  element 
of  the  entire  array.  Table  1  shows  the  overall  dimensions  of  the  AN/FPS-118  Antennas 
and  Figures  4  and  5  show  detailed  dimensions  of  the  transmit  and  receive  antennas, 
respectively. 
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TABLE  1. 

ANTENNA  CHARACTERISTICS 
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12 
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Each  segment  has  an  in-coverage  backscatter  sounding  antenna  system.  Two 
additional  out-of-coverage  sounder  antennas  are  provided  for  the  East  Coast  Radar.  One 
is  directed  due  west  and  the  other  pointed  at  geomagnetic  north.  The  sounder  transmit 
antenna  is  a  two-element  array  of  Log  Periodic  elements.  The  in-coverage  antennas  have  a 
750-foot  long  groundscreen  and  the  out-of-coverage  antennas  have  only  enough 
groundscrcen  to  stabilize  the  input  impedance.  The  in-coverage  receive  sounding  beams 
are  formed  using  24  high,  24  mid,  and  24  low  band  radar  elements.  The  element  signals 
are  obtained  by  a  3-dB  power  split  after  the  first  amplification  stage  in  the  receivers.  An 
analog  beamformer  for  each  group  of  24  elements  provides  eight  beams  steered  in  7.5- 
degree  increments  to  fixed  angular  positions  which  covers  a  60-degree  azimuth  angle. 
Full  time  delay  is  used  for  steering  to  obtain  frequency  independent  beam  pointing 
directions.  A  switching  matrix  provides  two  simultaneous  beams  from  any  of  the  three 
bands  which  can  be  routed  to  two  sounder  receivers.  An  out-of-coverage  sounder  receive 
antenna  can  also  be  connected  to  one  of  the  sounder  receivers.  Each  of  these  antennas 
consists  of  a  high  and  low  band  with  each  band  having  ten  elements.  The  elements  are 
positioned  on  the  backside  of  the  radar  antenna  backscreen  and  they  use  the  radar 
backscreen  as  a  reflector.  An  analog  beamformer  is  used  to  collimate  a  beam  in  the 
desired  direction,  depending  on  the  pointing  direction  of  the  companion  radar  antenna. 
2.0.  ARRAY  DESIGN 

The  total  5-  to  30-MHz  frequency  band  is  divided  into  six  sub-bands  for  the 
transmit  antenna.  The  frequency  bandwidth  ratio  of  each  band  is  1.348  (fh/fl)  The 
azimuth  clement  spacing  is  0.62  X  at  the  high  end  of  each  band,  and  each  band  has  twelve 
elements.  This  spacing  is  selected  to  place  the  azimuth  grating  lobe  oeak  well  into 
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imaginary  space  at  a  maximum  scan  of  27  degrees.  This  is  necessary  with  a  short  aperture 
(wide  grating  lobe)  to  minimize  the  grating  lobe  impact  on  average  sidelobe  levels.  The 
broadside  azimuth  beamwidth  for  uniform  illumination  varies  from  7.0  to  9.6  degrees 
over  the  band.  A  maximum  scan  angle  of  27  degrees  is  chosen  to  place  the  outer  -3-dB 
point  of  the  beam  beyond  the  30-degree  coverage  needed  at  the  high  band  end  of  each 
subarray. 

The  element  height  above  the  groundplane  for  bands  A  through  D  is  0.41  X  at  the 
high  end  of  each  band.  The  height  for  bands  E  &  F  is  0.44  X  at  the  high  end  of  each 
band.  These  heights  and  the  element  orientation  are  selected  to  provide  the  required 
elevation  angular  coverage. 

The  receive  antenna  element  spacing  is  0.61  X  (20  feet)  at  30  MHz  and  is  selected 
to  place  the  first  null  of  the  grating  lobe  at  the  edge  of  imaginary  space  for  a  maximum 
scan  of  36  degrees.  This  maximum  scan  angle  is  chosen  to  permit  depointing  the  receive 
antenna  6  degrees  relative  to  the  pointing  angle  of  the  transmit  antenna,  while  still 
providing  a  60-degrce  angular  coverage  from  each  antenna  pair.  The  reasons  for 
depointing  the  antennas  are:  1)  to  correct  for  parallax  arising  due  to  the  transmit  and 
receive  antenna  physical  separation  and,  2)  to  ensure  that  the  transmit  and  receive 
backiobcs  arc  never  in  alignment.  This  improves  the  round-trip  backlobe  level  and 
provides  backlobe  auroral  clutter  suppression.  The  antenna  has  246  total  elements,  of 
which  only  82  are  used  at  any  given  time.  The  high  band  uses  82  contiguous  elements  at 
the  center  of  the  array;  in  the  mid  band,  82  even-numbered  central  elements  are  employed 
and,  in  the  low  band,  every  third  clement  is  used.  The  frequency  bands  are  5.0  to  9.7 
MHz,  9.7  to  15.0  MHz  and  15.0  to  30.0  MHz.  The  azimuth  beamwidth  is  held  to  2.5 
degrees  across  the  frequency  band  (when  not  aperture  limited)  by  adjusting  the  number  of 
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elements  used  and  the  Taylor  weighting  on  the  aperture.  At  the  lower  end  of  the  high  and 
low  bands  the  beamwidth  increases  to  a  maximum  of  3  degrees.  This  weighting,  along 
with  the  error  sidelobes,  provides  the  required  sidelobe  suppression  to  allow  operation  in 
an  auroral  environment.  For  higher  azimuth  resolution  the  aperture  weight  can  be  changed 
to  provide  a  narrower  beam. 

3.0.  BACKSCREEN  AND  GROUNDSCREEN  DESIGN 

The  backscreen  for  the  antennas  is  constructed  of  copper-coated  steel  wire.  The 
wire  diameter  and  spacing  is  selected  to  reduce  the  leakage  through  the  screen  to  a  level 
which  is  consistent  with  the  required  backlobe  level.  The  backscreen  height  is  chosen  to 
control  backlobe  levels  due  to  diffraction  at  the  screen  truncation  edge.  Connected  mesh 
construction  is  used  to  control  the  vibrational  characteristics  of  the  backscreen  rather  than 
unconnected  wires  which  vibrate  at  high  frequencies  and  may  affect  the  subclutter 
visibility  of  the  radar.  Three  mesh  openings  are  used  for  the  transmit  antenna.  They  are: 
10"  by  10"  for  band  F;  14"  by  14"  for  bands  D  and  E;  and  22"  by  22"  for  bands  A,  B, 
and  C.  The  receive  antenna  uses  a  20"  by  10”  mesh  opening  with  the  10"  spacing  in  the 
horizontal  direction.  A  galvanized  steel  structure  consisting  of  towers  and  trusses 
supports  the  backscreen.  The  backscreen  material  is  fabricated  in  rolls  of  various  widths. 
The  rolls  of  material  are  thermite-welded  together  at  the  site  into  panels.  These  are  then 
tensioned  and  thermite-welded  to  the  support  structure. 

The  steel  structures  are  fabricated  into  the  largest  fully  welded  sections  which  can 
be  galvanized  and  shipped  to  site.  During  assembly,  all  metal-to-mctal  joints  in  the 
structure  are  coated  with  a  conducting,  anti-corrosive  material  to  prevent  any  noise  and/or 
spurious  frequency  generation  from  the  joints  over  the  life  of  the  system.  The  structures 
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arc  galvanized  to  minimize  non-linear  effects  due  to  the  iron  content  of  the  steel,  and  to 
prevent  corrosion. 

The  groundscreen  for  the  antennas  is  laid  on  the  ground  to  provide  a  reflecting 
surface  to  control  elevation  beamshape  and  ground  losses.  The  screen  material  is  an 
aluminum-coated  steel  wire  which  is  compatible  with  normal  soils.  The  groundscreen  is 
constructed  of  14  gauge  wire  welded  into  a  6"  by  6"  mesh.  The  material  is  fabricated  in  6- 
foot  wide  rolls,  380  feet  in  length.  The  rolls  are  installed  in  strips,  normal  to  the  antenna 
line.  The  strips  are  spaced  8"  apart,  and  are  spot-welded  together  every  4  feet  for  the 
transmit  antenna  and  every  3  feet  for  the  receive  antenna. 

A  computer  program  was  developed  to  calculate  the  elevation  pattern  performance 
of  the  antennas.  The  analysis  includes  the  diffraction  from  the  backscreen  top  edge,  the 
groundscreen  front  edge,  the  leakage  through  the  screen,  and  the  effects  of  the  ground 
surrounding  the  antenna.  A  flat  earth  is  assumed  for  the  analysis.  Good  correlation  has 
been  obtained  between  measured  scale  model  patterns  and  patterns  calculated  for  the  scale 
model  using  the  computer  program.  Measured  azimuth  patterns  for  the  ECRS  transmit 
antennas,  taken  at  several  elevation  angles,  show  good  correlation  with  calculated  data. 
The  program  integrates  the  patterns  to  calculate  the  directive  gain.  Figures  6  through  8 
show  the  transmit  elevation  coverage,  receive  elevation  coverage,  and  the  antenna  directive 
gain,  respectively. 
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Figure  6.  Transmit  Antenna  Elevation  Coverage 
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Figure  7.  Receive  Antenna  Elevation  Coverage 
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Figure  8.  Antenna  Directive  Gain 

4.0.  ELEMENT  DESIGN 

The  transmit  element  must  provide  a  good  impedance  match  (low  VSWR)  to  the 
transmitter  to  obtain  efficient  power  transfer  to  space.  Each  array  element  has  a  broadband 
transmitter  which  supplies  80  kW  of  power  for  VSWRs  of  less  than  1.6.  At  a  VSWR  of 
2.0,  the  power  output  falls  to  approximately  40  kW.  To  provide  low  VSWRs  over  the 
frequency  band  and  scan  angle,  a  broad  band  dipole  is  chosen  for  the  transmit  element. 

The  dipole  design  uses  a  bazooka  balun  to  provide  the  transition  from  a  50  ohm, 
6-1/8"  diameter  coaxial  feed  line  to  the  balanced  current  on  the  dipole  arms.  A  square 
aluminum  truss  structure  provides  the  element  support  and  also  forms  the  X/4  shorted  line 
section  for  the  balun.  The  radiating  arms  of  the  dipole  are  triangular  aluminum  truss 


structures,  which  are  isolated  from  the  aluminum  support  by  a  dielectric  structure.  The  6- 
1/8"  coaxial  input  line  proceeds  through  the  aluminum  support  to  the  element  feed  point. 
Each  arm  is  fed  through  an  open  circuited  3-1/8"  coaxial  line  located  in  the  arm.  This 
provides  a  series  capacitance,  which  is  used  in  conjunction  with  a  two-section  impedance 
transformer  in  the  6-1/8"  input  line,  to  match  the  active  dipole  impedance  to  50  ohms. 
Figure  9  shows  dimensional  data  and  the  configuration  of  the  transmit  elements. 
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Figure  9.  Transmit  Element  Dimensional  Data 

Scale  models  of  the  elements  were  constructed  without  the  impedance 
transformers  to  measure  the  single  clement  feed  point  impedance.  Seven  element  scale 
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models  were  constructed  of  bands  C  and  E  to  measure  the  mutual  coupling  coefficients  in 
the  array  environment  for  the  two  different  polarizations.  A  computer  program  was 
developed  to  calculate  from  the  measured  data  the  active  input  impedance  of  each  element 
versus  scan  angle  and  frequency.  An  impedance  matching  program,  using  an  iterative 
procedure,  calculates  the  parameters  of  a  two-section  transformer  which  minimizes  the 
input  VSWR  for  all  scan  angles,  frequencies,  and  element  positions  in  each  subarray.  The 
band  F  VSWR  is  optimized  from  22.25  to  28  MHz  due  to  operational  frequency 
restrictions  from  28  to  30  MHz.  The  scale  model  VSWR  predictions  are  shown  in  Table 
2.  Full-scale  VSWR  measurements  compare  very  well  with  the  scale  model  predictions. 
Figure  10  shows  a  comparison  of  measured  and  calculated  data  at  broadside  scan  and  10.6 
MHz  for  the  band  C  sub-array.  Figure  11  shows  the  F  band  subarray  and  construction 
details  of  the  transmit  elements. 

Each  element  in  bands  A  through  E  is  cantilevered  from  a  backscreen  support 
tower.  For  band  F,  three  elements  are  cantilevered  from  a  truss  which  is  supported  by 
towers.  The  element  6  1/8"  transmission  line  is  brought  down  the  rear  of  the  towers  to 
approximately  7  feet  from  the  ground.  A  bi-directional  coupler  is  provided  for  monitoring 
element  performance. 


TABLE  2 

TRANSMIT  ANTENNA  ACTIVE  VSWR 
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Figure  11.  Band  F  Transmit  Elements 


124 


The  receive  antenna  element  is  a  17.8-foot  long  monopole  fabricated  from 
galvanized  steel.  The  element  has  an  18"  by  18"  by  18”  triangular  cross  section  and  is 
isolated  from  the  ground  by  three  fiberglass  dielectric  insulators.  The  element  is  fed 
against  the  groundscreen  using  an  end  terminal  assembly  which  has  a  "Type  N"  coaxial 
input  connector.  A  surge  protector  is  used  to  prevent  receiver  damage  if  lightning  strikes 
the  antenna  structure.  Figure  12  shows  the  receive  element. 


Figure  12.  Receive  Antenna  Element 
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The  receive  element  must  operate  over  the  total  5-  to  30  MHz  frequency  band. 
The  element  length  and  distance  from  the  backscreen  are  chosen  to  be  as  large  as  can  be 
used  without  developing  an  appreciable  reduction  in  the  element  gain  in  the  broadside 
direction.  This  results  in  an  electrical  element  height  of  0.54  X  above  the  ground  plane 
and  a  0.28  X  (9.2  feet)  electrical  distance  to  the  backscreen  at  30  MHz.  These  dimensions 


provide  an  acceptable  loss  for  operation  at  5  MHz.  The  return  loss  of  a  typical  element  in 
a  loaded  array  environment  is  shown  in  Figure  13.  This  variation  in  return  loss  and  the 
resulting  mismatch  loss  is  compatible  with  the  external  noise  characteristics  at  HF 


frequencies. 
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Figure  13  Receive  Element  Return  Loss 


The  complete  array  is  filled  with  elements  to  maintain  impedance  similarity  for  all 
elements  in  the  three  bands.  Unused  elements  are  terminated  in  50-ohm  loads.  This  is 
necessary  to  meet  element  similarity  error  requirements  for  sidelobe  control. 

5.0.  FEED  SYSTEM 

The  transmit  antenna  feed  system  is  designed  for  an  input  power  of  100  kW  CW 
and  a  VSWR  of  2  for  each  element.  Semi-flexible  50-ohm  coaxial  transmission  lines  are 
used  to  connect  the  transmitter  outputs  to  the  element  inputs.  Bands  A  &  B  use  3" 
diameter  line;  bands  C  and  D  use  4”  diameter  line,  while  bands  E  &  F  use  5”  diameter 
line.  An  adapter  is  used  to  interface  with  the  6-1/8"  rigid  line  element  input.  A  total  of 
61,000  feet  of  coaxial  line  is  used  for  the  transmit  antenna.  Each  line  is  long  enough  to 
extend  from  the  element  to  its  transmitter  terminal.  The  electrical  length  of  the  lines  are 
measured  during  installation  and  the  data  used  in  a  beamforming  and  steering  algorithm  to 
collimate  the  beam  at  the  selected  scan  angle. 

The  receive  feed  lines  are  commercially  available  0.412”  diameter  coaxial  cables 
which  are  used  in  the  cable  T.V.  industry.  The  cable  characteristic  impedance  is  75  ohms 
and  164  cables  are  used  to  connect  the  active  elements  to  the  receivers.  Equal  length 
cables  are  used  to  all  the  active  elements  which  are  enclosed  by  the  mid-band  aperture. 
The  outboard  14  low-band  elements  on  each  end  of  the  antenna  use  just  enough  cable  to 
connect  the  elements  to  the  receiver.  The  cables  are  installed  and  buried  in  a  trench,  5-foot 
deep,  which  provides  a  maximum  temperature  variation  of  1 1  degrees  C  year  round  and 
limits  phase  changes  versus  temperature  to  an  acceptable  level.  The  cable  lengths  are 
measured  and  trimmed  to  a  tolerance  of  3.2  degrees.  These  measured  phase  differences 
are  then  used  as  a  phase  correction  in  the  beamforming  and  steering  algorithm  in  the 
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receive  digital  beamformer.  A  total  of  316,000  feet  of  cable  is  used  for  the  receive 
antenna. 

6.0.  Structural  Design 

The  antennas  are  designed  to  survive  a  peak  wind  velocity  of  90  knots  at  a  30- 
foot  elevation  without  ice  loading.  With  2  inches  of  rime  ice  or  1 .5  inches  of  clear  ice  on 
major  support  structures,  the  antennas  must  survive  a  50-knot  peak  wind,  acting  at  a  30- 
foot  elevation.  The  backscreen  and  elements  are  not  considered  major  support  structures, 
and  0.5  inches  of  clear  ice  is  used  as  the  ice  loading  requirement  on  these  structures. 

A  finite  element  model  of  the  transmit  antenna  is  created  for  each  band  using  the 
ANSYS  Engineering  Analysis  System.  The  structural  members  are  represented  as  three- 
dimensional  elastic  beams.  The  beam  elements  are  capable  of  developing  tension, 
compression,  and  bending  stresses.  All  of  the  joints  in  the  structure  arc  considered  rigid. 
The  loads  are  applied  to  the  structure  as  nodal  forces. 

The  wind  force  on  the  structure  is  calculated  by  taking  the  given  wind  velocity  as 
the  basic  wind  speed  measured  at  30  feet  with  the  gust  and  importance  factors  included.  A 
velocity  pressure  profile  is  then  developed  as  a  function  of  height  using  the  appropriate 
exposure  coefficient  in  accordance  with  American  National  Standards  Institute  (ANSI) 
A58.1  -  1982.  The  total  static  wind  force  on  the  structure  is  calculated  as  a  function  of  the 
velocity  pressure,  the  force  coefficients  for  square  towers,  and  the  projected  area  of  the 
structure.  The  wind  is  considered  to  be  approaching  from  the  front,  rear,  and  quartering 
directions.  The  structure  is  designed  to  withstand  the  most  severe  of  these  cases. 

The  wind  force  on  the  backscreen  is  accounted  for  through  the  use  of  three- 
dimensional  shell  elements.  The  shells  arc  stretched  across  the  area  of  the  structure  where 
the  backscreen  is  located.  A  pressure  is  then  applied  to  the  shells  equal  to  the  wind  force 
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on  the  backscreen.  The  density  of  the  shell  elements  is  adjusted  to  account  for  the  weight 
of  the  backscreen  and  the  modulus  of  elasticity  is  set  to  near  zero  in  order  to  avoid  the 
backscreen  contributing  any  stiffness  to  the  structure. 

The  ice  accumulation  is  taken  into  account  by  adjusting  the  density  of  the 
members  to  include  the  ice  and  by  increasing  the  projected  area  of  each  member  by  twice 
the  amount  of  radial  ice  present.  The  force  coefficients  are  recalculated  for  the  members 
because  the  members  are  considered  to  have  rounded  comers  for  iced  conditions.  New 
wind  loads  are  then  calculated  based  on  the  adjusted  values  for  the  projected  areas  and 
force  coefficients  of  the  members. 

The  stresses  calculated  are  compared  with  the  allowable  sties  es  prescribed  in  the 
American  Institute  of  Steel  Construction  (AISC)  Specification  for  the  design  of  structural 
steel  buildings.  The  calculated  stress  for  any  member  is  not  allowed  to  exceed  the 
allowable  stress  stated  in  the  AISC  Code  for  Tension,  Compression,  or  Bending  Stresses. 

A  structural  dynamic  analysis  was  performed  to  determine  the  mechanical 
resonant  frequencies  of  the  structures.  A  wind  gust  profile  is  assumed  and  used  to  excite 
the  structures  and  determine  their  vibrational  characteristics.  A  maximum  vibration 
amplitude  of  0. 1 85  inches  is  allocated  to  the  antenna  from  the  Subclutter  Visibility  Budget 
for  spectral  noise  density.  This  is  in  addition  to  the  initial  displacement  due  to  wind 
velocities  (including  gusts)  of  40  knots  at  30  feet  above  ground  level. 

7.0  SUMMARY 

Six  transmit  and  receive  antenna  pairs  have  been  fabricated  and  installed,  three  for 
the  ECRS  and  three  for  the  WCRS.  The  antennas  have  been  integrated  with  the  rest  of  the 
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radar  equipment  and  are  providing  aircraft  detection  capabilities  for  ranges  of  500  to  2000 
nautical  miles. 
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Herbert  M.  Aumann  and  Frank  G.  Willwerth 
Massachusetts  Institute  of  Technology,  Lincoln  Laboratory 
Lexington,  MA  02173 

ABSTRACT 

The  limitations  of  phased  array  sidelobe  performance  were  explored  with  a 
32-element  linear  array  operating  at  L-band.  With  conventional  far-held  cal¬ 
ibrations,  the  average  sidelobe  level  of  the  array  was  about  -40  dB.  Since  the 
array  contained  T/R  modules  with  12-bit  attenuators  and  12-bit  phase  shifters, 
considerably  lower  sidelobe  performance  had  been  expected. 

Initially,  array  performance  was  thought  to  be  limited  by  inadequate  cali¬ 
brations.  An  examination  of  individual  array  element  patterns  showed  a  mir¬ 
ror  symmetric  ripple  which  could  only  be  attributed  to  edge  effects  in  a  small 
array.  Simulations  indicated  that  more  precise  calibrations  would  not  have 
significantly  lowered  the  sidelobe  level. 

An  adaptive  calibration  technique  was  developed  which  iteratively  adjusted 
the  attenuator  and  phaser  commands  to  create  nulls  in  the  antenna  pattern 
in  the  direction,  of  the  nulls  of  a  theoretical  antenna  patterns.  With  adaptive 
calibrations,  < t>e  average  sidelobe  level  could  be  lowered  to  -60  dB.  Array 
patterns  with  deep  nulls  in  specified  directions  can  also  be  implemented. 
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1.0  INTRODUCTION 


A  phased  array  antenna  with  low  sidelobes  has  advantages  for  both  clutter 
and  interference  suppression.  However,  besides  of  measurement  range  and 
multipath  effects,  a  variety  of  causes  may  make  it  difficult  to  achieve  low 
sidelobe  performance. 

Quantization  effects  are  the  most  often  claimed  to  limit  array  sidelobe 
performance.  By  using  attenuator  and  phase  shifters  with  12-bits  of  control, 
quantization  errors  should  become  insignificant. 

Obviously  it  would  be  desirable  to  build  phased  array  components  with 
sufficient  tight  tolerances  to  make  further  alignment  unnecessary.  In  practice, 
amplitude  and  phase  differences  in  the  T/R  modules  and  transmission  lines 
of  the  beamforming  network  are  almost  unavoidable.  Array  calibrations  can 
be  used  to  compensate  for  these  differences  and  relate  digital  commands  to  an 
actually  applied  excitation  function. 

Usually,  array  element  patterns  are  assumed  to  be  uniform.  Mutual  cou¬ 
pling  and  edge  effects  may  distort  array  element  patterns,  thus  compromising 
the  array  calibrations  and  excitation  function  calculations. 

Rather  than  attempting  to  minimize  the  various  error  sources,  the  phased 
array  calibration  technique  described  in  this  paper  concentrates  on  achieving 
the  best  array  pattern  under  the  given  conditions. 


2.0  ANTENNA  RANGE  AND  PHASED  ARRAY  CHARACTERISTICS 

2.1  Antenna  Test  Range 

Phased  array  calibrations  and  antenna  pattern  measurements  were  made 
with  the  linear  array  in  Figure  1  on  a  100  ft.  ground  reflection  antenna  range. 
The  range  geometry  is  shown  in  Figure  2.  Since  this  range  did  not  meet  the 
far-field  criterion  for  the  test  array,  focusing  of  the  array  was  necessary. 

2.2  Phased  Array  Characteristics 

The  experimental  array  operated  at  L-band  and  was  about  15  ft.  by  2  ft.  in 
size.  It  contained  36  identical  waveguide  radiating  elements.  Two  elements  on 
either  end  of  the  array  were  terminated  to  reduce  edge  effects.  The  remaining 
32  elements  were  connected  to  transmit/receive  modules. 

Half-height,  open-ended  waveguide  radiating  elements  were  chosen  primar¬ 
ily  because  they  were  inexpensive  and  they  had  a  broad,  aspect  angle  insen¬ 
sitive  pattern.  Some  typical  elements  patterns  are  shown  in  Figure  3.  In  the 
center  of  the  array,  the  element  patterns  were  nearly  identical.  However  at  the 
ends  of  the  array,  the  element  patterns  exhibited  a  mirror  symmetric  ripple 
which  was  attributed  to  edge  diffraction.  Ordinary  array  calibrations  cannot 
correct  for  this  ripple. 

Good  sidelobe  performance  in  a  small  array  requires  transmit/receive  mod¬ 
ules  with  very  precise  attenuation  and  phase  control.  Each  of  the  modules  used 
in  the  test  array  had  12-bits  of  amplitude  and  phase  control.  This  was  achieved 
by  using  voltage  controlled  devices  operating  at  an  intermediate  frequency 
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[1],  Typical  characteristics  are  shown  in  Figure  4  and  Figure  5.  A  nomi¬ 
nal  quadratic  approximation  adequately  described  the  attenuator  the  phaser 
characteristics.  There  is  a  weak  attenuator  and  phaser  interaction  which  is 
not  easily  accounted  for  by  a  closed-form  calibration. 

2.3  Array  Performance  with  Standard  Calibrations 

Standard  array  calibrations  were  carried  out  with  the  array  at  broadside. 
That  is,  all  elements,  except  the  element  under  test,  were  terminated,  and 
the  characteristics  of  the  element  under  test  were  determined.  Using  such 
calibrations,  Figure  6  shows  an  antenna  pattern  of  the  test  array  with  a  40 
dB  Chebychev  taper.  Although  the  measured  average  sidelobe  level  of  -39  dB 
was  considered  quite  good,  with  12-bit  quantization  the  average  sidelobe  level 
should  have  been  -45  dB  for  a  40  dB  Chebychev  taper. 

One  way  of  determining  the  limits  of  array  sidelobe  performance  is  to 
measure  the  average  sidelobe  level  as  increasingly  more  severe  array  tapers  are 
commanded.  Figure  7  shows  that  the  array  sidelobe  performance  was  limited 
to  about  -40  dB. 

A  simulation  was  carried  out  using  actually  measured  element  patterns, 
such  as  shown  in  Figure  3,  and  perfect  calibrations.  Of  course,  the  measured 
element  patterns  included  range  multipath  as  well  as  small  array  effects.  The 
resulting  array  sidelobe  performance  was  only  a  few  dB  better  than  the  per¬ 
formance  measured  on  the  antenna  range.  This  indicates  that  the  array  per¬ 
formance  with  standard  calibrations  was  about  as  good  as  could  be  expected. 
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3.0  ADAPTIVE  CALIBRATION  TECHNIQUE 


A  technique  for  calibrating  a  phased  array  from  the  far-held  pattern  has 
been  has  been  described  by  Patton  [2].  Actually,  the  process  consisted  of  first 
calculating  the  far-held  pattern  of  an  array  by  Fourier  transformation  of  planar 
near-held  measurements.  The  array  element  excitations  were  then  computed 
by  inverse  Fourier  transformation  of  the  far-held  pattern  back  to  the  aperture 
plane. 

The  array  calibration  technique  described  here  is  based  on  iteratively  form¬ 
ing  nulls  in  the  antenna  pattern  [3].  The  calibrations  are  considered  ’’adaptive” 
in  that  the  measured  antenna  pattern  is  adjusted  to  match  a  theoretical  model. 
When  the  specihed  nulls  correspond  to  the  nulls  of  an  array  with  uniform  il¬ 
lumination,  the  technique  yields  the  insertion  attenuation  and  phase  of  the 
array  elements.  By  specifying  other  null  directions,  different  array  illumina¬ 
tion  functions  can  be  implemented. 

3.1.  Theory 

The  excitation  function  of  a  phased  array  can  be  calculated 
served  far-held  pattern  by  the  direct  matrix  inversion  method  [4]. 
of  a  linear  A^-element  array  in  the  direction  6  is  given  by 

p(o)  =  Ewnzn(0)  (1) 

n=l 


from  an  ob- 
The  pattern 
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En{On)  ,2tt 

zn{8)  =  — 7T~  expj—  {Hn  mod  A), 
itn  ^ 


where  wn  are  the  complex  excitation  coefficients  and  En(9n)  are  the  element 
patterns.  Usually  it  is  assumed  that  all  array  elements  are  identical  and  the 
range  geometry  is  such  that 


El(0l)  =  E2(d2)  =  ---  =  En(6n). 


(2) 


By  measuring  the  antenna  pattern  in  N  directions,  a  set  of  linear  equations 
can  be  written  and  solved  for  the  N  excitation  coefficients.  These  equations 
can  be  expressed  in  matrix  form  as 


W  =  Z1  P 


(3) 


where 
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p  = 


1 


\  P(Oi)  1 

P(02) 

.  P(On)  . 

In  theory,  measurements  of  the  antenna  pattern  in  any  N  directions  could 
be  chosen  to  determine  the  array  excitation  function.  In  practice,  the  quality 
of  results  obtained  from  a  single  ’’back  transformation”  has  not  been  very  sat¬ 
isfactory.  The  problem  is  that  the  assumption  of  identical  element  patterns 
is  of  limited  validity  and  the  range  geometry  is  not  sufficiently  well  known. 
However,  an  iterative  solution,  with  careful  selection  of  the  measurement  di¬ 
rections,  has  been  successful. 

3.2  Measurement  Selection 

Any  measured  antenna  pattern  will  obviously  be  affected  by  the  specific 
range  geometry.  Thus,  Figure  2  should  be  used  to  carefully  calculate  the  range 
to  each  element,  particularly  to  include  focusing  corrections.  Equation  (1)  also 
indicates  that  differences  in  the  amplitude  and  phase  of  individual  element  pat¬ 
terns  and  errors  in  the  excitation  function  have  similar  effects  on  the  antenna 
pattern.  This  interaction  can  be  minimized  by  making  measurements  in  the 
directions  of  the  pattern  nulls. 

For  commonly  used  array  tapers,  the  excitation  function  can  be  written 
in  factored  form.  Consequently,  the  null  locations  of  the  far-field  pattern  are 
explicitly  known  [5].  In  fact,  the  null  locations  were  probably  used  to  specify 
the  array  sidelobe  level  in  the  first  place.  For  an  Af-element  array,  there  can 
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be  at  most  N  —  1  nulls  in  the  antenna  pattern.  Ideally,  if  the  main  lobe  is  at 
0P  and  the  nulls  are  at  0\,02,  ■  •  •  ,^/v-i  then  the  normalized  pattern  values  are 


r  p(oP)  ' 

'  1  ‘ 

m) 

0 

.  P{0N- .)  . 

.  0 

and  the  matrix  inversion  in  Equation  (3)  will  return  the  theoretical  excitation 
function  WT 


As  a  result  of  the  assumption  made  in  Equation  (2),  the  nulls  in  the  antenna 
pattern  are  entirely  determined  by  the  array  factor.  The  above  calculation  will 
therefore  not  be  affected  by  the  shape  of  the  element  patterns. 


Since  the  array  under  test  and  the  measurement  range  are  not  perfect,  the 
measured  antenna  response  will  not  be  zero  at  the  null  locations,  that  is 


r  p(oP)  1 

"  1  ' 

P(0i) 

— 

.  P(0N-i)  . 

.  C/V-l  . 

However,  Equation  (5)  and  Equation  (3)  can  be  used  to  obtain  an  estimate 
W [■;  of  the  actually  applied  excitation  function.  In  the  following  calibration 
procedure,  attenuation  and  phase  commands  will  be  adjusted  until  differences 
between  the  estimated  and  theoretical  excitation  function  are  a  minimum. 
When  this  is  accomplished,  the  nulls  of  the  measured  and  theoretical  antenna 
pattern  have  been  matched. 


1 
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The  nulls  of  an  antenna  pattern  are  usually  of  no  particular  interest.  How¬ 
ever.  by  forcing  the  antenna  pattern  to  have  nulls  in  the  proper  directions,  a 
well  formed  antenna  pattern  is  virtually  assured. 


It  should  be  noted  that  the  process  is  mathematically  identical  to  adaptive 
nulling  by  sample  matrix  inversion  [6].  That  is  the  array  is  required  to  receive  a 
desired  signal  from  a  given  direction,  and  it  is  required  to  suppresses  undesired 
signals  in  N  —  1  carefully  specified  other  directions. 

3.3  Calibration  Procedure 


Usually,  the  first  step  in  the  array  calibrations  is  to  correct  for  array  element 
insertion  attenuation  and  phase  differences.  Initially,  an  excitation  function 
Wc  is  commanded  without  corrections  to  yield  a  uniform  illumination  Wy 


Wc  =  WT  = 


'  1  ' 
1 


(6) 


Li 


For  example,  the  resulting  antenna  pattern  of  the  uncalibrated  test  array  is 
shown  in  Figure  8.  The  amplitude  and  phase  of  the  antenna  pattern  is  mea¬ 
sured  at  broadside  and  at  the  null  locations  corresponding  to  a  uniform  illu¬ 
mination.  From  Equation  (3)  an  estimate  of  the  actually  applied  excitation 
function  W e(i)  is  obtained.  This  estimate  is  then  used  to  update  the  com¬ 
manded  excitation  function  on  the  next  iteration  as  follows 


139 


Wc(.  +  i)  =  wc(;)[J^-|  (7) 

vv  e{i)  J 

where  k  <  1  is  a  gain  constant.  The  actual  commands  have  to  be  expressed 
in  counts  based  on  nominal  attenuator  and  phaser  characteristics.  The  above 
process  is  repeated  until 

WE{i)  *  WT,  (8) 

and  there  is  no  further  improvement  in  the  null  depths.  For  example,  Figure  9 
illustrates  that  after  4  iterations  a  nearly  perfect  pattern  expected  for  an  array 
with  uniform  illumination  has  been  achieved. 

Under  the  above  conditions,  the  array  element  insertion  attenuation  and 
phase  are  given  by  the  reciprocal  of  Wq-  Shown  in  Figure  10  and  Figure  11 
are  the  insertion  attenuation  and  phase  as  determined  by  standard  calibrations 
and  as  measured  by  adaptive  calibrations.  Standard  array  calibrations  are 
usually  carried  out  with  a  source  at  broadside.  It  was  shown  in  Figure  3 
that  the  element  patterns  are  not  entirely  identical  at  all  angles.  Adaptive 
calibrations  include  element  pattern  differences  and  therefore  yield  slightly 
different  calibration  results. 

As  with  any  feedback  algorithm,  there  is  concern  with  convergence.  For 
example,  an  instability  may  occur  if  the  Z-matrix  does  not  adequately  model 
the  array  and  range  geometry.  Another  cause  for  instability  may  be  a  lack  of 
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monotonicity  in  the  attenuator  and  phaser  characteristics. 

3.4  Array  Performance  with  Adaptive  Calibrations 

3.4.1  Sidelobe  Suppression  by  Chebychev  Taper 

The  calibration  technique  is  not  limited  to  an  array  with  uniform  illumina¬ 
tion.  For  example,  shown  in  Figure  12  is  the  adaptively  determined  excitation 
function  required  to  implement  a  40  dB  Chebychev  taper.  A  parabolic  phase 
distortion  necessary  to  phase  focusing  is  clearly  evident.  Although  the  dif¬ 
ferences  between  standard  and  adaptively  determined  excitation  function  are 
usually  quite  small,  the  resultant  array  pattern,  in  Figure  13,  represents  a 
considerable  improvement  over  the  corresponding  pattern  with  standard  cali¬ 
brations  shown  in  Figure  6. 

With  adaptive  calibrations  exceedingly  low  sidelobe  performance  could  be 
demonstrated.  Figure  14  shows  an  antenna  pattern  with  a  60  dB  Chebychev 
taper.  At  that  level,  the  40  dB  dynamic  range  of  the  attenuators  becomes  the 
limiting  factor. 

Shown  in  Figure  15  is  the  sidelobe  performance  with  adaptive  calibrations. 
With  adaptive  calibrations,  the  array  performance  comes  much  closes  to  the 
theoretical  array  performance  expected  from  an  array  with  12- bits  of  amplitude 
and  phase  control. 

3.4.2  Sidelobe  Suppression  by  Null  Placement 

In  the  above  section  the  antenna  pattern  null  directions  were  chosen  to 
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correspond  to  standard  Chebychev  tapers.  This  resulted  in  very  predictable 
array  performance.  However,  uniformly  low  sidelobes  may  not  be  required  in 
all  applications.  Rather  than  accepting  the  loss  in  efficiency  and  broadening  of 
the  mainlobe,  selective  null  placement  in  the  direction  of  a  known  interference 
source  may  be  the  preferred  approach. 

For  example,  shown  in  Figure  16  is  an  adaptively  calibrated  array  pattern 
with  a  30  dB  Chebychev  taper  and  an  additional  null  at  -30  degrees.  Again, 
sidelobes  in  the  direction  of  the  interference  have  been  suppressed  by  60  dB. 

4.0  SUMMARY 

An  adaptive  calibration  technique  for  a  phased  array  was  described.  Since 
the  algorithm  derived  array  excitation  function  directly  from  the  measured 
far-field  pattern,  no  build-in  test  capability  is  required. 

The  technique  is  applicable  not  only  to  implement  standard  array  tapers, 
but  can  also  be  used  to  implement  antenna  patterns  with  deep  nulls  in  specified 
directions. 
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Figure  1:  Experimental  Phased  Array 


Figure  2:  Antenna  Test  Range 
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Figure  3:  Typical  Array  Element  Patterns 
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Figure  4:  Typical  Attenuator-Characteristics 
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.Figure  6:  Antenna  Pattern  with  Standard  Calibrations 


Figure  7:  Array  Sidelobe  Performance  with  Standard  Calibrations 
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Figure  8:  Uncalibrated  Array  Pattern 


Figure  10:  Comparison  of  Insertion  Attenuation 


Figure  11:  Comparison  of  Insertion  Phase 
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o  THEORY  40  dB  CHEBYCHEV 

□  ADAPTIVE  CALIBRATIONS  FREQUENCY  =  1.3  GHz 


Figure  12:  Excitation  Function  for  a  40  dB  Chebychev  Taper 
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Figure  12:  Array  Pattern  with  10  dB  Chebychev  Taper 
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Figure  14:  Array  Pattern  with  60  dB  Chebychev  Taper 


Figure  15:  Array  Sidelobe  Performance  with  Adaptive  Calibrations 
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Figure  16:  Adaptively  Calibrated  Antenna  Pattern  with  Null 
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Applications  of  Self-Steered  Phased  Arrays 


Dean  A.  Paschen 

Ball  Communication  Systems  Division 
Broomfield,  CO  80038 


ABSTRACT 


Self-steered  phased  arrays  were  studied  in  the  early 
1960's  to  aid  in  communications  with  non-stationary  low 
earth  orbit  satellites.  Introduction  and  widespread 
use  of  geosynchronous  satellites  eliminated  most  of  the 
necessity  for  tracking,  and  active  transponders  reduced 
the  need  for  large  earth  stations.  Consequently, 
interest  in  self-steered  arrays  declined.  In  spite  of 
this  decline,  several  uses  for  self-steered  array 
techniques  still  remain.  The  beneficial  applications 
range  from  lightweight,  portable,  quick  erect  antenna 
systems  to  fixed  installation,  ultra  high  gain  ground 
stations.  Antenna  systems  which  operate  in  a  disturbed 
atmosphere  also  benefit  from  the  automatic  phasing 
inherent  in  a  self-steered  array.  The  purpose  of  this 
paper  is  to  discuss  some  of  these  applications  along 
with  methods  of  implementing  the  self-steered  array. 
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1.0  INTRODUCTION  AND  HISTORY 

Self-steered  arrays  were  first  investigated  in  the 
early  1960's  for  use  with  low-earth  orbit  satellites. 
Two  types  of  self-steered  arrays  were  studied  for  two 
different  problems.  Retrodirective  or  Van  Atta  arrays 
were  studied  for  use  in  satellite  communications  to 
minimize  the  problems  of  achieving  high  transmit  and 
receive  gain  at  the  satellite  without  tracking  the 
received  signal.  The  other  use  of  self-steered  array 
technology  was  the  approach  presented  by  Eberle1.  Four 
30  foot  diameter  reflector  antennas  were  combined 
coherently  to  achieve  the  same  performance  as  a  60  foot 
diameter  single  parabolic  reflector.  The  array  was 
designed  for  use  with  Echo  I  and  Echo  II  passive 
satellites,  and  this  configuration  is  shown  in  Figure 
1-1.  Combination  of  the  signals  provided  maximum  array 
gain  in  the  direction  of  the  signal  source.  Advantages 
of  the  approach  were  that  combination  of  several 
smaller  antennas  was  more  cost  effective  than  a  single 
larger  antenna,  that  the  acquisition  of  the  satellite 
was  easier  due  to  the  larger  beamwidth  of  the  smaller 
antennas,  and  that  the  system  compensated  for  phase 
variations  due  to  atmospheric  inhomogeneities. 
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The  introduction  of  geosynchronous  satellites  using 
earth  coverage  antennas  eliminated  the  need  for 
tracking  on  the  satellite,  and  the  need  for 
retrodirective  arrays  disappeared.  The  use  of  active 
transponders  in  these  same  satellites  began  to  reduce 
the  requirement  for  very  large  reflector  antennas,  and 
the  study  of  self-steered  arrays  declined. 


2.0  SELF-STEERED  ARRAYS  PROVIDE  BENEFITS 

The  benefits  which  motivated  the  study  of  self-steered 
arrays  in  the  past  still  exist  for  large  antenna 
systems.  Those  benefits  are  lower  cost,  simplified 
acquisition  and  tracking  due  to  a  wider  field-of-view 
(FOV) ,  and  compensation  for  atmospheric  variations  from 
element  to  element.  Some  additional  uses  which  have 
been  investigated  more  recently  are  highly  portable 
satellite  ground  terminals,  large  baseline  arrays  where 
the  source  may  not  be  in  the  far  field,  and  arrays 
which  are  resistant  to  multipath  fading. 

The  cost  benefit  of  self-steered  arrays  is  based  on  the 

property  that  very  large  reflector  antennas  have  a  cost 

.  2 

versus  area  which  is  exponential  .  A  study  by  Potter 
et  al.  for  deep  space  dishes  found  that  cost  was 
roughly  proportional  to  area  to  the  1.4  power  for 

3 

constant  surface  tolerance  .  Self-steering  concepts 
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have  a  higher  initial  cost  due  to  the  electronics,  but 
the  increase  with  total  aperture  is  linear.  Figure  2-1 
shows  this  relationship  with  the  assumption  that  the 
electronics  cost  is  equal  to  the  cost  of  a  single  unit 
of  aperture  area.  The  cost  benefit  for  large  total 
aperture  areas  is  obvious. 

A  self-steered  array  is  made  up  of  elements  which  are 

combined  in  the  manner  of  a  phased  array.  The  elements 

may  be  any  antenna  type  although  parabolic  reflectors, 

planar  fixed  beam  arrays,  and  phased  arrays  have  been 

the  primary  types  considered.  One  advantage  of 

self-steering  is  that  since  each  element  or  subarray 

has  a  wider  beamwidth  than  the  composite  array,  the 

required  pointing  accuracy  for  each  element  is  reduced. 

Also,  the  FOV  of  the  system  is  actually  the  FOV  of  the 

* 

individual  elements,  so  that  acquisition  of  the  signal 
is  simplified. 

The  benefit  of  compensation  for  atmospheric 
inhomogeneities  which  was  mentioned  in  the  earliest 
work  is  important  for  systems  designed  to  operate  in 
nuclear  battlefields,  as  the  atmospheric  scintillation 
is  a  major  concern  for  communications  links.  In  this 
type  of  operation,  phase  variations  may  occur  which 
cause  the  planar  wavefront  to  be  distorted  such  that 
standard  antenna  techniques  suffer  signal  reductions. 
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Based  on  the  statistics  of  the  environment,  a  single 
aperture  may  be  broken  into  pieces  so  that  the  phase 
variation  across  each  piece  is  relatively  small.  The 
self -steering  technique  can  then  be  used  to  compensate 
for  the  larger  variations  from  aperture  to  aperture. 

One  of  the  major  incentives  for  the  investigation  of 
self-steered  phased  arrays  at  Ball  Communications 
Systems  Division  (BCSD)  was  the  potential  for  use  in 
highly  portable  satellite  ground  station  antennas. 
Figure  2-2  shows  a  conceptual  drawing  of  a  self-steered 
array  for  use  with  the  Defense  Satellite  Communications 
System  (DSCS) .  The  goal  is  to  provide  the  equivalent 
performance  of  a  20  foot  diameter  reflector  antenna  and 
allow  for  storage  in  the  aisle  of  a  S-250  or  larger 
shelter.  Without  self-steering,  the  panels  which  make 
up  the  antenna  system  would  have  to  be  positioned  very 
accurately,  and  the  cables,  filters,  amplifiers,  and 
other  components  would  require  phase  matched 
construction.  All  of  these  requirements  are  eliminated 
in  a  self-steered  system,  and  as  a  result,  the 
deployment  time  and  cost  of  this  type  of  antenna  array 
are  reduced  significantly. 

If  a  very  large  baseline  array  were  constructed,  the 

possibility  would  exist  for  a  signal  of  interest  to  be 

in  the  near  field  region  of  the  array.  Under  this 
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scenario,  a  self-steered  array  would  continue  to 
provide  maximum  gain  by  effectively  compensating  for 
the  incident  spherical  phase  front. 

A  common  problem  for  antennas  operating  at  low 
elevation  angles  or  in  urban  environments  is  multipath 
fading.  A  self-steered  antenna  array  will  provide  the 
best  phase  at  each  element  for  maximizing  the  total 
signal  to  the  receiver.  Brennan  showed  that  the 
self-phased  array  is  resistant  to  multipath  fading,  and 
that  under  certain  conditions  may  achieve  total  gains 
which  are  higher  than  that  obtained  for  the  same  array 

4 

in  an  isolated  environment  .  A  good  example  would  be 

an  array  operated  in  the  standard  manner  versus  a 

self-phased  array  operated  in  front  of  a  reflecting 

surface.  Figure  2-3  shows  how  the  capture  of  the 

energy  from  the  surface  may  provide  higher  gain  than 

for  the  array  operated  independently.  One  application 

which  takes  advantage  of  the  multipath  rejection 

benefit,  especially  in  urban  areas,  is  the  one 

developed  by  Blaupunkt  for  use  in  high  performance  FM 

5 

systems  for  automobiles  .  The  system  uses  four  foil 
antennas  in  the  corners  of  the  car  body  and  combines 
the  antennas  for  maximum  signal  using  self-steered 
array  techniques. 
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3.0  SELF-STEERED  ARRAY  TECHNIQUES 

The  motivation  for  self-steering  of  phased  arrays 
becomes  apparent  when  some  of  the  practical  aspects  are 
considered.  A  conventional  phased  array  is  shown  in 
Figure  3-1.  Here  the  phases  are  calculated  under  the 
assumption  that  the  position  of  the  elements  are  known, 
that  the  component  insertion  phases  are  also  known,  and 
that  the  incoming  signal  has  a  wavefront  which  is 
planar.  Figure  3-2  shows  a  view  of  an  array  where  all 
of  these  characteristics  are  violated.  In  this  case, 
severe  reductions  in  the  signal  level  to  the  receiver 
may  result.  One  type  of  processing  technique  which  is 
analogous,  is  the  code  acquisition  problem  in 
spread-spectrum  systems.  Figure  3-3  shows  how  the 
autocorrelation  of  the  spreading  sequence  with  a 
replica  changes  with  the  time  offset.  Maximum 
autocorrelation  indicates  time  alignment.  The  same 
process  can  be  used  in  an  antenna  array  by  using 
signals  from  separate  elements  as  the  comparison 
signals.  The  signal  at  each  element  is  identical  in 
content  with  only  displacement  in  phase  (time) .  The 
autocorrelation  of  the  two  signals  as  shown  in  Figure 
3-4  provides  the  necessary  information  to  compensate 
for  the  phase  difference  using  phase  shifters. 
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All  methods  of  self-steering  must  perform  two  basic 

functions.  First,  the  relative  phase  between  the 

elements  must  be  measured,  and  second,  the  relative 

phase  difference  must  be  eliminated.  Figure  3-5  shows 

one  of  the  earliest  methods  of  performing  the  phase 

alignment  necessary  for  self-steering.  Each  element  of 

the  array  performs  a  downconversion  to  intermediate 

frequency  (IF)  such  that  the  IF  is  phase-locked  to  an 

internal  reference  frequency.  The  phase-locked  loop 

performs  both  the  measurement  and  phase  difference 

elimination.  Since  the  output  phases  of  the  IF  outputs 

will  all  have  phase  relationship  to  the  reference 

oscillator,  the  relative  phase  between  the  IF  outputs 

will  be  zero.  This  approach  also  has  the  ability  to 

compensate  for  differential  frequency  shifts  between 

elements  due  to  doppler  shift  variations.  One 

potential  disadvantage  of  this  approach  is  that  each 

element  becomes  a  duplicate  receiver  front  end  with  the 

combined  output  at  the  IF  frequency.  An  alternative 

approach  which  provides  a  combined  output  at  the  radio 

frequency  (RF)  is  shown  in  Figure  3-6.  In  this 

example,  signal  from  each  element  is  coupled  to  both 

sides  of  a  phase  detector  or  autocorrelator.  One 

element  is  chosen  as  the  reference,  and  each  of  the 

other  elements  is  compared  in  sequence  to  determine  the 

change  in  the  phase  shifter  setting  to  provide  coherent 

signal  combination.  The  approach  shown  allows  any 
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element  to  be  the  reference  to  minimize  single  point 
failures.  Additionally,  a  single  element  can  be  fed 
simultaneously  to  both  sides  of  the  phase  detector  for 
calibration  purposes. 

A  self-steered  array  built  from  aperture  antenna 

elements  will  have  radiation  patterns  which  are 

dominated  by  the  characteristics  of  the  individual 

apertures.  The  sidelobes  of  the  composite  array  will 

fall  within  the  sidelobes  of  the  aperture  pattern,  and 

grating  lobe  formation  which  may  occur  will  fall  within 

the  main  beam.  The  total  antenna  gain  is  independent 

of  the  composite  pattern  and  depends  only  on  the  total 

aperture  area  which  is  combined.  To  illustrate  this 

principle.  Figure  3-7  shows  an  eight  element  subarray 

with  0.75  A  spacing.  Figure  3-8  shows  the  array  factor 

for  eight  of  these  subarrays  with  periodic  spacing  of 

10  A.  Figure  3-9  shows  the  same  array  factor  as  Figure 

3-8  with  a  normally  distributed  random  error  in  the 

position  having  mean  of  zero  and  standard  deviation  of 

two  wavelengths.  Practical  systems  will  be  of  the 

random  type.  Figure  3-10  and  3-11  show  the  product  of 

the  array  factor  and  the  subarray  element  factor  with 

the  element  factor  plotted  for  comparison.  In  both 

cases,  the  maximum  composite  pattern  falls  within  the 

envelope  of  the  subarray  pattern.  However,  the 

randomly  positioned  array  will  always  tend  to  provide 
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lower  peak  sidelobe  and  grating  lobe  levels  due  to  the 
reduced  gain  of  auxiliary  peaks  in  the  array  factor  for 
the  subarray  panels.  Total  energy  in  the  sidelobe 
region  is  the  same  in  either  case;  the  random 
positioning  of  the  panels  simply  reduces  the  peak 
levels . 

CONCLUSION 

Self-steered  phased  arrays  have  several  benefits  which 
should  find  renewed  interest  in  certain  applications. 
The  potential  for  cost  reduction  in  large  antennas  is 
certainly  strong  motivation  in  an  era  of  ever 
tightening  budgets.  Simpler  acquisition  and  tracking 
also  offer  cost  reduction  benefits.  Interest  in 
operation  of  large  antennas  in  a  disturbed  atmosphere 
or  in  the  near  field  encourage  the  study  of 
self-phasing  concepts.  Portable  operation  is  a  major 
new  area  for  the  development  and  use  of  self-steered 
phased  array  techniques.  Additionally,  resistance  to 
multipath  fading  encourages  self-steered  array  use  in 
low  elevation  or  urban  environments.  In  fact,  the 
combination  of  auto-tracking  with  multipath  rejection 
has  already  encouraged  use  in  commercial  mobile 
communications.  The  self -steered  phased  array  provides 
increased  flexibility  in  solving  complex  user 

requirements  for  many  types  of  applications. 
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Figure  1-1.  Four  thirty-foot  reflectors  combine  for 
Echo  satellite  communications. 
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2-1.  Reflector  cost  increases  faster  than  the 
cost  of  self-steered  phased  arrays. 
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Figure  2-2.  Self-steered  array  for  DSCS  provides 
increased  mobility  with  lower  volume. 


Large  structure 


Figure  2-3.  Additional  energy  may  be  captured  using 

self-steering  concepts. 
t6: 


Figure  3-2.  Errors  cause  gain  reduction  in 
conventional  phased  arrays. 
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Offset  (a) 


Figure  3-3.  Autocorrelation  allows  determination  of 
time  offsets  in  spread-spectrum  systems 


Figure  3-4.  comparison  of  element  phases  provides 
information  to  set  phase  shifters. 
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Figure  3-6.  Alternate  method  sets  phases  at  RF. 
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Subarray  pattern  provides  an  envelope  for 

the  composite  array  radiation  pattern. 
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Figure  3-7 


SEPTEMBER  1990 


ARRAY  THINNING 
USING  THE 

IMAGE  ELEMENT  ANTENNA 


PRESENTED  AT  THE  1990  SYMPOSIUM  ON  ANTENNA  APPLICATIONS 
ROBERT  ALLERTON  PARK,  ILLINOIS 


by 

JOE  KOBUS 
ROBERT  SHILLINGBURG 
and 

RON  KIELMEYER 


MOTOROLA  INC. 

TACTICAL  ELECTRONICS  DIVISION 
GOVERNMENT  ELECTRONICS  GROUP 
SCOTTSDALE,  ARIZONA 


172 


ABSTRACT 


Current  image-antenna  array  technology  is  presented.  The  dielectric-slab 
image  radiator,  discovered  by  Von  Trentini,  was  first  incorporated  into  thinned 
slot  arrays  at  Motorola’s  Government  Electronics  Group  (ll.S.  Patent  No. 
3990078).  This  paper  presents  some  previously  unconsidered  factors  affecting 
image-array  performance.  An  excessively  high  image-slab  permittivity  can, 
owing  to  enhanced  inter-slot  coupling,  compromise  the  efficiency  and  bandwidth 
of  even  a  highly  thinned  array.  Optimum  performance  (maximum  gain,  minimum 
sidelobes,  and  maximum  bandwidth)  in  general  demands  active  impedance 
matching  with  Clavin-pin  compensated  slots.  ThU  is  illustrated  for  a  specific 
hardware  application.  A  simple  analytical  extension  of  the  pioneering  work  of 
Von  Trentini  [1]  and  Immell  [2]  is  outlined  for  analysis.  In  addition,  a  fi«ld- 
theoretical  solution  is  sketched  which,  a’^tough  not  yet  fully  exploited,  should 
ultimately  yield  a  sufficiently  comprehensive  analytical  model  for  image  arrays. 
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1.  INTRODUCTION 

Von  Trentini  [1]  showed  that  the  gain  of  a  single  radiating  element  could 
be  increased  to  arbitrarily  large  values  simply  by  placing  the  element  between 
two  reflecting  planes,  one  of  which  is  a  semi-refiector  and  the  other  a  perfect 
reflector.  This  configuration,  using  a  slot  as  the  radiating  element  and  a 
dielectric  slab  as  the  semi-reflector,  is  displayed  in  Figure  1.  Bellee  et  al  [31 
called  it  an  Image  Antenna  since  the  gain  enhancement  can  be  explained  in 
terms  of  image  theory.  Immell  and  Sasser  [2]  introduced  directional  impedances 
into  Von  Trentini’s  model  in  order  to  calculate  the  reflection  coefficient  of  the 
partially  reflecting  dielectric  sheet.  Jackson  [41  recently  used  a  transmission  line 
analogy  to  compute  the  gain  enhancement  of  printed  circuit  antennas  radiating 
through  partially  reflecting  sheets. 

Subsequent  to  the  pioneering  efforts  [11  -  [3],  modifications  to  optical  and 
transmission  line  analyses  have  been  used  as  aids  to  image  array  design  at 
Motorola  (see  for  example  [11]).  While  providing  physical  insights  and  first- 
order  design  relationships,  they  do  not,  however,  yield  adequate  models  for 
inter-element  coupling.  Mutual  coupling  becomes  important  in  this  type  of  gain 
enhancement  scheme  for  thinned  arrays  when  the  permittivity  of  the  partially 
reflecting  dielectric  screen  becomes  significantly  large. 

Section  2  gives  an  outline  of  the  truncated  image  element  model,  useful 
as  a  simplified  approach  to  array  design.  Section  3  sketches  a  field  theoretical 
solution  which  will  ultimately  yield  an  adequate  mutual  coupling  model. 
Performance  of  a  specific  hardware  design  is  discussed  In  Section  4. 

2.  TRUNCATED  IMAGE-ELEMENT  MODEL 

Figure  2  presents  the  geometry  basic  to  understanding  broadside  array 
design.  Fundamental  theory  dictates  that  no  more  gain  can  be  realized  from  a 
finite  aperture  tftan 

G*=4  tta/x1  •  <’> 

where  A  is  Hie  total  area.  If  A  is  divided  Into  N  equal  S  x  S  sub-apertures,  then 
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clearly  for  (1)  to  hold  we  must  have 

Gj*S —  ^  A  /  N  >  (2) 

for  each  sub-aperture.  Attempts  to  design  an  array  for  which  element  gains  are 
forced  to  be  greater  than  G3  will  be  frustrated  by  mutual  coupling.  On  the  other 
hand,  each  gain  must  be  nominally  no  less  than  that  given  by  (2),  else  the 
composite  array  gain  will  fall  short  of  (1),  i.e.,  excessive  power  will  be  radiated 
into  grating  lobes. 


Refer  to  Figure  1.  The  analysis  of  Von  Trentini  applies  to  dielectric  slabs 
of  infinite  extent.  Clearly  this  is  not  appropriate  for  arrays.  An  approach  is  to 
assume  mutual  coupling  effectively  truncates  each  sub-aperture  in  Figure  2  to 
an  area  S  x  S.  Thus  we  need  to  establish  a  relationship  for  the  gain 
enhancement  relative  to  its  feed  slot  of  an  image  antenna  which  is 
constrained  to  a  finite  aperture.  This  is  done  by  limiting  the  number  M  of  the 
equivalent  "endfire  images"  in  Von  Trentini’s  model  (Figure  1).  Analysis  yields 

^ = o  -  k  f,  k)  ■  (3> 

where  M  (an  integer)  is  the  total  number  of  sources, 

c  fo  £*\  —  LNf)  +SIN*(%€)  ,  (4, 


F(?,a  = 


Sinh*(-j  '-Nf)  -f-  SlN^-jS) 

^  j cos  e  -tr  +  y 

X  =  free  space  wavelength 
0  =  angle  off  normal  to  dielectric  slab 

Sp  =  spacing  of  the  slab  from  the 
groundplane 
=  magnitude  of  the  slab 
reflection  coefficient 
^  =  phase  of  the  slab 
reflection  coefficient. 
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If  F0  is  the  band  center  frequency,  then  it  can  be  shown  [2]  that  maximum 
gain  occurs  at  d  =  0  (Figure  1)  when 


Sp  ~  Xo / x 


(6) 


and  the  dielectric  sheet  thickness  is 


r  = 


,  (7) 

where  €  is  the  dielectric  constant.  For  any  frequency  F  near  F0  it  then  follows 
from  an  analysis  of  directional  impedances  that 


?  =  f(Fo)  = 


6-1 


€+  I 


and 


T'l-tr 


)  *  m-  ~  £)] 


(8) 


(9) 


Clavin-compensated  slots  [5J,  (6],  [7],  [8],  [9],  are  used  to  minimize 
radiation  along  the  groundpiane.  They  have  been  found  to  reduce  inter-element 
coupling  as  well  as  Brewster-angle  and  grating-lobe  radiation.  For  Clavin 
compensation  it  is  easily  shown  that  a  slot  has  7.8  dBI  gain.  Thus  the  total 
gain  of  an  image  radiator  becomes 

<50=  l°LO&lo(Si)  +  7.8  +ao  LoG|0Ccose) 

9 

where  G*  is  given  by  (3). 


(10) 


Refer  to  Figure  3.  Calculation  of  field  radiation  patterns  for  an  array  of 
N  (here  assumed  Identical)  image  antenna  elements  is  made  from 


N  t 


e  =  >  E.e* 


(11) 


where 


■*- — I 
Y)z 


Ee-Ee,Cn)f^'  $<©) 


(12) 


Here  f(e)  is  the  feed  slot  pattern  sans  dielectric,  and  GJs  given  by  (3);  E0(n)  are 
the  amplitude  weightings  in  the  corporate  (combiner)  network.  In  terms  of  the 
Figure  3  coordinate  system  we  also  have 

s l  =[MRM]s,Ne  cosC 

= c°s /  +  Y(t>)s|n<^]s|ns  •  <13) 

Theoretical  image  antenna  E-plane  patterns  are  displayed  in  Figure  4. 
These  patterns  are  computed  at  a  frequency  for  which  the  image  slab  is  X / 4 
thick  and  spaced  ?V2  from  the  ground  plane.  The  slot  illumination  pattern  is 
assumed  to  obey 

which  requires  Clavin  pins.  Observe  that  the  untruncated  patterns  shown  in 
Figure  4  display  considerably  more  variation  with  €  than  do  the  patterns 
truncated  at  M  =  15,  a  value  corresponding  to  typical  image  element-to-element 
spacings  in  practical  broadside  arrays.  This  variation  is  easily  understood  by 
referring  again  to  Figure  t.  The  effective  amplitudes  of  successive  virtual 
Images  are  given  by 

V(< v*)  —  ECT  j0™*  (15) 

Thus,  taking  (8)  into  account,  the  V(m)  are  more  rapidly  attenuated  for  lower 
£  slabs.  Consequently,  the  effect  of  removing  images  is  more  dramatic  for 
higher  dielectric  constants  than  for  lower.  The  physics  of  the  situation  is  that 
less  power  is  radiated  per  "bounce"  from  higher  £  sheets  than  from  lower.  It 
is  also  why,  given  a  fixed  average  spacing,  inter-element  coupling  becomes  a 
greater  problem  for  higher  £  slabs;  more  power  remains  to  be  coupled  to 
neighbors  in  broadside  arrays. 

Comparison  of  E-plane  patterns  computed  for  a  typical  24-element  array 
are  presented  in  Figure  5.  Observe  that  patterns  obtained  with  untruncated 
(incorrect)  elements  exhibit  lower  first  (3dB)  and  second  side  lobes  than  do 
patterns  for  truncated  elements.  Recorded  patterns  of  arrays  larger  than  two 
image  elements  per  quadrant  commonly  display  higher  skMobes  than  those 
calculated  using  untruncated  elements.  Given  Figure  4,  these  higher  sidelobes 


f(e)=  coS© 


result  from  the  broader  truncated  image  element  patterns. 


In  conjunction  with  the  truncation  model,  Motorola  has  pursued  an 
intensive  laboratory  effort  to  understand  array  thinning.  The  paramount 
techniques  applied  in  this  work  were  measurement  and  analysis  of  active 
element  patterns,  gain,  and  impedance.  An  isolated  corporate  feed  network, 
Figure  6,  with  adjustable  amplitude  and  phase  weights  was  used  as  a 
measurement  tool  for  verification  of  low-sidelobe  image  array  designs  of  up  to 
28  elements.  Standard  array  theory  was  used  to  calculate  sum  and  difference 
patterns;  first-order  impact  of  coupling  was  approximated  by  truncation  of 
individual  element  (sub)  apertures.  This  model  yields  array  gain  in  excellent 
agreement  with  measurement  when  the  average  element  gain  (2)  is  assumed 
to  be  equal  to  (10),  but  only  for  sufficiently  low  €  image  slabs.  The  truncation 
model  also  yields  patterns  in  reasonable  agreement  with  those  recorded  for 
arrays,  but  only  for  a  number  (M&>15)  of  image  elements  different  from  that 
(M  t"?  8)  required  for  gain  computations.  Clearly  this  simplified  approach  is 
useful,  but  also  incorrect,  especially  since,  in  any  case,  all  elements  in  a  small 
array  hardly  have  identical  patterns  or  gain.  Consequently  a  field-theoretical 
solution  to  the  problem  of  thinned  arrays  whose  elements  are  modified  with  a 
dielectric  sheet  is  desired.  An  effort  in  this  direction  is  outlined  in  Section  3. 

As  an  example,  however,  of  the  research  findings  in  the  context  of  the 
truncation  model,  consider  the  24-element  array  of  Figure  7.  The  gain 
calculation  is  superimposed  on  the  Figure  8  recorded  data  for  image  element 
No.  1.  Agreement  is  good  if  one  recognizes  that  these  data  include  impedance 
mismatch  loss  and  stripline  feed  loss  as  functions  of  frequency  whereas  the 
computations  do  not.  On  the  other  hand,  measurements  for  a  dielectric  sheet 
with  relative  permittivity  of  37  fall  considerably  short  of  similar  calculations. 
Evidently  £  =37  is  too  high  whereas  £  =25  is  probably  close  to  optimum  for  the 
average  element  spacing  in  this  particular  array.  Further  data  support  the 
conclusion  that  a  relative  permittivity  of  37  results  in  tighter  slot-to-slot  coupling 
(Figure  9)  than  does  £  =  25.  The  truncation  model  clearly  does  not  adequately 
account  for  mutual  coupling. 
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Element  gains  for  the  two  permittivities  are  presented  in  Figures  10-11. 
Embedded  gain  was  recorded  at  an  SMA  coaxial  connector  near  each  slot 
radiator.  Series  stripline  [10]  is  used  to  couple  to  each  slot  in  the  groundplane 
opposite  an  image  slab  supported  by  a  hexagonal-cell  honeycomb  layer. 

As  expected,  inner  radiators  display  less  gain  than  outer  ones.  Since  they 
are  virtually  surrounded,  inner  elements  are  much  more  affected  by  coupling 
than  are  outer  elements,  exposing  another  shortcoming  of  the  truncation  model. 
Furthermore,  these  data  constitute  additional  evidence  in  support  of  the 
conclusion  that  lowered  permittivities  yield  less  coupling. 

Table  1  reveals  the  impact  of  image-slab  permittivity  on  aperture  efficiency. 
These  values  are  computed  with  the  aid  of  the  basic  array  relationship 

0-=  0cm 

n 

valid  for  generally  dissimilar  element  gains  G0(n).  Observe  that  when  all 
radiators  are  equally  weighted,  (16)  yields 


which  if  all  are  identical  reduces  to  (2).  Recorded  power  (dBI)  gains,  given  in 
Figures  (10)  and  (11),  respectively,  were  converted  to  voltage  ratio  gams 
^G»(n)  and  substituted  into  (17).  Table  1  gains  determined  in  this  manner  were 
also  compared  with  corporate-fed  gains  recorded  using  the  universal  test 
arrangement  of  Figure  6.  After  correction  for  feed  network  insertion  loss,  the 
two  independently  determined  array  gains  are  virtually  identical  (typically  within 
±  .3dB). 

Figure  12  shows  E-plane  patterns  recorded  for  a  24-element  array  in  the 
Figure  6  universal  feed  network.  Comparison  with  theoretical  patterns  obtained 
using  (10)-(13)  together  with  (3)-(5)  is  excellent. 
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3.  FIELD  THEORETICAL  ANALYSIS 

This  section  presents  a  field  theoretical  analysis  of  a  single 
Clavin  pin  compensated  slot  in  an  infinite  ground  plane  with  M 
dielectric  plates  stacked  in  the  z  direction.  The  image  antenna  is  a 
special  case  for  which  M=3.  Normally,  erM  would  be  air  but  the 

analysis  which  follows  is  kept  general  at  this  point.  Figure  13 
shows  the  specific  case  of  M=3. 

The  single  Clavin  pin  compensated  slot  solution  is  obtained  by 
solving  two  seperate  problems.  The  first  problem  is  that  of  a  single 
narrow  slot  enhanced  with  the  dielectric  slabs.  For  this  problem  a 
solution  to  the  fields  is  found  in  terms  of  the  assumed  Electric  field 
distribution  in  the  slot.  The  second  problem  is  that  of  a  monopole 
antenna  protruding  into  eri  in  the  z  direction.  A  solution  for  the 

fields  is  found  in  terms  of  an  assumed  current  distribution  on  the 
monopoles.  A  moment  method  solution  is  then  applied  to  the  problem 
to  determine  the  current  distribution  along  the  monopoles  in  terms 
of  the  assumed  Electric  field  distribution  of  the  slot. 

The  solution  to  the  first  problem  can  be  derived  from  the 
Magnetic  (rih)  and  Electric  (ne)  Hertzian  potentials  [12]  which  have 

single  components  in  the  direction  n  nal  to  the  surface  of  the 
dielectric  (i.e.  the  z  direction).  The  equations  which  must  be  solved 
and  their  solutions  are  shown  in  Table  2,  3  and  4. 


Table  2 

Starting  Equations  for  the  Image  Antenna 
Problem  -  Slot  without  Clavin  Pins 
v2  nh+k2  rih=o  v2  ne+k2  ne=o 

E=-jwp  vxrih  E=k2ne+  vvne 

H=k2rih+  vvrih  H=j(oe  vxne 

n(x,y,z)«  [2ic]-2  J<r  j*  n(a,p,z)e'j(ax+Py)da5p. 

n(a,p,z)=j*9?  JsR  n(x,y,z)ej(ax+Py)3x3y. 


(18) 

(19) 

(20) 

(21) 

(22) 


f 
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Table  3 

Solution  to  the  Equations  in  Regions  i*M 


nhi=  Aj{cosh[yj(z-Zj.-|  )]+vJ/hjsinh[Yj(z-Zj_i)]]}  (23) 

nei=Bi{vFeiCOsh[Yj(z-zj.1))]+sinh[Yj(z-Zj.i)]}  (24) 

Yj2=a2+P2-kj2  (25) 

Ai={jxi-1/|xj}{Cj-1+4/hj_1Si.i}Aj-i  (26) 

Bj={Yi-1/Yi}{Cj.-|  +4/ei-i  Sj.-j  }Bj-i  (27) 

C  i==cos  h  [Yi(zj-z  j_  -) )  ]  (28) 

S  j=sinh[Yj(zj-Zj_i )]  (29) 

T  i-1  =Sj-i/Cj-i  (30) 

^hi-1={^hi-ZhiTj.i}/{Zhi-Tj-i'Fhj}  (31) 

^  ei-1  ={^ei"^eiT  j--j  }/{Z©i-Tj_ivIx  gj)  (32) 

ZhHH^i-lHYi-i/Yj}  (33) 

Zei  =  {ei-1/ei}{Yi/Yj_1 }  (34) 

Ai={(jP)Ex(Ss)-jaEy(Ss)}/{jcam(a2+|32)}  (35) 

B  i  =-{jaEx(Ss)+jpEy(Ss)}/{Y1  (a2+p2)}  (36) 

Table  4 

Solution  to  the  Equations  in  Regions  i=M 
OhM=  Am©xp[-ymz]  (37) 

neM=  BMexphrMz]  (38) 

AM={^M-1/^M}{CM-l+vFh(M-1)SM-l}exppTMZM-l]AM-1  (39) 

BM=-{YM-1/YM}{CM-l+vFe(M-1)SM-l}exP[YMzM-l]BM-1  (40) 

^hM-lHTM-lZhM+1}/{ZhM+T(M-i)}  (41) 

VJ/eM-1=*(TM-lZeM+1}/{ZeM+T(M-l)}  (42) 


In  the  Tables  above,  /$?  indicates  integration  over  ±c*  and  Ex(Ss) 
and  Ey(Ss)  are  the  Fourier  transforms  of  the  fields  in  the  slot. 

The  far  field  patterns  can  be  obtained  from  the  total 
transformed  electric  field  in  the  Mth  dielectric  region  given  by 

E(x,y,z)=[27i]'2/<R  Jgt{{-pa(Qh-«e}x+{a2Qh+P2Qe}y 

+{a2+p2Hjp)/YMnez}Ey(Ss)(a2+p2)-1ej(ax+Py)3a9p.  (43) 
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«h=n[{Ci.i+vPh(i-l)Sj-i}]exp[YMZM-l] 

^e=n[{Cj-i  +^6(1-1  )Sj-i}]exp[YMZM-l] 


(44) 

(45) 


The  use  of  asymptotic  techniques  as  outlined  by  Bleistein  [13]  to 
evaluate  this  integral  provides  an  explicit  formula  for  the  far  field 
pattern  of  the  configuration.  Inserting  the  standard  conversions 
from  rectangular  to  spherical  coordinates  into  (43),  an  equation  of 


the  following  form  results 

E(p,0,4>)|p^oO=±kMCOs(e)[2^p3-1g(ao,po)exp[jp4(ao>(5o)-j7i/2]  (46) 

where  a0  and  p0  are 

a0=±kMSin(e)cos(<i>)  (47) 

po=±kMsin(0)sin(<}>).  (48) 

The  function  §(a0/p0)  is  given  by 

§(aO/Po)=Ey(Ss)  jftesin(<)>)/cos(6)0+cos(<|>)fth<!>}  (49) 

Inserting  these  results  into  (46)  we  obtain  for  the  far  field 
E(p,0,<|))|p_^Oo--jkM[27cp]-lEy(Ss) 

|Qesin(<]))0+cos(<|))cos(0)flh(t)}exp[-jpk|yi].  (50) 


A  similar  route  can  be  taken  to  obtain  the  field  patterns  of  the 
monopole  radiators  when  they  are  excited  by  a  current  distribution 
with  a  single  component  in  the  z  direction.  For  this  type  of 
excitation  nh=0  so  only  ne  is  considered.  Table  5  shows  the  starting 

equations  for  the  regions  where  the  monopoles  do  exist. 
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Table  5 

Starting  Equations  for  the  Image  Antenna 
Problem  -  Monopole  Radiators 


« 


v2  ne+k2  ne=(1/jcoei)J 
03  =  jtoe  {jpnex-janey} 

I  =  {k2+a2/0Z2  }ne  z+jaa/aznex+jpa/azney. 

nei=(-1/ja)eiYi)  lo.h  coshfyi z<3{^et c°sh[yi z>]+ 

sinh[yi  z>]}Jz(z0)dz0 

z<  =  min(zo  ,2) 
z>=  max(zo  ,  z) 

Bi  =(-1/jtoeiY-|  )J0jh  cosh[yiz0]Jz(zo)dzo 
Bi=exp[yjZj-i]/(ja)eiYi) 

n[{Ci-i+'Pe(i-1)Si-l}]Jo,h  cosh[y1z0]Jz(zo)dzo 


(51) 

(52) 

(53) 


(54) 

(55) 

(56) 

(57) 


(58) 


In  Table  5,  Jo.h  indicates  the  integral  from  0  to  h  .  The  Bi  in  Table 

5  can  be  used  in  (45)  instead  of  (36). 

In  order  complete  the  solution,  the  moment  method  must  be 
applied  to  the  system  so  that  the  sum  of  the  fields  from  the  slot  and 
the  scattered  fields  from  the  Clavin  pins  is  zero  on  the  surface  of 
the  pin.  To  this  end  we  use  basis  functions  of  the  form 


Jz1(zoMaplfx(z)  (59) 

where 

Jz2(zoMa|32fx(z)  (60) 

fx(z)=Xn  Anx6(z-zon)  x-1,2  (61) 

fapl=fapexp(jpd)  (62) 

f<xp2«f<xpexp(-j0d).  (63) 

fap=(wp/2)2sinc(pwp/2)sinc(awp/2)  (64) 
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Here  d  is  the  distance  from  the  slot  to  the  Clavin  pin.  Galerkin's 
method  is  used  to  reduce  the  integral  equation  to  a  matrix  which  in 
turn  is  solved  for  the  amplitudes  Anx-  Due  to  the  shifting  properties 
of  the  Fourier  transform,  this  method  can  easily  be  extended  to 
include  all  of  the  slots  in  an  array  and  all  of  the  Clavin  pins,  yielding 
a  much  bigger  matrix,  but  all  the  coupling  between  the  slots  and 
Clavin  pins  would  then  be  accounted  for. 

Figure  14  shows  the  gain  enhancement  in  the  E-plane  as  a 
function  of  angle  for  an  uncompensated  slot  with  an  image  plate 
dielectric  constant  of  25  along  with  the  predicted  sidelobe 

reduction  due  to  the  Clavin  pins.  Also  shown  in  Figure  14  is  a  * 

comparison  between  this  theory  and  that  of  Immell  and  Sasser  [2]. 

As  can  be  seen  from  the  figures,  the  two  theories  yield  identical 
results. 

Experimental  evidence  shows  the  existence  of  waveguide  type 
propagation  modes  between  the  ground  plane  and  the  high  dielectric 
constant  image  plate.  In  order  for  these  modes  to  exist,  the  integral 
(21)  has  residuals  at  values  of  a  and  p  which  would  be  determined 
by  the  functions  in  (23)-(25)  .  The  residuals  could  be  obtained  via 
graphical  methods  or  some  other  numerical  technique  (however  the 
authors  have  not  yet  completed  this  task).  In  order  to  meet  the 
constraints  imposed  by  the  continuity  of  the  tangential  fields  at  the 
boundarys  between  each  of  the  dielectric  plates,  a  and  p  must  be 
equal  in  all  dielectric  regions.  Suppose  that  residual  values  for  a 
and  p  were  found  such  that  a2+p2  =k2res-  Since  we  are  integrating 
along  the  real  axis,  only  real  values  of  a  and  p  are  expected. 

Furthermore  in  order  for  the  resulting  solution  to  be  physically 
meaningful,  k2air<k2res  <k2max  where  kmax  is  the  wavenumber  of 
the  highest  dielectric  constant  in  the  system  (Er2  in  Figure  13). 

Since  the  wavenumbers  are  inversely  proportional  to  the  wavelength 
kair>^res  >^max-  Therefore,  if  erl  in  Figure  13  is  equal  to  one, 
wavequide  propagation  in  region  1  would  be  expected  to  occur,  as  it 
does  experimentally. 
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4.  DESIGN  AND  PERFORMANCE  OF  A  SPECIFIC  IMAGE  ARRAY 

A  high  gain,  low  sum  and  difference  sidelobe,  deep  null  antenna  was 
needed  for  a  Motorola  seeker  application.  Analysis  revealed  this  could  be 
achieved  with  an  optimum  image  array  of  9  elements  configured  within  a  5 
wavelength  diameter  aperture.  The  antenna  (Figure  15)  consists  of  an  image 
plate  with  a  dielectric  constant  of  16,  a  honeycomb  spacer,  a  metal  ground  plane 
and  a  triplate  assembly  containing  the  sum  and  difference  channel  stripline  feed 
network.  Fasteners  (not  shown)  and  bonding  films  were  used  to  produce  a 
rugged  design.  Coordination  with  the  receiver  and  gimbal  designers  was 
required  to  obtain  an  assembly  that  functions  properly  within  the  limited  volume 
available. 

The  number  of  slots  and  the  dielectric  constant  were  selected  using  the 
truncated  image-element  model  discussed  above.  The  antenna  schematic  is 
shown  in  Figure  16  and  the  stripline  feed  circuit  in  Figure  17.  Signals 
intercepted  at  each  radiating  slot  pass  along  a  triplate  section.  The  two  slots  in 
each  quadrant  are  combined  In  an  isolated  power  divider  and  these  outputs  are 
combined  in  a  ring  hybrid  comparator  to  form  sum,  E-plane  difference,  and  H- 
plane  difference  channels.  The  center  slot  is  used  only  in  the  sum  channel. 

The  reason  why  the  center  element  is  limited  to  the  sum  channel  is 
explained  by  noting  aperture  amplitude  (as  well  as  phase)  distributions  for 
optimum  sum  patterns  are  considerably  unlike  those  for  optimum  difference 
patterns  [14].  For  given  minor-lobe  requirements  and  a  fixed  aperture,  an 
optimum  sum-aperture  distribution  yields  maximum  gain,  and  an  optimum 
difference-aperture  distribution  provides  maximum  tracking  sensitivity.  Now  the 
analysis  of  Bayliss  [14]  is  limited  to  line-source  arrays.  However,  it  can  be 
generalized  to  circular  apertures.  The  mathematics  is  somewhat  complicated  for 
image  arrays.  Near-optimum  difference  patterns  can  be  realized  by  not  using 
elements  near  the  center  of  the  array,  while  near-optimum  sum  patterns  require 
a  center  element.  The  separable  sum  and  difference  channel  aperture 
distributions  provided  by  Motorola’s  design  are  thus  its  key  performance 
advantage.  Motorola  has  successfully  Incorporated  this  technique  into  several 
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image  array  designs. 

Because  of  high  mutual  coupling  among  the  elements,  active  matching  is 
necessary  to  achieve  the  performance  required.  To  illustrate  the  pattern 
distortion  which  can  arise  from  electromagnetic  coupling  between  the  center  and 
the  outer  elements,  the  array  of  nine  described  above  was  assembled  with  a 
separate  waveguide-to-coax  adapter  feeding  each  slot.  Waveguide  was  selected 
for  simplicity  of  testing.  The  center  element  waveguide  adapter  was  connected 
to  a  receiver  while  the  outer  element  adapters  were  terminated  in  phase  shifters 
and  shorts.  A  pattern  of  the  center  slot  was  taken  with  all  of  the  outer  element 
phase  shifters  set  to  the  same  position;  this  was  repeated  for  a  total  of  4  phase 
shifter  positions.  The  results  at  a  specific  frequency  are  shown  in  Figure  18, 
which  clearly  demonstrates  considerable  variation  in  gain  for  the  center  element. 

To  achieve  the  required  low  sidelobes  the  center  element  is  weighted 
heaviest.  Consequently,  substantial  loss  of  array  gain  resuits  if  the  center 
element  gain  is  reduced  by  excessive  coupling  to  the  outer  elements.  Additional 
reduction  to  array  gain  can  occur  if  this  coupling  also  distorts  element  phase 
relationships.  Furthermore,  to  achieve  the  required  sidelobes  in  the  presence 
of  high  coupling  requires  greater  amplitude  taper  than  without  coupling.  Thus 
active  impedance  matching  is  necessary  to  achieve  the  performance 
requirements. 

Considerations  for  a  good  design  are: 

e  Broadbanding  demands  active  impedance  matching  at  the  radiators. 

e  Active  impedance  matching  must  take  into  account  amplitude  weighting 
of  the  elements 

Table  6  presents  a  performance  summary  of  the  specific  image  array. 

TABLE  6 

ANTENNA  PERFORMANCE 


PARAMETER 

GOAL 

TYPICAL 

WORST 

GAIN 

CNTR  FRFQ 

>21.0  DBI 

+22.0 

+22.0 

FREQ  BAND 

>21.0  DBI 

+21.9 

+21.8 

SUM 

E-PL 

>21  DB 

23 

22 

SIDELOBES 

H-PL 

>21  DB 

24 

22 

DIFF 

E-PL 

>22  DB 

25 

24 

SIDELOBES 

H-PL 

>22  DB 

28 

27 

NULL  DEPTH 

>27  DB 

35 

33 
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RELATIVE  POWER  (dB) 


FIGURE  4.  COMPUTED  E-PLANE  PATTERNS  OF  AN  ISOLATED 
CLAY I N- COMPENSATED  IMAGE- SLOT  ANTENNA 


UNTRUNCATED 

(INFINITE) 

APERTURE 


APERTURE 

TRUNCATED 

AT  M=1S 
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9 (DEGREES) 

(b)  TRUNCATED  AT  M=1S 

IT  CURE  5.  THEORETICAL  E- PLANE  PATTERNS 
OP  ARRAY  OE  24  UNIFORMLY 
WEIGHTED  IDENTICAL  IMAGE 
ELEMENTS  WITH  CLAVIN  COMPENSATION 


IMAGE  RADIATORS 


FIGURE  6.  FIXTURE  AND  NETWORK  SCHEMATIC 

FOR  ACTIVE  RETURN- LOSS  AND  PATTERN 
MEASUREMENTS  OF  24-ELEMENT  ARRAYS 
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FIGURE  7.  2 4- ELEMENT  IMAGE  ARRAY  SLOT  LOCATIONS 


TABLE  1.  GAIN  AND  APERTURE  EFFICIENCY  OF  24-ELEMENT  CLAVIN- 
COMPENSATED  IMAGE  ARRAY 


GA 

1_ Lli! 

€-37 

FREQUENCY 

20L°o(^-) 

G* 

G-GA 

EFF 

G* 

G-GA 

EFF 

>o\  X  1 

(dBI) 

(dBI) 

(dBI) 

(%) 

(dBI) 

(dBI) 

(%) 

FL 

29.3 

28.4 

-.9 

81 

28.0 

-1.3 

74 

FC 

29.4 

29.0 

-.4 

91 

27.6 

-1.8 

66 

FH 

29.5 

28.8 

-.7 

85 

27.2 

-2.3 

59 

*  CALCULATED  FROM  (17)  USING  DATA  OF  FIGURES  10  AND  11. 


193 


COUPLING  (dB) 


-196  -84  84  196 

F  -  Fq  (MHz) 


FIGURE  8.  BORESIGHT  GAIN  OF  IMAGE  RADIATOR  NO.  1 
IN  24-ELEMENT  ARRAY 


S/a  -  SPACING  BETWEEN  ELEMENTS 


FIGURE  9.  RECORDED  COUPLING  BETWEEN 
PAIRS  OF  IMAGE  RADIATORS 
IN  24 -ELEMENT  ARRAY 
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PATTERNS  OE  ARRAY  OF  24  CLAVIN  COMPENSATED  IMAGE  ELEMENTS 


Clavin 

Pin 


K 


Clavin 

Pin 


Figure 


Enhanced 

Gain 

(dB) 


Figure 


Side  View 


O 


O 


y 


b)  Top  View  of  Slot  and  Clavin  Pins 


1.3,  Illustration  of  the  single  Clavin  pin  compensated  slot. 


(degrees ) 

E-Plane  with  pins 
E-Plane  w/o  pins 
E-Plane  Immell  5  Sasser 

14.  Predicted  E-plane  patterns  with  and  without  Clavin  pins 
and  comparison  of  this  theory  with  that  of  Immell  and 
S  asscr. 
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FIGURE  15.  NINE-ELEMENT  IMAGE  ARRAY  ASSEMBLY 


RELATIVE  POWER  (dB) 


9  (DEGREES) 


I;I CURE  18.  CENTER  IMAGE  -  ELEMENT  GAIN 
AS  A  FUNCTION  OF  OUTER 
ELEMENT  TERMINATION 
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DISTRIBUTED  BEAMSTEERING  CONTROL 
OF  MONOLITHIC  PHASED  ARRAYS 

S.F.  NATI ,  G.T.  COKINOS,  D.K.  LEVIS 
GENERAL  ELECTRIC  CO. 

ELECTRONICS  LABORATORY 
P.O.  BOX  4840 
SYRACUSE,  NEV  YORK  13221 

1.0  ABSTRACT 

As  satellite  communication  systems  move  to  higher  frequencies, 
a  great  deal  of  effort  is  being  directed  at  monolithic  phased 
arrays.1  These  arrays  combine  printed-circuit  antenna  technology 
with  monolithic  amplifiers  and  phase  shifters.  While  significant 
advances  in  these  areas  are  helping  to  make  monolithic  arrays  a 
reality,  array  control  still  presents  a  challenge. 

Phased  arrays  traditionally  have  been  controlled  by  a  single 
central  computer  that  calculates  a  phase  setting  for  each  array 
element.  This  architecture  has  worked  well  for  arrays  operating  at 
lover  frequencies,  having  a  small  number  of  elements  and  modest 
update  rates.  Distributed  control  of  phased  arrays  has  been  used 
for  radar  applications2  and  offers  a  solution  to  the  problems  of 
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array  control  for  high  frequency,  monolithic  arrays.  Since  each 
element,  or  small  group  of  elements,  has  its  own  control,  the 
bandwidth  and  number  of  connections  to  the  central  controller  is 
reduced.  In  addition,  the  computation  hardware  is  distributed, 
allowing  the  beamsteering  calculations  to  occur  in  parallel, 
providing  the  means  for  higher  update  rates.  Compensation  for 
variation  in  element  performance  become  possible  with  a  distributed 
control  scheme.  Distributed  control  architectures  also  offer  the 
potential  for  integration  with  the  microwave  components  which  would 
result  in  a  truly  monolithic  array. 

This  paper  will  discuss  the  design  of  a  distributed 
beamsteering  controller  (BSC)  for  high  frequency  phased  array 
applications.  The  controller  is  designed  using  AT&T's  digital  GaAs 
HEMT  technology.  A  review  of  the  architecture  and  design  for 
digital  HEMT  implementation  will  be  presented.  Areas  of  future 
development  will  also  be  discussed. 


2.0  PROBLEM 

The  use  of  higher  operating  frequencies  for  satellite 
communications  systems  is  reducing  the  area  available  for  the 
antenna  and  related  circuit  elements  such  as  amplifiers  and  phase 
shifters.  Figures  la  and  lb  show  a  block  diagram  of  a  typical 
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antenna  element.  Figure  la  shows  a  transmit  element  while  Figure 
lb  shows  a  corresponding  receive  element. 


Figure  la.  Transmit  Module 


Figure  1b.  Receive  Module 


Figure  L  shows  the  available  area  as  a  function  of  frequency, 
assuming  a  uniform  element  spacing  of  one  half  the  wavelength. 

This  decrease  in  element  spacing,  in  addition  to  reducing  circuit 
area,  also  limits  the  number  of  physical  connections  that  can  be 
made  to  the  antenna  element  and  its  control  circuitry.  Connections 


must  be  provided  for  RF  distribution,  power,  and  control. 


In  addition  to  the  smaller  element  spacing,  the  desire  for 
increased  array  performance  is  causing  the  number  of  elements  per 
array  to  increase.  This  combination  has  an  effect  on  the  number, 
speed  and  interconnect  density  of  the  required  control  signals.  An 
architecture  that  minimizes  the  number  of  connections  has  a 
distinct  advantage  as  operating  frequencies  increase. 

The  need  for  highly  manueverable  platforms  to  maintain 
communications  results  in  a  need  to  update  the  array  beam  pointing 
direction  at  higher  rates.  In  addition,  the  use  of  techniques  for 
improving  security,  such  as  spread  spectrum  communications,  place 
even  greater  demands  on  the  beamsteering  control  function. 
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3.0  CONTROL  ARCHITECTURES 


There  are  several  approaches  to  calculating  and  distributing 
the  phase  setting  for  each  element  in  a  phased  array.  The  most 
simplistic  approach  is  to  compute  the  phase  setting  for  each 
element  in  a  central  computer  and  apply  the  phase  settings  for  the 
entire  array  simultaneously.  This  approach  proves  to  be 
impractical  due  to  the  enormous  amount  of  interconnect  that  would 
be  needed  for  an  array  of  any  significant  size.  Additionally,  the 
amount  of  time  needed  to  perform  the  calculations  for  the  entire 
array  in  a  serial  fashion  would  severely  limit  the  update  rate  of 
the  array. 

Ar.  alternative  approach  is  to  compute  the  phase  settings  for 
each  row  and  column  and  add  the  results  at  each  element.  This 
approach  reduces  the  number  of  interconnects  required  but  is 
restricted  to  arrays  with  uniform  element  spacing.  The  limiting 
factor  in  this  approach  is  the  time  needed  for  calculation  and 
distribution  of  each  row  and  column  phase  setting.  As  the  array 
size  increases,  the  update  rate  will  become  limited.  This  approach 
also  makes  it  difficult  to  provide  compensation  for  element 
variat ions . 

A  fully  distributed  approach,  where  the  beamsteering 
calculation  is  performed  at  the  element  or  subarray  level,  has  many 
desirable  features.  By  broadcasting  the  desired  antenna  direction 
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cosines  from  a  central  array  controller  to  a  beamsteering 
controller  located  at  each  element,  the  number  of  interconnects  can 
be  reduced  significantly.  Each  beamstaering  controller  will 
calculate  the  phase  shifter  settings  in  parallel,  simultaneously, 
reducing  the  amount  of  time  needed  for  computation.  In  this  type 
of  approach,  the  communication  and  computation  time  is  basically 
independent  of  the  number  of  elements  in  the  array.  Variations  in 
element  performance  can  be  compensated  by  storing  a  correction 
memory  in  the  beamsteering  controller  and  applying  these 
corrections  on  an  element  by  element  basis. 

4.0  BSC  ARCHITECTURE 

The  overall  array  control  concept  is  shown  in  Figure  3.  Each 
beamsteering  controller  controls  a  subarray  of  the  entire  array. 

The  beamsteering  controller  receives  commands  from  a  central  array 
controller  and  based  on  the  command  performs  the  required 
operations  for  each  element  in  the  subarray.  These  commands  may  be 
broadcast  to  all  BSCs  simultaneously,  as  in  a  pointing  angle 
update,  or  may  be  addressed  to  an  individual  BSC,  as  in  a 
modification  to  a  compensation  table  entry.  The  BSC  is  also 
capable  of  performing  a  self  test  function  and  reporting  it's 
status  to  the  array  controller. 
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Figure  3.  Control  Concept 

A  timeline  of  a  typical  beam  update  sequence  is  shown  in 
Figure  4.  The  sequence  begins  with  the  array  controller 
broadcasting  to  each  BSC  the  direction  cosines  of  the  overall  array 
pointing  angle.  Upon  receiving  this  command  each  BSC  will 
calculate  the  required  phase  shifter  settings  for  each  element  of 
its  subarray.  These  settings  will  be  stored  within  the  BSC  until 
an  update  command  is  received.  Next,  the  array  controller  selects 
one  BSC  for  testing.  The  selected  BSC  receives  the  test  command 
and  compares  the  results  of  the  previously  calculated  phase  shifter 
settings  with  those  supplied  by  the  array  controller.  The  array 
controller  then  broadcasts  the  update  command  to  all  the  BSCs  and 
the  new  phase  shifter  settings  are  applied  to  each  phase  shifter  in 
the  array.  The  BSC  that  was  selected  for  testing  then  transmits 
the  test  results  to  the  array  controller. 
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Figure  4.  Beam  Steering  Timeline 


Figure  5  shows  a  functional  block  diagram  of  the  beamsteering 
controller.  The  main  functions  include  a  single  wire  interface  to 
the  main  array  computer,  command  decoding  logic,  beamsteeering 
calculation  logic,  a  compensation  table,  built  in  test  logic,  and 
timing  and  control. 


to  phase  shtttks 


Figure  5.  BSC  Functional  Block  Diagram 
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5.0  BSC  DEVELOPMENT 


The  number  of  elements  that  each  BSC  will  control  is  based  on 
tradeoffs  with  regard  to  the  required  system  update  rate,  the 
processing  capability  of  each  BSC,  along  with  area,  power  and 
interconnect  considerations.  It  was  decided,  based  on  the  above 
tradeoffs,  to  design  the  BSC  to  control  up  to  sixteen  elements. 

The  BSC  receives  its  c  mmands  from  the  array  controller  over  a 
single  bidirectional  interface.  The  bidirectional  decoder  that  is 
used  to  interface  between  the  array  controller  and  the  BSCs 
provides  the  capability  for  broadcast  or  individual  addressing. 

The  BSC  instruction  set  is  shown  in  Figure  6.  The  CALCULATE, 
UPDATE,  and  TEST  commands  are  used  during  each  beam  update,  while 
the  LOAD  COMPENSATION  TABLE,  and  the  LOAD  POSITION  VECTOR  commands 
are  used  during  array  initialization  or  when  changes  in  parameters, 
such  as  temperature,  require  the  values  in  the  compensation  table 
to  be  adjusted.  Each  instruction  is  a  fixed  length  and  is  divided 
into  a  three  bit  command  field  and  a  twenty-four  bit  data  field.  A 
single  bit  is  provided  for  parity. 
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Figure  6.  BSC  Instruction  Set 


The  CALCULATE  command  contains  the  direction  cosines  for  the 
desired  antenna  pointing  direction.  This  command  is  normally 
broadcast  to  each  BSC  simultaneously.  The  BSC  calculates  the  phase 
shifter  settings  for  its  subarray  based  on  these  direction  cosines 
and  the  position  of  the  BSC  in  the  array. 

The  UPDATE  command  is  also  normally  a  broadcast  command.  The 
BSC  will  update  its  phase  shifter  settings  upon  receipt  of  this 
command.  The  data  that  has  been  calculated  during  the  previous 
CALCULATE  command  is  transferred  to  the  phase  shifter  control  lines 
for  each  element  in  the  subarray. 
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The  TEST  command  is  used  to  perform  a  functional  test  on  one 
BSC  during  each  beam  update  period.  The  data  that  has  been 
computed  as  a  result  of  the  CALCULATE  command  is  compared  on  a  bit 
by  bit  basis  with  the  data  computed  by  the  array  controller  for  the 
designated  BSC.  The  BSC  will  report  the  result  of  the  comparison 
back  to  the  array  controller. 

The  LOAD  COMPENSATION  TABLE  command  contains  the  information 
to  be  loaded  into  the  compensation  table.  This  provides  the 
ability  to  adjust  the  phase  shifter  settings  for  each  subarray 
based  on  data  obtained  during  an  initial  calibration  sequence. 

The  LOAD  POSITION  VECTOR  command  gives  each  BSC  its  position 
within  the  array.  This  information  is  used  during  the  CALCULATE 
command . 

5.1  CALCULATIONS  PERFORMED 

The  beamsteering  calculation  that  is  implemented  in  the  BSC  is 
shown  in  Figure  7.  The  BSC  assumes  that  the  array  is  planar  and 
has  uniform  element  spacing  within  the  subarray.  Although  the  BSC 
could  be  extended  to  the  more  general  case,  these  assumptions  allow 
the  BSC  to  use  two  multiplies  per  oeam  update,  minimizing  the 
amount  of  time  needed  for  computation.  In  addition,  the  internal 
logic  for  the  beamsteering  calculation  is  implemented  using  a 
pipelined  bit  serial  approach. 
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Figure  7.  Phase  Shift  Command  Calculation 


This  approach  is  highly  efficient  in  area  and  power.  A 
compensation  table  provides  the  ability  to  adjust  to  varying  phase 
shifter  performance  between  subacrayc.  i igure  8  shows  an  example 
of  an  imperfect  phase  shifter  and  the  benefits  that  can  be  obtained 
through  the  use  of  the  compensation  table. 
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Figure  8.  Effects  of  Compensation 
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5.2  FAULT  MONITORING 


The  approach  that  is  used  for  fault  monitoring  at  the  subarray 
level  is  to  test  a  single  BSC  during  each  beam  update.  This  allows 
the  entire  array  of  BSCs  to  be  tested  in  a  round-robin  fashion.  A 
BSC  is  tested  by  comparing  its  calculated  phase  shifter  data  with 
data  that  has  been  calculated  by  the  array  controller  for  the  same 
pointing  angle.  This  will  test  a  significant  portion  of  the  BSC, 
including  the  array  interface,  command  decode  logic,  calculation 
logic,  position  registers,  and  compensation  table.  In  addition, 
parity  and  message  format  checks  are  performed  on  the  incoming 
command  sequence.  These  errors,  along  with  a  test  status  are 
reported  back  to  the  array  controller  at  the  conclusion  of  the  test 
sequence . 

6.0  FUTURE  DEVELOPMENT 

The  ultimate  goal  of  incorporating  the  BSC  with  the  phase 
shifters,  amplifiers  and  antenna  in  a  monolithic  array  requires 
that  the  BSC  be  capable  of  being  integrated  into  a  microwave 
compatible  process.  For  analog  components  at  EHF,  it  is 
anticipated  that  a  gallium  arsenide  HEMT  process  will  be  used. 

With  this  in  mind,  the  BSC  was  designed  to  be  implemented  in  AT&T's 
digital  GaAs  HEMT  process.  This  process,  or  similar  processes, 
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although  not  intended  for  implementing  analog  components  at  EHF 
frequencies  at  the  present  time,  could  conceivably  be  modified  to 
do  so.  A  process  using  selective  molecular  beam  epitaxy,  would 
provide  the  capability  to  produce  devices  optimized  for  analog  or 
digital  functions  on  the  same  substrate.  This  would  be  the  next 
step  in  producing  a  truly  monolithic  array. 
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ABSTRACT 

An  important  concern  in  the  design  of  phased  array  antennas 
is  the  effect  of  mutual  coupling  between  elements  on  the  array 
performance.  For  large  periodic  arrays,  a  common  approach  is  to 
use  an  infinite  array  formulation.  This  approach  does  not  include 
edge  effects  and  is  valid  only  for  a  uniform,  periodic  excitation 
of  the  array.  Alternatively,  the  element-by-element  approach  can 
be  used,  but  this  may  be  impractical  for  large  arrays.  In  this 
paper  a  method  for  the  analysis  of  finite,  periodic  arrays  of 
arbitrary  size  and  excitation  is  discussed.  The  method  is  based 
upon  an  efficient  solution  of  the  extended  infinite  periodic 
array.  In  essence,  the  infinite  array  is  excited  with  each  of  the 
coefficients  of  the  discrete  Fourier  transform  of  the  desired 
finite  array  excitation.  The  finite  array  behavior  is  then 
obtained  by  superposition.  After  initial  characterization  of  the 
infinite  array,  finite  arrays  with  arbitrary  excitation  and 

number  of  elements  can  be  solved  by  simple  FFT  operations  on  the 
excitation  and  impedance  matrices.  Finally,  by  extracting  an 


approximate  open-circuit  impedance  matrix  for  the  array,  the 
solution  for  free  or  forced  voltage  excitation  is  also  obtained. 


1.0  INTRODUCTION 

An  important  criteria  in  the  design  of  any  phased  array 
antenna  is  the  behavior  of  the  active  element  impedance  as  a 

function  of  beam  scan  angle.  For  large  arrays  with  a  single  beam, 
this  can  often  be  sufficiently  characterized  either 
experimentally,  by  constructing  a  small  section  of  the  array,  or 
theoretically,  using  infinite  array  theory  [1]. 

This  paper  is  concerned  with  the  analysis  of  finite  but 
periodic  phased  arrays  of  arbitrary  size  and  excitation,  as  might 

be  used  for  steering  multiple  beams  and  nulls.  The  analysis  of 
this  type  of  array  is  more  difficult,  particularly  for  relatively 
large  arrays  where  the  element-by-element  approach  can  become 
impractical  due  to  the  large  size  of  the  impedance  matrices. 

The  method  presented  here  for  the  analysis  of  finite  arrays 

is  based  upon  an  efficient  solution  for  the  active  array  impedance 
of  the  infinite  array  created  by  continuation  of  the  finite 
periodic  array.  This  approach  allows  arrays  of  essentially 
arbitrary  size  to  be  analyzed.  Importantly,  storage  requirements 
are  proportional  to  either  the  square  of  the  number  of  basis 
functions  on  each  element  or  the  total  number  of  elements  squared 
(whichever  is  greater,)  not  the  product  of  the  two,  as  in  the 

element  by  element  approach. 
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The  method  is  simple  to  implement  because  it  requires  only 
the  computation  of  the  impedance  of  an  infinite  array  of  the  same 
elements.  A  matrix  of  infinite  array  scan  impedances  is  computed 
over  one  period  of  scan  (visible  and  invisible)  space.  This 
matrix  is  multiplied  by  the  discrete  Fourier  transform  of  the 
desired  current  excitation  to  obtain  the  terminal  voltages  in  the 
finite  array,  from  which  the  active  element  impedances  are 
obtained.  In  addition,  the  element  impedances  can  be  obtained  for 
the  case  of  voltage  excitation,  either  free  or  forced  [1]. 

A  brief  formulation  of  the  method  is  given  in  Section  2. 
Computational  efficiency  is  briefly  discussed  in  Section  3,  and 
some  examples  are  considered  in  Section  4.  Comparisons  are  made  to 
element-by-element  results  for  arrays  in  free  space  and  microstrip 
dipole  arrays. 

2.0  THEORY 

It  is  well  known  that  the  impedance  of  a  single  element  can 
be  obtained  from  integration  of  the  impedance  of  an  infinite  array 
of  the  same  elements  [2].  Additionally,  Ishimaru  and  co-workers 
have  shown  that  the  impedance  of  a  finite  array  of  current  excited 
elements  can  be  expressed  as  a  convolution  integral  of  the 

Fourier  transform  of  the  array  aperture  current  distribution  and 
the  infinite  array  active  impedance  as  a  function  of  scan  angle. 
A  limitation  in  that  work  was  slow  numerical  convergence  in  the 
convolution  integrals.  The  formulation  here  is  entirely  discrete, 
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and  may  be  more  straightforward  conceptually,  but  the  concepts 
are  very  similar  to  those  presented  by  Ishimaru  [3].  The  technique 
is  founded  in  the  basic  principles  of  the  array  scanning  method 
presented  by  Munk  and  Burrell  [2] 

Of  particular  interest  here  is  a  finite  array  of  microstrip 
dipoles,  a  section  of  which  is  illustrated  in  Fig.  1.  The  array 

is  periodic,  represented  by  lattice  vectors  St  and  S2>  The 

dipoles  have  width  a  and  length  b,  and  are  printed  on  a  substrate 
of  thickness  t  and  relative  permittivity  For  generality,  the 

solution  allows  a  superstrate  layer  of  arbitrary  thickness  and 
relative  permittivity,  and  loss  tangents  in  both  the  superstrate 

and  substrate.  The  method  under  consideration  is  based  upon  the 
formulation  of  the  active  array  impedance  for  a  infinite,  periodic 
extension  of  such  an  array. 

2.1  Infinite  Array  Impedance  vs.  Scan  angle:  Z  0?Q). 

The  active  impedance  of  the  infinite  array  can  be  found  using 
well  known  methods  [1,4, 5, 6].  In  this  work,  a  moment  method 
solution  was  used  which  uses  subsectional  basis  functions  with  a 
linear  dependence  (triangle)  in  the  y  direction  (length)  and  an 
edge  condition  in  the  x  direction  (width).  The  edge  condition  is 
not  necessary  for  very  thin  elements,  but  can  give  improved 

results  for  wide  strips  when  no  x  directed  basis  functions  are 

employed.  The  details  of  the  moment  method  solution  used  here  can 
be  found  elsewhere  [5,6] 
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For  purposes  of  discussion  only,  the  form  of  the  infinite 
array  array  impedance  that  obtains  from  a  single  mode 
approximation  to  the  current  is  considered  here: 


z  =  x  F  F  G  (K  )  J  (k  )  (i) 

in  A  L  L  yy  mo  y  mn 
m  n 

^  AA 

A  is  the  area  of  one  unit  cell  in  the  array,  G  is  the  yy 

yy 

component  of  the  spectral-domain  dyadic  periodic  Green’s  function 

A# 

for  the  problem  under  consideration,  and  J  is  the  Fourier 
transform  of  the  single  mode  current  approximation.  The 

summations  are  from  -oo  to  +  oo.  Equation  1  is  a  summation  of 
Floquet  modes  defined  by  the  vector  K  , 

mn 

-»  ->  -> 

K  -  fi  +  fin  (2a) 

mn  mn  0  ' 

where: 

A.  =  +  X  and  0„  =  +  "X  <2b> 


The  vector  represents  the  vector  from  the  (0,0)  Floquet  mode 

to  the  (m,n)  Floquet  mode.  and  fi2  are  lattice  vectors  for  the 

reciprocal  or  periodic  lattice,  as  shown  in  Fig.  2  for  an 

arbitrary  array  periodicity.  For  a  rectangular  array,  fi  ~  In  IS 

and  fi2  =  2 nlS2.  fiQ  is  an  offset  of  each  Floquet  mode  to  allow 
beam  scanning,  as  discussed  below.  As  indicated  by  Eq.  1,  the 
active  array  impedance  is  obtained  by  summing  Floquet  modes  over 
the  reciprocal  lattice;  hence  the  active  array  impedance  is 
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periodic  over  the  reciprocal  lattice  just  as  the  array  structure 
is  periodic  over  the  spatial  lattice. 

To  scan  the  main  beam  of  the  array,  the  voltage  source  at 
element  (m,n),  denoted  V(Sffln)  is  adjusted  to  give 

“4 

V(S  )  -  V(S  )e~^o  -(mS1  +  nS2).  The  constants  mA  and  n  are 

mn  00  0  0 

essentially  the  direction  cosines,  which  for  a  rectangular  array 
become  (following  the  conventions  in  [6]): 

mQ  =  -S1  sin0  cos <f>  (3a) 

n  =  -S  sin#  sin^  (3b) 

0  2 

where  S  =  |S  [.  One  can  alternatively  excite  element  (m,n)  with 

a  current  source  I(Smn)  =  I(SQO)e^o  (mS,  +  nS2)  The  active  array 

impedance  is  the  same  in  either  case  for  an  infinite  array. 

For  notational  convenience  the  active  array  impedance  of  the 
infinite  array  is  denoted  Z(/?0),  to  indicate  that  it  is  a  function 

of  fi  ,  and  that  it  can  be  considered  a  spectral  domain  quantity. 

O 

2.2  Finite  Array  Impedance :  Arbitrary  Current  Excitation. 

There  are  many  applications  where  one  would  like  to  define 

excitation  coefficients  for  each  element  of  a  finite  array.  For 

an  N  x  N  element  array,  this  defmes  a  discrete  excitation 

1 

function  Ieil(S  ),  which  gives  the  excitation  current  at  each  of 
the  N  N  elements  of  the  array,  and  where  S  =  mS,  +  nS  .  Having 
defined  an  array  geometry  and  the  desired  excitation  currents,  one 
can  analyze  the  array  using  the  standard  element-by-element 
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approach.  However,  one  can  proceed  in  an  alternative  manner  as 


follows.  If  the  x  Ny  array  under  consideration  is  periodic, 

one  can  consider  the  hypothetical  extension  of  the  array  to  an 
infinite  structure.  The  goal  is  to  synthesize  the  excitation  of 
the  finite  array  by  a  sum  of  appropriately  weighted,  uniform 
periodic  excitations  of  the  infinite  structure.  Those  elements 

with  zero  current  can  be  considered  open-circuited.  For  resonant 
length  elements,  these  elements  have  negligible  effect,  allowing 
the  finite  to  be  approximately  recovered.  For  elements  for 

which  p*  circuiting  does  not  render  them  negligible,  the  current 
would  need  to  be  forced  to  zero  at  additional  points. 

Unfortunately  any  finite  array  excitation  is  not 

band-limited,  and  cannot  be  synthesized  exactly  with  a  finite  set 
of  periodic  excitations  of  the  infinite  array.  However,  one  can 

obtain  approximate  results  as  follows.  Given  an  excitation 
function  for  the  finite  array,  one  computes  the  discrete  Fourier 
transform  (DFT)  using  LK  samples  of  the  array  excitation: 


I  ex(j3*n) 


DFT 


ns]k> 


(4) 


where  the  DFT  pair  for  a  discrete,  periodic  function  defined  by 

■) 

LxK  samples  on  the  2-D  lattice  defined  by  given  by  [7]: 
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L  -  1  K  -  1 


f  (/?““)  =  DFT  |f(Su)j  =  l  [  f(Slk)  e'j27l(ml/L  +  nk/K)  (5a) 


1=0  k  =0 


L-l  K-l 


f(Smn)  =  IDFT  ^  [  [  f(/?““)  <J27r(ml/L  +  nk/K)(5b) 


1=0  k  =*0 


where 


0™  =  m  p{IL  +  n  fi2/K  0  m  ^  L-l,  0  £  n  £  K-l  (6) 


and  IDFT  denotes  inverse  discrete  Fourier  transform. 


LK,  the 


number  of  sampling  points,  will  typically  be  much  larger  than 
N  N  ,  the  number  of  elements  in  the  array.  Thus  Ie*(S,  J  will  be 

x  y  Ik 

zero  for  all  of  the  LK  samples  except  for  the  N  ^  which 
correspond  to  elements  of  the  finite  array. 

Once  the  DFT  of  the  excitation  currents  is  calculated,  the 
DFT  of  the  terminal  voltages  in  the  finite  array  can  be  determined 


from: 


A/  ->  A/  -)  A/  -> 

V(Pmn)  =  Z(Pma) 


The  element  voltages  for  the  N  x  N  array,  V(S  ),  can  be 

x  y  mn 

obtained  by  inverse  DFT: 


V(S  )  =  IDFT  (  V(/Tn)  \ 

mn  I  0  1 


This  IDFT  returns  LK  unique  element  voltages,  which  include  the 
voltages  for  the  original  x  Ny  array.  All  other  coefficients 


should  be  approximately  zero.  Finally,  the  input  impedance  of 
the  element  (m,n)  in  the  finite  array  is  obtained  from: 

-*  V(S  ) 

Z(S  )  =  - (9) 

Iex(S  ) 

mo 

The  approach  described  here  can  be  viewed  as  a  superposition 
technique.  The  DFT  of  the  desired  array  excitation  gives  a  finite 
set  of  uniform  periodic  excitations  of  the  infinite  array.  One 
must  keep  in  mind  that  the  DFT  of  the  excitation  current  is  a 
periodic  function,  with  periods  fi  and  fiy  It  is  defined  at  the 
LK  sample  points  defined  by  the  vectors  This  set  of  LK 

excitations,  when  summed,  recovers  a  periodic  version  of  the 
desired  excitation  function,  periodic  over  an  array  of  LxK 
elements.  In  essence,  it  creates  an  "array  of  arrays";  an 
infinite,  periodic  array  consisting  of  repetitions  of  the  desired 
finite  array.  It  is  necessary  that  this  second  periodicity  be 
large  enough  so  that  these  repeated  arrays  would  not  affect  the 
current  distribution  on  the  original  x  Ny  element  array.  This 
problem  will  be  examined  by  comparison  of  results  with  the 
element-by-element  approach. 


2.3  Open  Circuit  Impedance  Matrix  and  Mutual  Impedances . 

It  is  of  interest  to  obtain  the  open  circuit  impedance  matrix 
elements  for  the  finite  array.  For  an  x  Ny  array,  element 
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(p,q)  of  the  N  N  by  N  N  open  circuit  impedance  matrix  of  the 

x  y  x  y 

array  is  as  given  by: 

V 

Zoc  -  — (10) 

N  I  I  =0,  q  p 

a 

where  p  and  q  denote  a  single  integer  denumeration  of  the  array 
elements.  For  elements  which  have  essentially  no  current  when  open 
circuited,  such  as  half-wave  dipoles,  one  can  obtain  the  open 
circuit  matrix  for  a  finite  array  by  choosing  an  excitation  matrix 


-»  1 ,  m  =  0 ,  n  =  0 

rx(s  )  = 

mo  0,  otherwise 


(11) 


In  this  case,  the  response  at  the  unexcited  elements  of  the  array 
are  given  by: 

V0C(S  )  =  IDFT  (  ZOO  1  (12) 

mn  I  0  1 

The  open  circuit  impedance  matrix  elements,  (the  mutual 
impedance  between  any  two  elements  in  the  array  with  all  other 
elements  open  circuited)  can  be  obtained  from  Eq.  9, 

Z  fS  -  SJ  =  V0C(S,  it  i )  (13) 

mv  mn  kr  |  m-k  |  ,  ]  n-1 1  '  ’ 


Note  that  there  are  only  unique  values  obtained  by  Eq.  13  for 

an  x  array.  Equation  12  returns  an  approximation  to  the 
open  circuit  element  voltage  the  finite  array;  the  location  in 
the  array  is  not  taken  into  account.  Thus,  while  edge  effects  are 
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included  in  the  excitation  matrix,  they  are  not  included  in  the 
open  circuit  matrix.  However,  as  stated  earlier,  this  is  a  good 

approximation  for  many  practical  elements. 

2.4  Finite  Array  Impedance:  Arbitrary  Voltage  Excitation. 

For  a  voltage  excited  array,  direct  application  of  the  above 
formulation  leads  to  impedance  with  all  unexcited  elements  short 

circuited  (or  terminated  in  the  generator  impedance).  However, 
the  input  impedance  for  the  elements  under  voltage  excitation 
with  all  unused  elements  open  circuited  can  be  obtained  using  this 
method  by  using  the  open-circuit  impedance  matrix.  Once  the 
open-circuit  impedance  matrix  is  found,  the  element  currents 
under  voltage  excitation  are  determined  by: 

[I]  =  <[Z0C]  +  [ZT])'[V]  (14) 

where  [Z^J  is  a  diagonal  matrix  which  contains  the  generator 
impedances  and  any  resonating  elements.  This  gives  the  element 

currents  under  "free"  excitation.  Setting  [ZT]  =  0  gives  the 
"forced"  voltage  excitation,  according  to  the  definition  in  [1]. 
In  general,  the  inverse  of  the  Z  matrix  is  not  found,  but  the 
matrix  equation  is  solved  for  the  unknown  currents.  The  input 

impedance  of  the  m,nth  element  is  obtained  from 

•) 

Z(S  )  =  V  /I  -  ZT  (16) 

ran  ran  ran  ran 

where  the  subscripts  m,n  indicate  the  dependence  of  all  quantities 
on  their  location  in  the  array,  given  by  S 
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3.0  NUMERICAL  CONSIDERATIONS 

Clearly  the  most  time  consuming  aspect  of  the  solution  is  the 
filling  of  the  Z(/?“D)  matrix.  Filling  this  matrix  requires 
sampling  the  infinite  array  active  impedance  over  one  period  of 
the  spectral  lattice.  For  a  rectangular  array  of  center  fed 

dipole,  symmetry  reduces  the  number  of  samples  by  a  factor  of 
four.  Even  with  symmetry,  a  256x256  DFT  requires  the  computation 
of  16,641  samples.  Fortunately,  the  impedance  matrix  is  filled 
only  once  for  a  given  (infinite)  array.  Having  filled  the 

impedance  matrix,  the  performance  of  any  subset  of  the  infinite 
array  can  be  evaluated  for  arbitrary  excitation  using  the  methods 
of  Sec.  2.  Still,  sampling  the  impedance  matrix  sufficiently  to 
allow  the  analysis  of  microstrip  arrays  can  be  quite  time 
consuming  if  not  efficiently  done.  Two  methods  have  been  used  to 
reduce  the  fill  time  significantly. 

First,  the  time  required  to  compute  each  sample  of  the 

impedance  matrix  can  be  reduced  significantly  by  applying 
acceleration  techniques  to  the  summations  of  Floquet  modes 
expansions.  Many  such  summations  can  be  required  to  compute  the 
moment  method  solution  for  the  array  impedance  at  each  scan  angle. 
The  method  used  here  requires  obtaining  an  asymptotic  expansion 
for  the  generalized  impedance  matrix  elements.  The  asymptotic 

form  is  subtracted  from  the  original  expression,  leaving  a  rapidly 
convergent  summation.  The  asymptotic  terms  can  be  directly 

accelerated  using  Poisson’s  summation  formula.  Details  of  the 
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method  are  given  in  [6]. 

The  use  of  this  method  results  in  a  rapidly  convergent 
series,  typically  giving  a  reduction  in  computation  time  of 

approximately  two  orders  of  magnitude  over  a  direct,  unaccelerated 
summation.  This  has  allowed  the  routine  use  of  256x256  or  512x512 
impedance  matrices.  High  sampling  densities  are  needed  for 
microstrip  problems,  where  surface  waves  can  lead  to  singularities 
in  the  impedance  vs.  scan  angle  surface. 

Second,  an  adaptive  sampling  routine,  similar  to  adaptive 

integration,  combined  with  a  birational  interpolation  technique, 
has  been  used  to  reduce  the  number  of  scan  angles  at  which  the 

array  impedance  is  calculated.  With  this  routine,  the  impedance 
surface  is  sampled  over  a  grid  of  17x17  points.  Next,  the 
impedance  is  calculated  using  a  birational  interpolation  technique 
at  the  center  of  each  of  the  16x16  subsections  defined  by  the 

initial  gridding.  If  the  interpolated  result  is  sufficiently 
close  to  the  actual  computed  value  at  that  point,  that  subsection 
is  designated  as  complete,  and  the  remaining  values  filled  by 

interpolation.  If  the  interpolated  value  does  not  meet  tolerance, 
that  subsection  is  subdivided,  and  the  testing  procedure 
continued.  A  birational  interpolation  routine  was  chosen  because 
of  its  performance  near  singularities  [8].  This  technique  is 
continued  over  each  of  the  16x16  subsections,  reducing  the  grid 
size  in  each  subsection  only  as  far  as  necessary  to  meet  a 
predefined  tolerance. 


229 


The  time  savings  obtained  by  the  adaptive  filling  algorithm 
could  be  improved  by  more  sophisticated  interpolation  algorithms. 

The  time  savings  is  also  highly  dependent  on  the  nature  of  the 
active  array  scan  impedance.  For  arrays  in  free  space  or  over 

ground  planes,  where  the  scan  impedance  is  relatively  smoothly 
varying,  it  can  provide  significant  savings.  For  microstrip 
arrays,  it  can  be  less  effective.  Table  1  shows  the  number  of 

impedance  samples  that  were  computed  (as  opposed  to  interpolated) 

and  the  relative  time  savings  compared  to  computing  each  sample 
required  to  fill  the  impedance  matrix  for  a  microstrip  dipole 
array.  The  array  had  element  width  a  =  0.0UQ,  length  b  = 

0.39Aq,  array  spacing  =  S2  =  0.5AQ,  substrate  thickness  t  = 

0.19Aq  and  =  2.55.  In  addition,  a  loss  tangent  of  d  =  0.005  in 
the  substrate  was  used  to  eliminate  singularities.  As  would  be 
expected,  the  effectiveness  of  the  method  improves  with  finer 
sample  grids. 

Using  these  methods,  fill  times  of  5  minutes  for  a  256x256 
DFT  are  typical  on  a  DECstation  5000. 

4.0  RESULTS  AND  DISCUS&ION 

The  method  has  been  tested  by  comparing  to  results  in  the 
literature  and  by  comparing  to  results  obtained  using  a  standard 
spatial  domain  moment-method  code  employing  subsectional 
sinusoidal  basis  and  testing  functions.  Three  examples  will  be 
considered. 
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TABLE  1 


Reduction  in  computation  time  resulting  from  adaptive  sampling. 
Savings  is  expressed  as  the  ratio  of  the  total  sample  points  to 
the  actual  number  of  calls  to  the  impedance  evaluation  routine  as 
a  function  of  the  DFT  resolution.  Example  is  for  a  microstrip 
dipole  array,  with  spacing  of  0.5A0  in  x  and  y,  element  length  b 
=  0.39Aq,  width  a  =  0.01AQ,  on  a  substrate  with  t  =  0.19Ao, 
£r  =  2.55,  <5  =  0.005. 


DFT  Total  Actual  Savings 

Resolution  Sample  Points  No.  Calls.  Factor 


64 

1,089 

455 

2.4 

128 

4,225 

1,067 

4.0 

256 

16,641 

2,419 

6.7 

512 

66,049 

5,098 

13.0 
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4.1  Mutual  Impedance  Calculations. 

The  mutual  impedance  between  parallel  and  collinear 
microstrip  dipoles  was  computed  using  the  present  method  and 
compared  to  results  computed  by  Alexopolous  and  Rana  [9].  The 
method  described  here  is  denoted  the  DFT  method  in  the  Figures. 
The  impedances  were  computed  for  dipoles  with  length  1  =  0.333Afl, 
S  =  S  =  0.5  A  ,  on  a  substrate  with  e  =  3.25  and  thickness 

12  0  r 

of  t  =  0.1016Aq.  No  superstrate  layer  was  included.  A  loss 
tangent  of  0.005  was  included  in  the  substrate  layer  to  eliminate 

As 

the  singularity  in  the  impedance  surface  Z(/?Q).  This  small  loss 
tangent  was  found  to  have  little  effect  on  the  numerical  results. 
Results  for  the  mutual  impedances  are  shown  in  Figs.  3a  and  3b, 
with  good  agreement. 

4.2  Finite  Arrays  of  Microstrip  Dipoles  with  Voltage  Excitation. 
Next,  the  impedance  of  a  19x19  array  of  microstrip  dipoles 

was  computed  versus  scan  angle  and  compared  to  results  calculated 
in  [10].  For  this  case,  1  =  0.39Aq,  Sj  =  S2  =  0.5Ad  and  the 
substrate  had  =  2.55  and  t  =  0.19Aq.  Again,  a  loss  tangent  of 
0.005  was  included  in  the  substrate  layer,  and  a  256x256  point  DFT 
was  used.  The  results  in  [10]  were  calculated  using  the 
elemcnt-by-element  approach,  using  a  sinusoidal  approximation  to 
the  current  on  each  element  (1  entire-domain  mode  was  used  on  each 
element).  The  results  from  the  present  method,  also  using  a 
single  mode  approximation  to  the  element  current,  and  the 
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element-by-element  method  approach  are  compared  in  Figs.  4a  and 
4b.  The  reflection  coefficient  of  the  center  element  in  each  row 

along  the  E-plane  is  shown  for  scan  angles  of  0  and  30  degrees 

(Fig.  4a)  and  45  and  60  degrees  (Fig.  4b). 

It  was  found  in  [10]  that  little  change  in  the  reflection 
coefficient  vs.  scan  angle  was  seen  when  additional  modes  were 
used  to  represent  the  current  on  the  elements.  This  indicates 
that  the  behavior  of  the  array  versus  scan  angle  is  well 

represented  using  a  single  mode  for  the  case  of  approximately 

resonant  length  dipoles.  However,  it  should  be  noted  the 
significant  error  in  the  actual  impedance  values  can  results  from 
this  approximation.  Using  a  single  mode  the  input  impedance  of  a 
single  isolated  dipole  was  found  in  [10]  to  be  Z  =  105  +  j72  Q. 
This  compares  well  to  the  value  of  103  4-  j65  Q  found  using  the  DFT 
method  and  a  single  mode.  However,  using  19  subsectional  modes 
with  the  DFT  method  resulted  in  an  impedance  of  154  +  j38  Q,  which 
compares  well  to  the  result  of  154  +  j31  Q  obtained  by  Michalski 
and  Zheng  using  a  subsectional  basis  function  moment  method 
solution  with  51  unknowns  [11]. 

4.3  19  Element  Array  with  Two  Beams. 

As  a  final  example,  a  19  element  linear  array  of  y  directed 
0.5A  dipoles  in  free  space  was  considered.  The  dipole  elements 
were  spaced  0.5A  apart  along  the  x  axis.  Initially,  a 
Tschebyscheff  distribution  with  30  dB  sidelobes  was  selected.  To 
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demonstrate  the  ability  to  handle  arbitrary  excitation,  two  beams 
were  created  by  superimposing  on  the  original  excitation  the  same 
distribution  phased  to  produce  a  beam  at  45  degrees.  The  required 
generator  voltages  needed  to  create  the  current  excitations  for 
this  array  pattern  were  determined  using  the  methods  of  Section 
2.4,  with  ZT  =  50  +  jO  Q  at  each  generator.  The  results  are  shown 
in  Figs.  5-7.  Figures  5a  and  5b  compare  the  active  element 
impedances  obtained  with  current  generators  using  the  present 
method  (DFT)  and  using  a  standard  element-by-element  approach.  In 
Figs.  6a  and  6b,  the  active  element  impedances  obtained  with  both 
methods  using  voltage  excitation  are  compared.  In  both  cases  a 
generator  impedance  of  Z t  =  50  4-  jO  Q  was  used.  The  results 
compare  very  well  in  both  cases,  demonstrating  the  ability  of  the 
method  to  handle  both  voltage  and  current  excitations.  Finally, 
the  array  patterns  obtained  using  "uncorrected"  and  "corrected" 

voltage  excitations  are  shown  in  Fig.  7.  The  "uncorrected  feeds" 
simply  used  the  desired  current  taper  as  the  voltage  excitation 

taper.  This  method  gives  relatively  good  results  but  the  lobe  at 

45°  is  several  dB  below  the  ideal  pattern.  The  "corrected  feeds" 

were  obtained  by  determining  the  terminal  voltages  obtained  with  a 

current  generator  excitation  with  the  desired  taper.  These 

voltages  were  then  used  to  excited  the  array.  The  ideal  pattern, 
simply  the  array  factor  of  the  desired  element  currents  is  not 

shown,  but  it  is  indistinguishable  from  the  pattern  obtained  using 

the  "corrected  feeds". 
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5.0  CONCLUSION 

The  infinite  array  approach  for  analyzing  large  arrays  can  be 
adapted  to  the  efficient  analysis  of  finite  arrays  with  arbitrary 
excitation.  An  important  feature  of  the  method  is  that  at  least 
for  the  case  of  current  excitation,  the  storage  and  computational 
effort  numerical  effort  is  essentially  independent  of  the  number 
of  elements  in  the  array.  For  voltage  excitation,  the  method 
requires  only  the  storage  and  solution  of  an  array  of  order 
(N  xN  )2,  where  N  and  N  are  the  number  of  elements  in  the  array. 

x  y  x  y 

Fill  time  of  the  array  is  independent  of  the  number  of  elements. 

The  method  could  find  application  in  the  design  of  large 
phased  arrays  with  arbitrary  excitation,  where  the  active  element 
impedance  can  vary  significantly  over  the  array,  not  just  near  the 
edges.  Currently,  the  method  has  been  applied  with  good  results 
to  arrays  of  microstrip  dipoles.  The  method  works  for  both  free 
and  forced  excitation,  and  can  give  good  approximations  to  the 
mutual  impedance  between  array  elements. 
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Ground  Plane 


Figure  1.  A  section  of  a  periodic  array  of  microstrip  antennas 
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Figure  2.  Periodic  lattices  in  the  spatial  and  spectral  domains. 
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Figure  3.  Mutual  impedance  between  collinear  microstrip  dipoles.  Comparison 
of  results  using  DFT  method  and  results  in  [9].  a)  Dipoles 
collinear.  b)  Dipoles  parallel. 
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Figure  5.  Active  impedance  for  the  elements  in  a  19  element  array  current 

generator  excitation  selected  to  form  two  beams,  at  0  and  45 
degrees,  a)  Resistance,  b)  Reactance. 


Active  Imp.  with  Uncorrected  Voltage  Generators  Active  Imp.  with  Unconected  Voltage  Generators 
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Array  Pattern  — Corrected  Feeds 

Uncorrected  Feeds 

90* 


Figure  7.  Array  pattern  for  a  19  element  array,  with  voltage  generator 
excitation.  Taper  selected  to  form  two  beams,  at  0  and  45 
degrees.  Uncorrected  taper  is  identical  to  the  desired  currents, 

and  causes  errors  in  the  realized  pattern.  Corrected  taper  is 

the  excitation  which  gives  the  desired  element  currents. 

Generator  impedance  was  50  +  jO  £2. 
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MULTIMODE  PERFORMANCE  FROM  A  SINGLE  SLOTTED  ARRAY  ANTENNA 

by 

John  Cross,  Don  Collier  and  Len  Goldstone 
Norden  Systems,  Norwalk,  Ct. 

ABSTRACT: 

Today’s  military  aircraft  must  carry  out  various  missions,  some  of  them 
simultaneously.  Consequently,  antennas  must  be  designed  with  the  capability 
to  function  in  a  multiplicity  of  simultaneous  as  well  as  switchable  modes.  Slotted 
flat-plate  arrays,  with  their  unique  ability  to  divide  available  aperture  space 
among  several  antenna  requirements,  have  become  increasingly  popular. 
Several  recent  design  innovations  have  advanced  the  state  of  the  art  in 
attainable  operational  modes.  Taking  advantage  of  these  approaches,  one 
multimode  antenna  can  now  perform  air-to-air  pulse-doppler  monopulse, 
ground  mapping,  terrain  following,  terrain  clearance,  terrain  avoidance, 
weapons  transponder  tracking,  and  ground  target  tracking. 

Three  multimode  antenna  designs  will  be  presented,  ranging  from  the  simplest 
and  therefore  least  expensive,  to  the  most  complex,  namely  a  flat-plate  array 
design  containing  eight  separate  apertures  and  ten  RF  outputs.  Measured 
performance  data,  including  far-field  patterns  when  operating  in  the  various 
modes,  will  be  shown. 

1 .0  THE  LOW  COST,  MULTI-FUNCTION,  FLATPLATE  ARRAY  ANTENNA 
1.1  BACKGROUND:  The  Bandwidth  Tradeoff 

The  slots  on  any  flatplate  antenna  are  excited  by  waveguide  runs  located 
behind  the  front  plate  of  the  antenna  aperture.  RF  energy  travels  through  these 
waveguides  and  is  emitted  from  the  slots  whenever  a  component  of  H-field 
within  the  waveguide  is  parallel  to  the  slot  center  line.  The  slots  are  precisely 
spaced  to  insure  that  the  RF  from  each  slot  is  in  phase  with  that  of  its  neighbor 
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[1  ].  The  aperture  is  subdivided  into  resonant  subarrays  to  simplify  the  feed 
network.  This  simplification  limits  the  operational  frequency  band, 
however,  and  for  this  reason,  antennas  with  good  sidelobe  characteristics  at 
center  frequency  show  higher  sidelobes  and  VSWR  at  the  band  edges.  This 
deterioration  can  be  minimized  by  keeping  the  waveguide  runs  short.  As 
energy  travels  down  the  guide  from  the  RF  feed  port,  the  phase  error  of  the  RF 
increases  as  it  reaches  each  successive  slot.  Therefore,  if  the  antenna 
aperture  were  divided  into  small  subarrays,  and  each  subarray  were  fed  by  a 
separate  RF  feed  port,  the  phase  errors  could  be  kept  to  a  minimum.  However, 
this  method  is  expensive,  as  it  requires  a  complex  feed  network  attached  to 
the  back  of  the  antenna. 

Less  expensive  flatplate  antennas  must  therefore  be  designed  with  fewer 
subarrays  and  longer  resonant  waveguide  runs  in  each  subarray.  To 
compensate,  a  narrower  frequency  band  and/or  slightly  higher  sidelobe  levels 
are  tolerated. 

1.2  DESIRED  MODES  OF  OPERATION: 

Three  isic  modes  are  required  for  a  number  of  aircraft  radar  antenna 
applications:  a  pencil  beam  mode  for  long  range  target  location;  a  Cosecant- 
Squared  (CSC2 )  mode  for  mapping  and  ground  target  identification;  and  an 
elevation  monopulse  mode  for  Terrain  Following/Terrain  Avoidance  (TF/TA). 

1.3  DESIGN: 

Refer  to  Figure  1 .  In  this  "no-frills"  economy  design,  a  12X14  element  module 
was  chosen.  In  the  28  slot  elevation  plane,  there  was  sufficient  space  for  two 
modules,  and  in  the  azimuth  plane,  space  for  four.  This  limited  the  required 
number  of  RF  corporate  feed  ports  to  eight,  and  simplified  the  structure  of  the 
feed  manifold  on  the  back  of  the  array.  Equiphase  summing  the  eight  RF  ports 
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with  a  corporate  teed,  in  a  manner  similar  to  that  described  in  Reference  2, 
accomplishes  the  pencil  beam  performance. 

P.  Smith  of  Norden  Systems  developed  and  patented  a  means  to  create  a 
CSC2  elevation  beam  from  a  flatplate  aperture  by  1)  switching  the  power  to  only 
the  upper  half  the  radiating  aperture  and  2)  introducing  a  fixed  phase  shift 
into  the  central  row  of  the  active  aperture  half.[3]. 

Power  switching  is  accomplished  by  inclusion  of  an  RF  switch  in  the  feed 
manifold  as  shown  in  Figure  2.  In  more  complex  array  designs,  especially 
in  designs  requiring  low  rms  sidelobes  in  the  pencil  beam  mode,  the  central  row 
of  slots  has  its  own  feed  port,  with  a  phase  shift  section  that  can  be  switched  in 
or  out  as  required.  See  Figure  8,  for  example. 

In  the  simplified  design  (low  cost  array),  a  CSC2  pattern  is  obtained  by  placing  a 
fixed  phase  shift  in  the  central  row  of  aperture  slots.  See  Figure  3.  The 
resulting  deterioration  in  rms  sidelobe  level  in  the  pencil  beam  mode  was 
accepted  as  a  consequence.  See  Figure  4.  The  azimuth  beam  is  unaffected, 
as  can  be  observed  in  Figure  5. 

The  final  mode,  elevation  monopulse  for  TF/TA  operation,  becomes  possible 
with  the  inclusion  of  a  magic  tee  at  the  RF  input/output  port.  See  Figure  2.  The 
difference  pattern  is  taken  from  the  delta  port,  the  pencil  beam  pattern  from  the 
sum  port.  Figure  6  shows  the  resultant  sum  and  difference  patterns 
superimposed. 

2.0  A  MORE  COMPLEX  MULTIMODE  DESIGN  : 

The  layout  of  this  flatplate  antenna  is  shown  in  Figure  7.  The  design  contains 
not  only  an  Az/EI  monopulse  array,  switchable  to  form  a  ground  map 
antenna,  but  also  3  interferometers  and  a  guard  aperture.  The  guard  array 
serves  a  sidelobe  suppression  function  for  both  the  interferometers  and  the 
monopulse  array.  In  addition,  this  antenna  includes  two  BIT  (built  in  test) 


circuits  that  permit  a  reference  signal  to  be  injected  into  the  array  to  calibrate 
the  waveguide  runs  during  operation.  The  monopulse/ground  map  simplified 
diagram  is  shown  in  Figure  8.  The  ground  map  phase  shifters  and  RF  switch 
are  shown  as  well  as  a  T/R  circulator.  There  are  8  waveguide  ports  (1  guard, 
3  interferometer,  receive  sum,  azimuth  and  elevation  receive  difference,  and 
transmit  sum),  plus  2  coaxial  RF  injection  (BIT)  ports.  The 
interferometer/guard  BIT  diagram  is  shown  in  Figure  9. 

*  This  more  complex  antenna  is  subdivided  into  smaller  sized  modules  than  the 

simplified  array  previously  described,  to  achieve  a  wider  frequency 
bandpass  and  lower  pattern  sidelobes.  The  modules  are  basically  4X4  slot 
layouts.  The  array  diagram  illustrating  the  module  breakdown  is  shown  in 
Figure  7.  The  monopulse  array  has  60  modules.  Each  interferometer  has 
10  modules.  Each  subarray  or  module  is  fed  by  a  reactive  power  divider  that 
is  sensitive  to  module-to*module  impedance  mismatch.  Module  impedance 
thus  becomes  a  critical  design  parameter  if  the  desired  sidelobes  are  to  be 
achieved.  Any  mismatches  between  modules  lead  to  aperture  phase  and 
amplitude  errors  which  impact  on  the  desired  antenna  patterns.  Typical 
module  impedances  are  shown  in  Figure  10. 

The  complete  aperture  is  covered  by  a  layer  of  Kaptan  which  protects  the 
antenna  face  and  holds  waveguide  pressure.  The  ability  of  the  overall  design 
to  perform  within  specifications  depends  upon  manufacturing  tolerances, 
materials  properties,  and  the  quality  of  the  bonding  process.  The  array  is 
fabricated  largely  as  machined  and  dip  brazed  subassemblies  and  bonded  or 
bolted  in  various  manufacturing  stages  into  one  unit.  In-process  testing  and 
quality  control  are  vital  to  the  attainment  of  performance  goals.  Final 
adjustment  of  the  completed  unit  takes  place  on  the  antenna  range,  where  the 
8  monopulse  sections  are  trimmed  for  equal  phase,  using  shims  in  trombone 
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phase  shifters.  Pattern  performance  is  "fine  tuned"  by  adding  thin  layers  of 
dielectric  tape  over  the  slots  of  selected  modules,  when  required. 

Finally,  boresight  is  measured  and  the  array  mounting  pads  are  machined  to 
make  the  electrical  boresight  orthogonal  to  the  mounting  pads.  Typical  patterns 
are  shown  in  Figure  1 1 . 

3.0  A  COMPLEX  MULTIMODE  FLAT  PLATE  ARRAY: 

This  multimode  flat-plate  antenna  was  designed  for  a  radar  system  which  was 
to  be  integrated  with  a  fighter-bomber  aircraft.  The  projected  aperture  is 
approximated  by  a  truncated  ellipse  32"  wide  and  24"  high.  Contained  within 
this  area  are  5  separate  apertures;  an  azimuth  and  elevation  monopulse,  a  3- 
port  azimuth  interferometer,  a  2-port  elevation  interferometer,  a  single  port 
guard  and  a  single  port  weapons  uplink  for  a  total  of  10  RF  outputs.  All  of  the 
apertures  operate  over  a  3%  bandwidth  at  Ku  band  except  for  the  uplink,  which 
operates  at  X  band  and  is  used  in  transmit  only.  The  only  Ku  band  transmit 
channel  is  the  monopulse  sum,  which  has  the  capability  of  transmitting  either 
a  pencil  beam  or  shaped  beam  pattern  in  elevation  and  a  pencil  beam  in 
azimuth.  The  9  receive  channels,  namely  Sigma  (shaped  or  pencil  beam), 
Delta  Az  and  Delta  El  (for  monopulse),  Left,  Center  and  Right  (for  the  azimuth 
interferometer),  Up  and  Down  (for  the  elevation  interferometer)  and  Guard, 
are  connected  to  8  receivers  (Delta  Az  and  Delta  El  are  time  shared) 
mounted  on  the  back  of  the  array.  Each  microwave  receiver  consists  of  a 
transmit/receive  limiter  (TRL),  a  low  noise  amplifier  (LNA)  and  a  dual 
down-conversion  module  so  that  the  array  receive  outputs  are  at  the 
desired  IF  frequency.  Also  included  in  the  array  package  are  7  RF  injection 
ports  for  calibration  and  sensitivity  checks,  and  an  inertial  measurement  unit 
(IMU)  for  precise  antenna  orientation. 
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3.1  DESCRIPTION  OF  APERTURES  AND  FEED  NETWORKS: 

The  outline  of  the  monopulse,  azimuth  interferometer,  elevation 
interferometer,  guard  and  uplink  apertures  is  detailed  in  the  array  face-plate 
drawing  of  Figure  12.  The  individual  radiators  of  each  of  the  5  apertures  are 
horizontally  polarized  resonant  slots  which  are  grouped  into  subarrays  of 
resonant  waveguide  networks  and  fed  by  reactive  2-way  power  dividers.  It 
was  not  possible  to  use  4-port  power  dividers,  because  of  space  and  weight 
limitations. 

The  4-quadrant  monopulse  is  the  most  complicated  antenna  within  the  array, 
not  only  because  of  the  large  percentage  area  that  it  occupies,  but  also  due 
to  the  complexity  involved  in  the  feed  network  to  achieve  either  a  pencil  beam 
or  shaped  beam  in  elevation.  The  aperture  is  divided  into  4  main  quadrants 
with  16  subarrays  each,  combined  by  a  reactive  feed  network  and  a  comparator. 
The  center  rows  just  above  and  below  the  azimuth  center  line  have  a 
separate  feed  and  comparator.  The  switching  of  the  center  rows  provides 
elevation  beam  shape  diversity.  In  the  pencil  beam  mode,  the  entire  aperture  is 
illuminated  with  equal  phase  (Figure  14),  while  in  the  mapping  mode,  only  the 
lower  half  of  the  aperture  (PW  6  and  PW  7  of  Figure  15)  and  the  2  center 
rows  are  illuminated.  Optimization  of  the  shaped  beam  pattern  is  achieved  by 
advancing  the  phase  of  the  center  rows  relative  to  the  lower  quadrants. 
Transfer  switches  S2  and  S3  in  Figure  14  are  used  to  switch  in  the 
appropriate  phases  through  line  lengths  to  the  center  rows.  The  predicted 
shaped  beam  elevation  pattern  and  azimuth  sum  pattern  with  errors  are 
shown  in  Figures  16  and  17. 

The  feed  networks  for  the  remaining  apertures  are  similar  to  those  of  the 
monopulse  antenna,  but  have  no  comparator  assembly.  Resonant  waveguide 
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radiators  are  summed  through  a  teed  network  to  a  single  output  for  each  of 
these  arrays. 

3.2  MODES  OF  OPERATION: 

The  monopulse  aperture,  azimuth  and  elevation  interferometers  and  guard 
singly  or  in  combination  support  all  of  the  following  radar  modes.  Air-to- 
ground  modes  including  ground  mapping  (real  beam,  doppler  beam 
sharpening  and  synthetic  aperture),  beacon  and  GMTI  (pulse  doppler  and 
clutter  suppression  interferometry),  air-to-sea  modes  including  sea  surveillance 
and  ship  classification,  as  well  as  air-to-air  modes  including  air  surveillance 
and  air  combat,  terrain  clearance  including  terrain  following  and  terrain 
avoidance. 

3.3  RF  INJECTION: 

RF  injection  is  provided,  mainly  for  in-flight  calibration  of  the  antenna  system.  As 
shown  in  the  block  diagram  of  Figure  13,  monopulse  calibration  is  achieved 
by  single  quadrant  injection,  while  each  of  the  azimuth  and  elevation  channels 
as  well  as  the  guard  have  individual  injection  points.  In  addition,  each 
interferometer  injection  line  contains  a  variable  phase  trimmer  for  boresight 
alignment. 

4.0  CONCLUSIONS: 

Today's  complex  mission  requirements  make  flexibility  and  multimode 
capability  an  operational  imperative  in  modern  radar  systems.  To  support 
these  needs,  apertures  with  one  basic  antenna  pattern  are  no  longer  viable. 
Until  the  long-awaited  phased  array  becomes  economically  competitive,  the 
best  solution  to  date  has  been  the  multimode,  slotted  flat-plate  array  antenna. 
For  cost-effective  performance  in  a  variety  of  scenarios,  it  is  hard  to  beat. 
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Figure  1.  Front  Face  of  Low-Cost  Array  Antenna. 


Figure  2.  Simplified  Schematic  Showing  RF  Switch  in  CSC2  Mode  Position. 
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Figure  3.  CSC2  Antenna  Pattern  in  Elevation  Plane. 
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Figure  4.  Pencil  Beam  Pattern  in  Elevation  Plane. 
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Figure  5.  Azimuth  Beam  Pattern. 


6.  Elevation  Sum  and  Difference  Pattern 


Figure  7.  Front  Face  of  Array  Antenna. 


Figure  8.  Simplified  Diagram  of  Monopulse/Ground  Map  Array. 
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Figure  10.  Four  Typical  Module  Impedance  Curves. 
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Figure  1 1 .  Typical  Measured  Patterns. 
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Figure  12.  Complex  Flatplate  Array  Front  Face 


ANTENNA  FLAT  PLATE 


Figure  14. Monopulse  Waveguide 
Compatator  Schematic  Shown  in 
Monopulse  Sum  Mode. 
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Figure  15.Monopulse  Waveguide 
Comparator  Schematic  Shown  in 
CSC2  Mode. 
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Figure  17.  Simulated  Azimuth  Pencil  Beam  Pattern. 
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ABSTRACT 


This  paper  describes  a  new  adaptive  array  which  derives  a  reference  signal  from 
frequency  and  angle  of  arrival  information.  The  array  employs  a  Reference 
Acquisition  Loop  (RAL)  consisting  of  an  arbitrary  filter  and  a  nonlinear  directional 
gain  constraint  to  ensure  proper  operation.  The  array  is  more  flexible  since 
extremely  accurate  a  priori  information  is  not  required.  Simulation  results  are 
presented. 


1.  INTRODUCTION 

Adaptive  beamformers,  which  have  been  widely  discussed  and  applied  in  radar 
and  communication  systems1-4,  use  two  major  algorithms  (i.e.,  the  Least  Mean 
Square  (LMS)5  and  Applebaum  algorithms6).  For  proper  operation,  an  Applebaum 
array  needs  a  precise  steering  vector  to  point  the  beamformer  toward  the  desired 
signal  source;  the  LMS  algorithm  requires  a  reference  signal  which  is  highly 
correlated  with  the  information  signal  to  distinguish  the  unwanted  components  from 
the  desired  signals. 

To  circumvent  the  dependency  on  a  priori  knowledge,  several  studies  on  self- 
generated  reference  techniques  have  been  proposed.  These  studies  can  be  separated 
into  two  categories:  spread  spectrum  techniques7-10  and  'data-derived' 

techniques11.  For  both  techniques,  the  required  frequency  bandwidth  is  increased 
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and  an  accurate  desired  signal  characteristic  is  needed. 

We  propose  an  Autorecursive  Nonlinear  Adaptive  Array  (ARENA)  using  a 
Reference  Acquisition  Loop  (RAL)  to  extract  the  reference  signal.  This  array 
requires,  a  priori ,  the  frequency  and  Angle  Of  arrival  of  the  Desired  signal  (AOD). 
Thus,  the  proposed  array  is  a  more  general  system  for  most  applications. 

The  general  structure  and  the  basic  concept  of  the  array  are  described  in 
Section  2,  while  the  operation  is  explained  in  the  Section  3.  Computer  simulation 
results  are  presented  in  Section  4  and  the  Conclusion  is  in  Section  5. 

2.  FORMULATION 

Consider  an  m  element  half  wavelength  spaced  LMS  array  with  the  self- 
generated  reference  signal  as  shown  in  Fig.  1.  The  array  consists  of  two  parts:  the 
first  is  an  m  element  LMS  adaptive  processor  and  the  second  is  the  RAL.  Note 
that  the  reference  signal  is  derived  from  the  output  of  the  array.  The  RAL  must 
allow  the  desired  signal  component  at  the  output  of  the  array  to  pass  through  witn 
only  limited  distortion  (thus  retaining  correlation  between  the  desired  signal  and  the 
self-generated  reference  signal).  Also,  it  must  reduce  the  undesired  signal 
component  content. 

We  define  the  input  signal  vector,  X(k),  at  discrete  time  k,  to  be 


X(k)  =  Xd(k)  +  X,{k)  +  N(k ),  (1) 

where  Xd(k),  X{(k)  and  N(ife)  represent  the  input  desired  signal,  interference  and 
thermal  noise  vectors,  respectively.  Assume  Xd(k)  and  Xx(k)  are  not  correlated  to 
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each  other.  The  thermal  noise  vector  is  a  zero  mean  mutually  independent  random 
vector.  The  LMS  algorithm  is  used  to  adjust  the  weight  vector5,  W(A:), 


W(k  +  1)  =  W(k)  +  ±KVwE[\e(k)\2}  =  W(k)  +  KE{((k)X'(k))  ,  (2) 

where  k  is  the  gain  constant  of  the  adjustment  loop,  £■[•]  denotes  the  expected 
value  and  (•)  denotes  complex  conjugation.  e(k)  represents  the  system  residual 
signal  which  is  the  difference  between  the  self-generated  reference  signal  and  the 
system  output,  y(k),  i.e.,  €(fc)=G(y(fc)]-y (&),  where  <?[•]  denotes  the  function  of 
the  reference  signal  generator.  V  w  represents  the  gradient  with  respect  to  the 
weight  vector.  From  Eq.  (2),  when  the  weight  adjustment  loop  reaches  steady 
state,  t(k)  and  the  input  signal  are  orthogonal  to  each  other, 


£[«(*)X'(t)| 

=  WTWu)-X^k)W„)X‘(k)  +  (gi'(Uk)TW„)-xT(k)W„)x\k) 

+  (s,*(Af  T(k)W„)-NT[k)W„)X'(k)} 

=  0  ,  (3) 

where  Wt,  represents  the  steady  state  weight  vector,  9d{f),  y,(/)  and  gn(f)  are  the 
RAL  gains  of  the  desired  signal,  interference  and  the  thermal  noise,  respectively. 
The  superscript  (-)r  denotes  the  transposition.  For  simplicity,  we  assume  that  y^, 
gn  and  y,  are  constant  over  the  appropriate  frequency  band. 

When  both  the  desired  signal  and  interference  are  zero  mean  and  uncorrelated 
with  each  other,  Eq.  (3)  cam  be  further  decomposed  into 
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E\{gd-\)XTdWaax'd  +  (gi~l)X? WS,X‘  +  (ffn-l )NtN’}  =  0  . 


(4) 


To  ensure  that  the  steady  state  weight  vector  will  converge  and  suppress  the 
undesired  signal  component  the  RAL  is  designed  to  have  unity  gain  for  the  desired 
signal  but  non-unity  gain  for  the  interference  and  thermal  noise  components. 
Therefore,  Eq.  (4)  becomes 


[gi-l)E{X‘xJ'\W„  +  (g„-l)E[N‘NT\W..  =  0  .  (5) 


Obviously,  when  both  g,  and  gn  are  not  unity,  a  null  weight  vector  is  able  to  satisfy 
Eq.  (5)  which  is  a  trivial  solution.  Therefore,  a  particular  scheme  is  required  in  the 
RAL  to  avoid  this  while  minimizing  the  undesired  signal  power.  In  particular,  it  is 
desirous  to  have  gn,  y,<l  so  that  the  weight  adjusting  loop  will  have  a  faster 
convergent  characteristic  and  better  stability. 

3.  THE  ARENA  ARRAY 

A  type  of  adaptive  array  based  on  the  general  structure  shown  in  Fig.  1  and  a 
nonlinear  direction  constraint  in  the  RAL  is  proposed  as  shown  in  Fig.  2.  The  RAL 
consists  of  a  weight  constraint,  C(k),  and  a  filter,  F.  To  determine  the  design  of  the 
weight  constraint,  first  consider  an  output  signal  vector,  y{k), 


y (fc)  =  \Xd(k)+X,[k)+N{k))TW{k)  .  (6) 


Here 
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Xd{k)  =  ad(k)Sd  , 


(7) 


Xi(k)  =  £  a,(k)S ,  ,  (8) 

/= l 

and 

N(k)  =  [M*)>  n2(*)»  *m(*)]r  ,  (9) 

with  Sd  and  5/  (1=  l,2,...,p)  the  phase  arrival  vector  of  the  desired  and  interference 
signals,  respectively.  ad(k)  is  the  waveform  of  the  desired  narrowband  signal.  A 
continuous  (non-impulsive)  type  wide  band  interference,  Xt(k),  is  represented  by 
the  sum  of  narrow  band  waveforms  aj(&),  /=  1,  ...,  p.  n,(£),  »=1,  2,  •  •  *  ,  m 
represents  the  thermal  noise  component  at  i-th  element.  Substituting  Eqs.  (7)-(9) 
into  Eq.  (6) ,  the  array  output  is  given  by 

y(k)  =  a d(k)sf  tr(t)+£  o,(k)S?  W{k)+NT(k)  W(k)  .  (10) 

.  1  =  1 

We  define  the  output  of  the  RAL,  r(k),  to  be 

'(*)  =  <?(*)%(*)! . 

=  C(k)[gjad{k)Sj W(k)+£  „ a,(k))Sf  W(k)+!„NT(k)  1 V(*)|  ,  (11) 

1=1 
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where  F[-\  is  a  filter  function,  C(k)  is  a  directional  constraint,  and  gt  is  the  gain  seen 
by  a\.  Note  that  the  output  of  the  RAL  is  utilized  as  the  reference  signal  and  fed 
back  to  the  weight  adjusting  loop.  We  define  the  system  residue,  e(fc),  as 

<(*)  =  r(k)  -  y(k) 

/  =  1 

(12) 


In  order  to  maintain  a  constant  gain  in  the  designated  direction  and  also  avoid 
the  zero  weight  convergence,  we  have 


C(k)  = 


9 

gtWT(k)St  ’ 


(13) 


where  g  is  a  arbitrary  gain  constant.  Note  that  C(k)  is  a  nonlinear  constraint  and  is 
inversely  proportional  to  the  gain  of  the  weight  vector  in  the  AOD.  Substituting 
Eq.  (13)  into  Eq.  (12),  the  system  residue  now  becomes 


(14) 
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3.1  Optimum  Weight  Vector 


According  to  the  complex  LMS  algorithm12,  the  weight  vector  which  yields  the 
minimum  E\  |  e(k)  |2]  can  be  obtained  by  solving 


E[e(k)V  we(k)‘]  =  0  . 


(15) 


After  some  manipulations,  Eq.  (15)  can  be  rewritten  as 


E\t(k)Vw/(k)\ 


WW~  W'W  WWS'J~ "'Wl 

9dSdwik)  9d{S!iW{k)Y 

=  0  ,  (16) 


where  (•)*  denotes  complex  conjugate  transpose.  Note  that  the  first,  second,  and 
third  terms  in  Eq.(l6)  correspond  to  the  desired  signal,  interference,  and  noise 
contributions,  respectively. 

It  can  be  shown  that  the  solution  for  the  weight  vector  of  Eq.  (16)  when 
9i=9n= 0,  for  all  /,  and  0  is  the  optimal  Wiener  solution,  i.e.  W  =  $Zl$d  where 
$u  is  the  correlation  matrix  of  the  undesired  signals.  Otherwise,  W(fc)  will  be  a 
sub-optimal  weight  vector.  However,  both  optimal  and  sub-optimal  weight  vectors 
tend  to  annihilate  the  undesired  signal. 
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From  the  first  term  in  Eq.  (16)  we  see  that  S  J  W(k)=g.  This  condition,  drives 
the  weight  to  form  a  gain  g  at  the  AOD.  Combining  this  requirement  with  the 
constraint  of  Eq.  (13),  results  in  maintaining  unity  gain  for  the  desired  signal  in  the 
RAL.  Consequently,  this  structure  offers  a  more  flexible  design  of  the  filter  in  the 
RAL  and  better  system  stability. 

3.2  System  Realization 

To  realize  the  ARENA  structure,  the  gradient  of  the  residual  signal  versus  the 
weight  vector  is  required.  However,  this  is  not  generally  available  in  practice. 
According  to  the  orthogonality  principle13,  in  such  situations,  V  (ye(A:)  is  replaced  by 
X(k).  The  weight  adjusting  equation,  then,  is  written  as 


W(k+ 1)  =  W(k)  +  «[*(*)*'(*)]  ,  (17) 

where  k  is  a  step  size  constant.  Comparing  Eqs.  (17)  and  (2)  we  see  that  the 
weight  adjusting  equation  of  the  ARENA  array  is  equivalent  to  the  conventional 
LMS  array.  Thus,  the  steady  state  weight  vector  can  be  solved  in  the  same 
manner,  i.e. 


1=1 


Qd 


Sj  W{k) 
sjw{k) 


Si  +<7n( 


9n9 


9dSTdW{k) 


-!)*(*) 


=  0  . 


(18) 
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From  Eq.  (18)  we  find  that  this  practical  weight  adjusting  algorithm  has  a 
sii.  ilar  solution  as  Eq.  (16)  which  minimizes  the  residual  signal  power  by  letting  the 
weight  vector  approach  SjW(k)=g  and  Sf  W (k)=Q.  Note  we  must  ensure  that 
gi^9d  and  gn^Qi  since  if  gi—gd  it  means  the  RAL  has  completely  lost  its  ability  to 
distinguish  the  interference  and  the  desired  signal  or  the  interference  has  the  same 
characteristic  as  the  desired  signal,  the  interference  will  be  treated  as  the  desired 
signal  or  vice  versa.  A  matched  filter  can  be  used  to  accomplish  this.  For  a  wide 
band  noise,  a  white  noise  cancelor14  cam  be  utilized  in  the  RAL  to  decrease  the 
RAL  gain  for  noise  and,  consequently,  avoid  gn—gd • 

4.  SIMULATION  RESULTS 

In  this  section  we  look  at  several  examples  which  illustrate  the  results  derived 
above.  We  first  examine  a  case  v.\h  a  desired  signal  and  one  narrow  band 
interference  signal.  We  use  a  5  element  half  wavelength  spaced  ARENA  array. 
The  frequency  of  the  desired  signal  is  fd  and  the  interference  is  1.258 /j.  The  input 
Signal-to-Noise  Ratio  (SNR)  and  Interference-to-Noise  Ratio  (INR)  are  10  dB  and 
20  dB,  the  AOD  and  the  Angle  Of  arrival  of  the  Interference  (AOI)  are  15*  and 
30  ° ,  respectively.  We  assume  the  gain  of  the  RAL  for  the  desired  signal  is 
04=0  dB.  We  consider  two  cases  for  the  gain  of  the  interferer,  i.e.  noise  gain  is 
assumed  to  be  gn=—\dB.  In  the  computer  simulations,  the  initial  weight  values  aure 
all  set  to  zero  except  for  the  first  element  which  is  arbitrarily  set  to  0.01.  The 
output  SENR  of  the  aura y  is  illustrated  in  Fig.  3  for  both  cases.  We  see  that  the 
weights  converge  as  expected  and  drive  the  system  to  the  optimal  output 
performance;  however,  the  case  of  weaker  suppression  of  the  undesired  signal  by 
the  RAL  filter  results  in  the  poorer  output  SENR  performamce  as  indicated  by  the 
slower  convergence  and  greater  jitter. 
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The  major  innovation  of  the  proposed  array  is  the  employment  of  a  Nonlinear 
Direction  Constrainer  (NDC)  in  the  RAL.  The  NDC  prevents  the  weight  vector 
from  converging  to  zero  by  constraining  a  constant  gain  in  AOD.  You  may  wonder 
if  the  AOD  is  known  why  we  need  the  reference  signal?  With  a  precise  estimation 
of  the  AOD  an  Applebaum  array6  can  be  used  to  achieve  satisfactory  performance. 
Previous  study15  has  shown,  however,  that  the  Applebaum  array  is  extremely 
sensitive  to  look  direction  error.  The  nonlinear  direction  constraint  used  in  the 
RAL  is  a  soft  constraint  which  allows  for  some  error  without  causing  significant 
deterioration  on  the  system  performance.  Because  of  this  low  sensitivity  to  steering 
vector  error,  the  performance  of  the  ARENA  array  is  superior  to  the  Applebaum 
array  when  a  precise  steering  vector  is  not  available.  Simulations  have  been 
•'onducted  to  verify  the  results.  Consider  a  5  element  ARENA  array  with  half 
wavelength  spacing  at  the  center  frequency,  fa,  of  the  desired  signal.  The  AOD  is 
15  °  and  AOI  is  30  ° .  The  frequency  of  the  interferer  is  1.258/,*.  The  loop  gains  of 
the  RAL  are  y<*= 0  dB,  y*=— 20  dB  and  y„=— 1  dB.  The  input  SNR  and  INR  are  10 
dB  and  20  dB,  respectively.  The  output  SINRs  of  the  ARENA  and  Applebaum 
arrays  with  an  estimated  AOD  equal  of  25  *  are  illustrated  in  Fig.  5.  The 
simulation  results  indicate  that  the  ARENA  array  is  much  less  sensitive  to  the 
estimated  error  than  the  Applebaum  aoray. 

5.  CONCLUSION 

In  this  study  we  propose  an  adaptive  array  which  employs  a  direction 

constrained  Reference  Acquisition  Loop  (RAL)  to  extract  the  reference  signal  from 

the  array  output  without  using  spread  spectrum  or  coding  techniques.  The  a  priori 

knowledge  required  for  the  proposed  array  is  the  AOD  and  the  frequency  band  of 

the  desired  signal.  It  has  been  shown  that  the  proposed  array  is  less  sensitive  to  the 
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estimation  error  of  the  AOD.  The  theoretical  analysis  and  the  computer  results 
indicate  that  the  proposed  array  successfully  nulls  the  unwanted  signal  and 
maintains  a  constant  gain  in  the  arrival  angle  of  the  desired  signal. 
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Fig.  1  General  Struture  of  an  LMS  array  with 
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Fig.  5  The  output  SINK  of  the  ARENA  with  gd=0  dB,  g,=-20dB, 
gn=-ldB,  AOD=15  *,  AOI=30*,  and  estimated  AOD=25  * 
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1.0  ABSTRACT 


This  report  covers  the  work  done  under  Task  5  of 
Project  45194263,  entitled,  "Communications  Adaptive 
Array  Processor  Evaluation". 

Adaptive  antenna  technology  has  been  an  area  of 
interest  for  many  years.  The  progress  to  date  can  be 
characterized  by  extensive  research  on  algorithms  and 
processing  architectures,  development  of  a  large  number 
of  prototype  systems  and  fielding  of  a  few  systems. 

Adaptive  antenna  systems  normally  consist  of  two  or 
more  antenna  elements,  two  or  more  receivers,  amplitude 
and  phase  weighting  m  vorks  for  one  or  more  elements 
and  a  signal  combiner.  For  some  applications  the  size, 
weight,  power  consumption  and  cost  of  multiple 
receivers,  antenna  elements  and  associated  radio 
frequency  (RF)  cabling  can  be  limiting  factors.  In  spite 
of  these  limitations,  the  impressive  interference 
cancellation  capabilities  of  adaptive  antenna  systems  is 
highly  desirable.  Recent  advances  in  photonics 
technology  could  provide  solutions  to  the  limitations 
described  above.  RF  modulated  lasers,  optical  cabling 
and  optical  detectors  could  reduce  much  of  the  size, 
weight,  power  consumption  and  cost  of  conventional 
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adaptive  array  systems. 

Using  the  RADC/DC  RF  laboratory  anechoic  chamber  and 
the  "Flexible  Adaptive  Spatial  Signal  Processor"  (FASSP) 
this  effort  was  performed  to  established  the  feasibility 
of  using  optical  technology  to  implement  an  adaptive 
antenna  array  system.  This  report  covers  the  design, 
implementation  and  testing  of  an  optically  linked 
adaptive  antenna  array.  The  performance  of  the 
optically  linked  array  is  evaluated  by  comparison  to  a 
conventional  adaptive  array. 

2.0  IN-HOUSE  TEST  FACILITIES 

RADC/DC  has  an  adaptive  array  processing  test  bed. 
The  test  bed  shown  in  Figure  1  consists  of  a  rectangular 
anechoic  chamber,  a  flexible  adaptive  spatial  signal 
processor  (FASSP) ,  a  antenna  pattern  recorder,  various 
types  of  jammer/desired  signal  sources  and  satellite 
communication  simulation  and  analysis  programs. 

The  test  bed  simulation/analysis  computer  programs 
are  used  to  study  and  compare  adaptive  processing  system 
concepts,  techniques  and  algorithms.  This  provides  a 
fast  look  approach  to  determine  the  merit  and  feasibility 
of  a  concept.  If  the  results  show  promise  the  concept 
is  further  tested  using  real  signals  and  adaptive 
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Figure  1.  ADAPTIVE  ARRAY  PROCESSING  TESTBED 

processor  hardware  to  determine  the  actual  benefit 
attainable.  The  test  bed  is  reconf igurable  and  functions 
as  a  tool  to  support  the  development  of  methodologies 
for  comparing  and  evaluating  new  adaptive  array 
processing  algorithms,  architectures  and  techniques 
suitable  for  meeting  satellite  communications 
requirements  (such  as  those  for  the  Defense  Satellite 
Communications  System  [DSCS]). 

2.1  ANECHOIC  CHAMBER 

The  anechoic  test  chamber  shown  in  Figure  2  is  a 
rectangular  structure  40  ft.  long,  28  ft.  wide  and  18 
ft.  high2.  The  inner  chamber  is  isolated  from  RF  fields 
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Figure  2.  anechoic  chamber 

from  150  MHz  to  18  GHz  by  at  least  100  dB. 

It  has  a  six  foot  diameter  spherical  quiet  zone 
located  at  its  longitudinal  axis.  The  center  of  the 
quiet  zone  is  midway  between  the  ceiling  and  the  floor 
and  about  50  inches  from  the  tips  of  the  absorber  on  the 
back  wall.  The  receive  element  array  is  positioned  in 
the  center  of  the  quite  zone  to  minimize  reception  of 
all  reflected  signals.  All  chamber  walls,  ceilings  and 
floors,  except  walkways,  are  covered  completely  with 
energy  (RF)  absorbing  material.  A  Scientific  Atlanta 
model  5315C-5  antenna  positioner  is  installed  in  the 
chamber.  The  tip  of  the  model  tower  (which  supports  the 


array  elements)  is  located  in  the  center  of  the  "quiet 
zone".  The  chamber  is  wide  enough  and  has  provisions 
so  that  several  signal  sources  can  be  used 
simultaneously  at  the  front  wall  opposite  to  the  "quiet 
zone" . 

Six  feet  of  the  40  ft.  chamber  is  partitioned  off 
and  is  used  as  an  equipment  room  to  house  the  signal 
sources  and  antenna  positioner  controls.  Absorber 
panels  are  removable  to  allow  access  for  mounting 
signal/ jammer  antennas.  Signal  and  control  connections 
between  the  chamber  and  the  laboratory  equipment  (FASSP 
and  Scientific  Atlanta  2020  system)  are  provided  through 
bulkhead  feed  through  panels  at  each  end  of  the  chamber. 
All  functions  such  as  source  power,  frequency,  mode  and 
receiver  antenna  position  (pedestal  rotation)  are 
controlled  from  outside  of  the  chamber. 

2.2  FLEXIBLE  ADAPTIVE  SPATIAL  8IGNAL  PROCESSOR  (FASSP) 

The  FASSP4  test  bed  is  shown  in  Figure  3.  All 
adaptive  Spatial  processing  systems  consist  of  three 
generic  components: 

1.  An  array  of  receiving  elements  (spatial  array) 
to  provide  the  degrees  of  freedom  required  to  null  out 
a  number  of  directional  jamming  signals, 


280 


Figure  3.  FLEXIBLE  ADAPTIVE  SPATIAL  SIGNAL  PROCESSOR 

2 .  An  adaptive  processor  that  uses  the  signal 
samples  from  the  array  receive  sensors  to  compute  the 
adaptive  weights  that  produce  the  resultant  spatial 
response,  and 

3.  Weighting  networks  that  apply  the  adaptive 
weights  to  the  signals  from  the  individual  input 
channels . 

The  design  of  adaptive  spatial  processing  systems 
is  very  complicated  because  of  the  close  interaction 
among  these  three  basic  components.  Although  computer 
simulations  can  be  used  to  compare  the  performance  of 
adaptive  algorithms  and  techniques,  the  hardware 


implementation  effects  cannot  always  be  easily  modeled. 

With  this  in  mind,  RADC/DCCD  conceived  the  ideas  of 
a  Flexible  Adaptive  Spatial  Signal  Processor  (FASSP) , 
which  was  designed  and  fabricated  for  RADC  by  Syracuse 
Research  Corp.  The  FASSP  is  a  general  purpose 
flexible  hardware  adaptive  array  processor  system  that 
supports  the  integration/ test  of  adaptive  processing 
algorithms,  architectures,  techniques  and  real 
components . 

The  FASSP  system  was  fabricated  with  high 
performance  quality  components  and  consists  of  12  real 
RF  receivers  and  weighting  networks  that  are 
reconf igurable.  The  designer  uses  a  computer  based 
operating  system  to  select  the  adaptive  algorithms  and 
hardware  configuration  and  specify  the  necessary  system 
parameters.  Adaptive  processor  performance  can  then  be 
evaluated  against  RF  signal  (desired)  and  jamming 
signals  using  an  antenna  array  and  an  anechoic  chamber. 
3.0  CONVENTIONAL  ADAPTIVE  ARRAY 

The  conventional  adaptive  array  processing  system 
shown  in  Figure  4  was  used  as  the  baseline  system 
configuration.  The  receiver  antenna  was  a  line  array  of 
four  standard  gain  horns  separated  by  3/2  wavelength 
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Pigure  4.  CONVENTIONAL  ADAPTIVE  ARRAY 

spacing.  The  array  was  located  on  the  positioner  with 
the  elements  centered  in  the  chamber  quiet  zone.  The 
optically  linked  adaptive  array  configuration  is  shown 
in  Figure  5  and  is  described  in  Section  4. 

The  CW  and  wide  band  noise  jammer  antennas  were 
located  at  the  opposite  end  of  the  chamber.  One  jammer 
transmitter  element  was  positioned  broadside/boresight 
(zero  degrees  azimuth,  zero  degrees  elevation)  to  the 
array  and  a  CW  signal  was  applied.  Figure  6  is  a  plot 
of  the  array  steered  to  that  CW  source. 

The  array  response  plot  shows  that  the  main  beam  is 
centered  at  zero  degrees  azimuth  and  the  sidelobes  are 
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Figure  5.  OPTICALLY  LINKED  ADAPTIVE  ARRAY  CONFIGURATION 

located  at  plus  or  minus  fourteen  degrees  azimuth. 
The  second  jammer  transmit  element  was  positioned  at 
fourteen  degrees  centered  on  one  of  the  side  lobes. 

Nulling  tests  were  then  performed  using  single 
jammers  at  each  position. 

3 . 1  BASELINE  PERFORMANCE 

The  performance  measure  used  to  compare  the 
performance  of  the  optical  implementation  with  that  of 
the  standard  coaxial  implementation  was  cancellation 
ration  (CR)3. 

CR's  for  a  given  input  jamming  signal  are  measured 
by  comparing  the  adaptive  and  conventional  output  power 
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Figure  6.  BEAMFORMED  ARRAY  PATTERN 

of  the  adaptive  processor.  Usually,  CR  is  used  to 
characterize  the  performance  of  the  adaptive  processor 
against  different  jamming  scenarios.  It  is  a  direct 
measure  of  the  adaptive  processor  nulling  capability 
and  is  shown  in  Equation  (1)  where  Pa  is  the 

CR(dB)  =  101og(Pa/Pc)  (1) 

adaptive  processor  output  power  in  the  adaptive  mode  and 
f  Pc  is  the  adaptive  processor  output  power  in  the 

conventional  or  beamformed  mode.  Normally  adaptive 
processor  output  is  measured  in  dbm  and,  thus,  Equation 
(2)  is  used  to  compute  the  CR  ratio. 

CR(dB)  =*  Pa(dB)  -  Pc(dB)  (2) 
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4.0  THE  OPTICALLY  LINKED  ARRAY 


An  adaptive  array  usually  consists  of  a  number  of 
antenna  elements,  a  receiver  for  each  element  and  the 
RF/IF  cabling  and  connectors  necessary  to  connect  the 
system  together.  One  channel  of  such  a  system  is  shown 
in  Figure  7.  The  block  marked  “Flexible  Adaptive  Spatial 
Signal  Processor”  (FASSP)  indicates  the  equipment  that 
performs  the  adaptive  signal  processing  tasks.  Because 
of  the  losses  in  the  RF  cabling,  the  SHF  receivers  are 
co- located  with  the  antenna  elements  in  the  anechoic 
chamber.  Along  with  the  RF/IF  cabling,  the  receivers 


Figure  7.  CONVENTIONAL  ADAPTIVE  ARRAY  CONFIGURATION 

introduce  another  source  of  unwanted  RF/IF  radiation. 
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Shown  in  Figure  8  is  an  optically  linked/connected 
adaptive  antenna  array. 


Figure  8.  OPTICALLY  LINKED  ADAPTIVE  ARRAY 

Because  of  the  very  low  losses  in  the  optical  fiber 
it  is  possible  to  remove  the  SHF  receivers  from  the 
chamber  and  co-locate  them  with  the  FASSP.  The  optical 
linkage  also  eliminates  the  IF  cabling  from  t’.e  chamber. 

4 . 1  HARDWARE 

Component  selection  was  based  on  anechoic 
chamber  and  associated  equipment  interface  requirements. 
Since  the  antenna  system  was  designed  to  run  at  8.41  GHz, 
the  task  of  optical  link  selection  was  significantly 
simplified.  The  Ortel  TSL-1000  optical  transmitter  and 
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the  Ortel  RSL-25  optical  receiver  were  selected  based 
on  their  frequency  range  of  10  kHz  to  10  GHz. 

The  maximum  received  signal  likely  to  be  seen  at  the 
output  of  each  antenna  element  (feedhorn)  is  less  than 
-40  dBm.  Therefore,  50  dBm  gain,  narrow  band  (centered 
at  8  GHz)  Miteq  amplifiers  were  selected  to  boost  the 
signal.  This  will  supply  the  Ortel  transmitters  with 
power  near  their  maximum  RF  input  level  of  12  dBm  and 
thereby  allowing  us  to  maximize  the  usable  dynamic  range 
of  the  links. 

Other  devices  are  incorporated  in  each  elemental  leg 
of  the  antenna  system  (see  Figure  8)  .  These  include 
electrical  and  optical  isolators  to  prevent  impedance 
mismatch  reflections  from  the  feedhorn,  and  optical 
reflections  to  the  laser,  respectively,  50  meter  lengths 
of  fiber  and  phase  adjusting  mechanisms  to  time  match 
the  links  are  also  integral  parts  of  the  system1. 

4 . 2  FUNCTIONAL  DESCRIPTION 

Each  array  link  consists  of  the  equipment  described 
in  Section  4.1.  An  adaptive  array  is  composed  of  two 
or  more  links,  usually  one  for  each  antenna  element  in 
the  array.  In  terms  of  a  single  link  the  RF  signal 
arrives  at  the  array  and  is  received  by  the  antenna 


element,  in  this  case  an  SHF  feedhorn.  The  low  level 
signals  are  amplified  to  provide  an  appropriate  signal 
level  for  the  optical  transmitters.  The  optical 
transmitters  put  out  a  single  frequency  of  light  the 
amplitude  of  which  is  modulated  by  the  RF  input  signal. 
The  light  is  then  transmitted  along  an  optical  cable  to 
the  phase  adjusting  device  where  the  electrical  length 
of  the  optical  cable  can  be  adjusted  so  that  all  the 
links  in  the  array  can  be  time  and  phase  matched5.  The 
modulated  and  time  and  phase  matched  optical  signal 
then  travels  along  the  optical  cable  to  the  optical 
detectors  where  the  light  is  demodulated  and  the  RF 
signal  is  output  to  the  SHF  receivers  through 
conventional  RF  cabling. 

5.2  TEST  AND  EVALUATION 

The  optically  linked  array  was  tested  in  a  manner 
identical  to  the  conventional/baseline  array  as 
described  in  Section  3.0.  Data  was  collected  to  plot 
cancellation  ratio  for  three  separate  sensors3.  All  of 
the  figures  indicated  below  are  identical  in  form.  The 
horizontal  axis  (Pe)  is  the  ratio  of  signal  plus  noise 
to  noise  power  as  measured  at  the  antenna  element.  On 
the  vertical  axis  two  quantities  are  plotted.  The 
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quantity  (Pc)  is  referred  to  as  the 
conventional /unadapted  output  power  as  measured  at  the 
beamformer/residue  output  port  of  the  adaptive 
processor.  When  measuring  this  quantity  the  array  is 
beamformed  to  broadside  and  then  the  interference  signal 
is  activated  and  again  the  signal  plus  noise  to  noise 
ratio  is  measured.  The  quantity  (Pa)  is  referred  to 
as  output  power  after  the  adaptive  processor  is  allowed 
to  adapt.  Again  the  signal  plus  noise  to  noise  ratio 
is  measured  at  the  residue  output.  Using  these  graphs 
the  CR  in  dB  can  be  read  directly  as  the  difference 
between  the  two  plotted  lines  on  the  graphs. 

Shown  in  Figure  9a  is  a  plot  of  the  conventional  and 
adaptive  response.  Figure  9b  is  a  plot  of  the  optically 
linked  conventional  and  adaptive  response.  Both 
Figures  9a  and  9b  are  for  a  single  continuous  wave  (CW) 
signal  broadside  to  the  array. 

Shown  in  Figure  10a  is  a  plot  of  the  conventional  and 
adaptive  response.  Figure  10b  is  a  plot  of  the 
optically  linked  conventional  and  adaptive  response. 
Both  Figures  10a  and  10b  are  for  a  single  CW  source  in 
the  first  sidelobe  which  is  14  degrees  from  broadside 
for  the  antenna  array  in  use. 


Shown  in  Figure  11a  is  a  plot  of  the  conventional 
and  adaptive  response.  Figure  lib  a  plot  of  the 
optically  linked  conventional  and  adaptive  response. 
Both  Figures  11a  and  lib  are  for  a  wide  band  noise  (WBN) 
source  broadside  to  the  array. 

6.0  CONCLUSIONS  AND  RECOMMENDATIONS 

Analysis  of  the  data  presented  in  Section  5.0 
indicates  that  the  performance  of  the  optical  linked 
array  is  as  good  as  the  baseline  conventional  array  and 
is  in  some  cases  better.  During  the  design  and 
fabrication  of  the  optically  linked  array  it  was 
apparent  that  optical  technology  would  provide  reduced 
size  and  weight  over  the  conventionally  linked  array. 
The  RF  cabling  associated  with  a  conventional  link  is 
bulky  and  weighs  about  15  -  20  lbs.  as  compared  to  the 
optical  link  which  is  far  less  bulky  and  weighs  only  a 
few  pounds,  including  modulator  and  demodulator  These 
size  and  weight  improvements  are  obtained  with  no  loss 
in  performance  and,  in  some  cases,  with  even  better 
performance. 

Further  improvements  are  attainable  by  the  following 
means.  An  increase  in  the  dynamic  range  of  the  optical 
links  could  be  obtained  by  using  indirect  modulation  of 
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Figure  10a.  CONVENTIONAL  CW  SIDELOBE  CR 


•  Pc  -t-  pa 


Figure  10b.  OPTICAL  CW  SIDELOBE  CR 
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Figure  11a.  CONVENTIONAL  WBN  BROADSIDE  CR 


Figure  lib.  OPTICAL  WBN  BROADSIDE  CR 
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the  lasers  verses  the  direct  modulation  used  for  these 
tests.  Single  mode  fiber  optic  cable  should  be  used 
instead  of  multimode  to  prevent  cable  losses  due  to 
gradient  effects  when  the  fiber  is  bent  or  deformed  in 
any  way.  Using  a  single  laser  with  external  modulators 
for  each  antenna  element  would  ensure  system  coherence 
and  make  the  system  even  more  cost  effective.  The  single 
laser  technique  would  also  allow  all  the  elements  to  be 
multiplexed  on  to  a  single  fiber  that  could  then  be 
demultiplexed  at  the  adaptive  processor.  The 
detector (s)  could  be  RF  biased  to  provide  down 
conversion,  thereby  eliminating  the  FASSP  receivers. 
Elimination  of  the  FASSP  receivers  will  result  in  a 
reduction  in  size,  weight  and  power  of  a  complete 
adaptive  antenna  system.  These  means  of  improvement 
will  be  further  investigated  and  verified  in  future 
in-house  efforts.  This  RADC/DCCD  in-house  initial 
investigation  of  optical  technology,  as  applied  to 
adaptive  antenna  systems,  indicates  that  evolving 
optical  technology  can  and  will  provide  important 
enhancements  especially  for  mobile  platform  systems. 
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